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FOCUSING OF RADIO WAVES REFLECTED FROM A ROUGH 
CURVED IONOSPHERE 


By J. D. WuitEHEAD* 
[Manuscript received June 8, 1960] 


Summary 


The reflection of radio waves from a partially rough, curved ionosphere is considered. 
The relationship between the amplitude of the echo, A, and phase path P when the 
ionosphere moves overhead with a horizontal velocity V at a height h is the same as 
that for a smooth curved ionosphere, i.e. 
h @P 
A?a1]1— —~ —_, 
2V2 df 


although because of the different physical conditions the large increases in echo amplitude 
observed when reflection takes place from a smooth ionosphere are not expected for 
reflection from a rough ionosphere. 

A method of testing the relationship experimentally is suggested. 


I. INTRODUCTION 

It occasionally happens that the ionosphere behaves as a specular reflector 
of radio waves incident normally on it, so that for many purposes it may be 
treated as a smooth mirror. However, Munro (1950) and Bramley and Ross 
(1951) have shown that this smooth “‘ mirror ”’ is not plain but has ripples moving 
through it; the ripples having a wavelength of about 100 km and an amplitude 
of a few kilometres. As they move overhead, the ripples have the effect of 
producing large increases in the amplitude of the reflected radio waves: this 
effect is known as focusing and the author (Whitehead 1956) has shown that, 
for vertical incidence, the amplitude A and phase path P are related by the 


equation 
eee, Sees 
2 eoy2 art 


when h is the average height of the reflector and V is its horizontal velocity. 


However, it is more usual for the ionosphere to be a rough reflector, so that 
the echo fades and the amplitude has a Rayleigh distribution (Mitra 1949). It 
is of interest to enquire whether ‘“ focusing” is to be expected under these 
circumstances if the rough reflector is curved. It is indeed found that the mean 
- amplitude fluctuates slowly (Meadows and Moorat 1957), and it is of interest to 

see if there is a similar relationship to (1) between the mean amplitude and the 
phase path. It is the purpose of the present paper to derive this relationship. 
However, it is shown that the large increases in mean amplitude are not to be 
expected because of the difference in the physical situations. 


* Radio Research Board Laboratory, Electrical Engineering Department, University of 
Sydney. 
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IJ. DERIVATION OF THE EQUATION 

Briggs and Phillips (1950) have shown that, for radio wavelengths shorter 
than about 200m, the ionosphere can usually be treated as a partially rough 
reflector so that most of the power returned to a transmitter on the ground is 
contained within a cone of semi-angle 0, if the ionosphere is flat. 9, was found 
to be about 5°, so that the effective area of the ionosphere for reflection is that 
contained within a circle of radius 9h, about 10-30 km, small compared with 
the wavelength of the ripples. 

Now suppose the rough reflector is slightly curved concave downwards 
with a radius of curvature R>h, the curvature being in one dimension only so 
that it forms part of a cylinder. Power is now returned from a greater area of 
the ionosphere, an ellipse with semimajor axes 0, and oh, where ¢ is the angle 
an incident ray makes with a line drawn perpendicular to the reflecting surface 
from the transmitter when the specularly reflected ray makes an angle 20, to 
the incident ray in the plane of curvature of the ionosphere (Fig. 1). Thus the 
power returned is increased in proportion to 9/0). It is seen that 


sing sin 0) 


Rere—h 
therefore 
Bing h 
sin 0, Boh ' RW 


so that the square of the amplitude is given approximately by 
h 
a w1+——. 
op/Oy~1+ Wick 


This is identical with the relation derived for the amplitude of a wave 
reflected from a curved mirror: thus it follows that when the rough curved 
reflector moves with velocity V at right angles to the axis of the cylinder, 
“hap 
2Vs. de 


In this equation A? and d?P/dt? are to be interpreted as mean values taken 
over several rapidly fading cycles. 


tA oe 


IIT. Discussion 

Although the same equation has been derived for both rough and smooth 
ionospheres, there is an important physical difference. For a mirror reflection 
the principal contribution to the amplitude comes from the first few Heetie: 
zones: the amplitude is proportional to the area of the zones, or for cylindrical 
curvatures, to the length of the (elongated) zones. However, for a rough reflector 
it is the power or the square of the amplitude which is proportional to the area 
or length, of the effective reflecting part of the ionosphere. With a smooth flat 
reflector at a height of 200 km, the first Fresnel zone has a radius of about 5 km 
at a radio wavelength of 100 m: if the zone is elongated to a half length of 15 km 
small compared to the wavelength of a ripple, the amplitude is increased by : 
‘factor of three. For a rough ionosphere at the same height, for which 05=5°; 
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the radius of the effective area for reflection is 20 km: if the amplitude is to be 
increased by the same factor, the half length of the effective area has to be 
increased to 180 km, which is no longer small. Thus it is to be expected that 
equation (1) will apply to a rough ionosphere when the amplitude increases above 
normal are quite small: for large values of the right-hand side of equation (1), 
the amplitude will be less than that given by the equation. 


Fig. | 


IV. PROPOSED EXPERIMENTAL TEST 

The difficulty of testing the relationship is that, under rapid fading con- 
ditions, the changes in phase path are also rapid, and the usual method of 
measuring phase path (Findlay 1951) is unsatisfactory, partly because of the 
tedious work involved in counting wavelength changes from a confused film 
record and partly because the large and rapid changes in echo amplitude make it 
difficult to obtain a continuous record. 

The Findlay method consists of beating an echo from a pulse transmitter 
with a ¢.w. oscillator output phase locked to the transmitter at a frequency of 
about 50 ke/s greater than the transmitter frequency. The rectified output is 


differentiated and recorded on film. 
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It is proposed to overcome the difficulties associated with this method by 
feeding the rectified output into a phase-sensitive rectifier whose comparison 
phase is derived from a 50 ke/s oscillator phased with respect to the original 
transmitted pulse by the rotation of a condenser. The output from the phase- 
sensitive rectifier will be used to drive the condenser in such a direction as to 
reduce this output. The rotation of the condenser will be recorded: this gives 
the changes in phase path. To ensure that a sufficiently pronounced beat is 
produced at the receiver output, the output amplitude of the c.w. radio-frequency 
oscillator will be made proportional to the echo amplitude. 


This system of measurement should also enable us to distinguish between 
the rapid fading produced by small irregularities and that produced by a deep 
phase screen of the type discussed by Hewish (1951). 
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GEOMAGNETIC MICROPULSATIONS 
By GBA BLiis* 
[Manuscript received July 11, 1960] 


Summary 


This paper describes simultaneous observations of geomagnetic micropulsations 
at three places ranging from 28°S. to 51°S. geomagnetic latitude. It is shown that 
there is no observable change in the micropulsation period with latitude although there 
is a monotonic increase in the amplitude with latitude for all periods between 10 and 
100 sec. The interpretation of these results in terms of existing theories is discussed. 


I. INTRODUCTION 

The study of the micropulsations of the geomagnetic field has been intensified 
in recent years with the realization that these phenomena may provide evidence 
for the propagation of hydromagnetic waves in the outer atmosphere. Dungey 
(1954), for example, has suggested that micropulsations may be caused by 
standing hydromagnetic waves along the geomagnetic field lines, oscillating 
perpendicular to the meridian field. In the fundamental mode the period of 
these oscillations would increase with geomagnetic latitude as the length of the 
lines increased. Obayashi and Jacobs (1958) have reported that in middle and 
high latitudes the period of some micropulsations appears to increase with 
latitude and they calculated that the density of the outer atmosphere needed to 
support standing wave oscillations of the periods observed agreed fairly well 
with that obtained from the dispersion of whistling atmospherics (Storey 1953). 

Dungey (1954) has suggested also that interplanetary gas, flowing over the 
surface of the Chapman-Ferraro boundary separating the terrestrial and inter- 
planetary atmospheres, would generate waves in a similar way to the generation 
of waves on water by the wind; propagating downwards as Alfvén waves, they 
would reach the ionosphere in the auroral regions. If their effects were localized, 
micropulsations due to this mechanism would not be observed at low and middle 
latitudes. However, it has recently been shown by Bomke et al. (1960) that 
magnetic disturbances caused by the Argus high altitude atomic explosions 
appeared to be propagated horizontally in hydromagnetic ducts between the 
ionosphere and a height of 2500 km. Such propagation would also be likely to 
spread micropulsations. In this case there would be a change in amplitude with 
latitude but not in period. It is relevant to this mechanism that observations 
at the boundary of the geomagnetic field with the Pioneer I space probe have 
shown the existence of large field fluctuations in this region (Sonett, Judge, and 


Kelso 1959). 


* Upper Atmosphere Section, C.S.1.R.O., Camden, N.S.W.; present address: Department 
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A third alternative has been put forward by Lehnert (1956) and Maple 
(1959), who, following Holmberg (1951), have suggested intra-layer hydro- 
magnetic oscillations within the ionospheric # and F layers. The main objection 
to this hypothesis is that the high attenuation of hydromagnetic wave propagation 
under normal ionospheric conditions would appear to preclude such resonance 
effects (see, for example, Piddington 1959). 

For the systematic study of these theories it would be necessary to record 
micropulsations simultaneously over large parts of the world. Nevertheless, 
it is clear that useful information would be obtained from a more limited investi- 
gation by comparing the amplitude and period of micropulsations recorded 
simultaneously at places in the same geomagnetic meridian, particularly at middle 
and lower latitudes where the disturbances of the auroral zone are less pronounced. 


This paper presents the results of a series of observations over a range of 
geomagnetic latitudes extending from 28° S. to 51°S. for the period September 
1959 to April 1960. 


IT. OBSERVATIONS 

Micropulsations were recorded simultaneously at Townsville, Qld., Camden, 
N:S.W., and Hobart, Tas. These three places are at nearly the same geomagnetic 
longitude, but have geomagnetic latitudes of 28° 8., 42° S., and 51° S. respectively. 
Each recorder used a 4-turn horizontal pick-up loop with an average area of 
10 000 m2. The loops were connected through galvanometer-photocell amplifiers 
to pen recorders with a chart speed of 6in/hr. The frequency response of the 
amplifiers was flat from d.c. to 0-15 c/s, while the highest frequency which could 
be resolved by the recorders was 0-17 c/s. The deflection sensitivity of the three 
recorders for full scale deflection was normally 0-7, 0-5, and 0-3 y/sec 
respectively, and the amplifier noise was less than 0-005/see (1 y=10-° gauss). 
Each recorder was calibrated periodically, both with a calibration loop in the 
same trench as the pick-up loop and with a signal generator applied to the input 
terminals of the amplifier. The peak-to-peak amplitude of the third highest 
oscillation and the average period in 5-min intervals was measured, except when 
the period was greater than 50 sec, when 6 cycles were used. One hundred and 
nine individual occurrences of micropulsations occurring on 40 days during 
February, March, and April 1960 were used in the analysis ; 74 of these were 
in the day-time and 35 at night. Records obtained during magnetically active 
periods (,>6) were not used because of the difficulty of separating the many 
different frequency components without a spectrum analyser. 


Although a superficial inspection of the simultaneous records showed many 
occasions when there was an apparent large increase in period with latitude, 
closer examination showed that in all cases this was due to the superimposition 
of micropulsations of different periods but of different amplitudes at the three 
places. In general, those of longer period were stronger (compared with those 
of shorter period) at Hobart than at Townsville. This was shown clearly on 
occasions when bursts of relatively pure oscillations of different periods were 
recorded in succession. The periods were then always observed to be the same 
at all stations within the limits of measurement, which were estimated to be 
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+5 per cent. On the other hand, the change in amplitude was often pronounced, 
particularly for the longer periods. Figure 1 shows sample records of micro- 
pulsations at Townsville, Camden, and Hobart on April 13, 1960. The two 
bursts of micropulsations at A and B on these records had periods of 30 and 17 sec 
respectively. The relatively greater amplitude of the shorter micropulsations 
at Townsville should be noted. Figure 2 shows the average amplitudes of 
simultaneous micropulsations of different periods at the three places and Table 1, 
the average ratios of their amplitudes at Hobart and Townsville. 


TABLE 1 
THE RATIO BETWEEN MICROPULSATION AMPLITUDES AT HOBART AND TOWNSVILLE 


Period range (sec) ye 10-20 20-30 30-45 45-60 60-100 


Amplitude at Hobart 
Amplitude at Townsville 


1-7540-8 | 1-7540-5 | 4-041-5 | 4-141-25 | 4-041-4 


No. of observations .. 28 24 14 17 26 


III. Discussion 
It is clear from the observations that the variation of period with latitude 
was very small between 28°S. and 51°S. To see how this result compared with 
the predictions of the standing-wave theory, the fundamental period of oscillation 
of the geomagnetic field lines was calculated using recent data on the density of 
the outer atmosphere. For a line of length /, passing through a fully ionized 
medium of mass density p, the period is 


L 2ds 
r= vay 


where ds is the line element and V(s) the Alfvén velocity, 
V(s)=H(s)(4re(s)-*. (Dungey 1954) 


H(s) is the magnetic field intensity. In the lower parts of the outer atmosphere 
the density of neutral atoms is not negligible and it is necessary to modify the 
expression for the velocity. Piddington (1959) has shown that the propagation 
constant is then 

K=H-1(4rp)+4(w/2t)#(1 2), 


which describes waves travelling with velocity 
H(4m9)-#(2e7)* 


and with attenuation 
=H-1(4re)(w/2t)*. 


+ is the period needed to accelerate neutral atoms to the same velocity as the 
ions (see Piddington 1957, Appendix I). 


AMPLITUDE (y/SEC) 
+0:25+ =+0:35> 


*+*O0-2—e 


Fig. 
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1.—Records of micropulsations at Townsville, Camden, and Hobart on April 13, 
1960. The individual oscillations were resolved on the original records. 
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Fig. 2.—The observed variation in amplitude with geomagnetic 
latitude of micropulsations with different periods. 
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Fig. 4.—The variation of the Alfvén velocity along different 
geomagnetic field lines. 
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Fig. 5.—The variation of hydromagnetic wave velocity below 
1000 km height when the effect of neutral atoms is taken into account. 
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Fig. 6—The attenuation below 1000 km height for a wave 
frequency of 0-01 c/s when neutral atoms are taken into account. 
The attenuation is proportional to (frequency)?. , 
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The velocity of propagation along field lines ending at different latitudes 
was computed assuming a dipole form for the geomagnetic field and using the 
model of the outer atmosphere recently proposed by Johnson (1960) (Fig. 3). 
Figures 4 and 5 show the velocity profiles obtained. Figure 6 shows the high 
attenuation for propagation within the ionosphere when neutral atoms are taken 
into account, and Figure 7 shows the variation of the fundamental period with 
latitude. 

With this model the period would be 30 sec at 28°S. increasing to 55 see 
at 51°S. An earlier exponential model of the outer atmosphere with greater 
densities between 2000 and 10 000 km heights (Obayashi and Jacobs 1958) gave 
periods of 45 and 70sec at 28°S. and 51°S. respectively. 


300 
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Fig. 7.—The calculated variation of the fundamental period of 
oscillation of geomagnetic field lines using the model of Figure 3. 


The absence of any observed variation in the period over this range of 
latitudes would therefore appear to point to some other cause of micropulsations. 
The monotonically increasing amplitude with latitude, on the other hand, 
suggests strongly that they originate at higher latitudes and are propagated 
towards the equator, with greater attenuation at the longer periods. The rapid 
decrease in the propagation velocity of longitudinal hydromagnetic ERE below 
3000 km (Fig. 4) supports the idea of suitable horizontal hydromagnetic ducts 
such as Bomke et al. proposed. Although the fundamental periods of standing 
waves along high latitude field lines would be much greater than the periods 
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observed, it should be noted that the considerable variation in velocity along 
these lines would favour the generation of many higher order oscillations with 
periods not integral submultiples of the fundamental. 


As in the case of Dungey’s hypothesis of ripples in the Chapman-Ferraro 
surface, micropulsations due to higher order oscillations of the longer field lines 
would first appear in the auroral zone. Lower latitude observations could not 
distinguish between these two possibilities. However, very accurate measure- 
ments of the phase difference between micropulsations at slightly different 
latitudes should provide information about their horizontal propagation. 
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PHOTOMETRIC OBSERVATIONS OF 5577 A AND 6300 A AIRGLOW 
DURING THE I.G.Y. 


By R. A. DUNCAN* 
[Manuscript received July 25, 1960] 


Summary 

Twenty-one months’ observation of the airglow from near Sydney, Australia, 
shows that (a) aurorae are detected 10° above the southern horizon at 6300 A whenever 
the magnetic disturbance index (K) reaches 5 but K must reach 7 before detection is 
certain at 5577A; (b) the 6300 A zenith intensity increases rapidly with K once this 
equals or exceeds 4, but the 5577A zenith intensity is independent of magnetic 
disturbance ; (c) the zenith intensity of 5577 A tends to be a maximum at 03 hr local 
time ; (d) the zenith intensity of 6300 A drops rapidly from dusk till 01 hr and then 
rises till dawn. 


I. INTRODUCTION 

As part of the I.G.Y. programme a nightglow photometer, lent by the 
American National Bureau of Standards, has been operated at Camden, near 
Sydney (geographic lat. 34°S., geomagnetic lat. 42°S.), since March 1957. 

The intensity of the green oxygen emission at 5577 A was monitored on 
280 clear nights from March 1957 till July 1958. After July 1958 attention was 
concentrated mainly on the red oxygen emission at 6300 A. Sixty-six clear 
nights’ observation of the red nightglow were obtained to the end of 1958. 

In the present paper these observations are studied for dependence on 
magnetic and ionospheric parameters and for nocturnal behaviour. The results 
are compared with those of similar studies in the northern hemisphere. 


II. THE INSTRUMENT 

The photometer has been described by St. Amand (1955). It has high 
spectral purity, successfully eliminating the continuous background even during 
bright moonlight. It scans the sky in a series of horizontal circles at successive 
zenith distances of 80, 75, 70, 60, 40, and 0°, a complete sky scan taking about 
four minutes, and usually being repeated each quarter hour. 

The photometer sensitivity was calibrated by Roach (1958) against a portable 
standard photometer, thus permitting comparison with nightglow photometric 
_ observations in other parts of the world. 


III. PHOTOMETRIC OBSERVATIONS OF AURORAS 
Auroras show clearly on the photometric records as pronounced brightening 
to the south or as structures such as ares and rays (Duncan 1959). Figure 1 
shows the probability of detecting an auroral form, as a function of the magnetic 


* Upper Atmosphere Section, C.S.I.R.O., Camden, N.S.W. 
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disturbance index (K). It will be seen that the red line is much more prone to 
auroral excitation than the green. Auroras are always detected at 6300 A 
when the magnetic index K reaches 5 but AK must reach 7 before we can be sure 
of detecting an aurora at 5577 A. This accords with the common observation 
(e.g. Seaton 1956) that middle latitude auroras are predominantly red. 


Q 
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GREEN (5577 A RED (6300 & 
MAR. 57-JULY 58 JULY-DEC 58 


80 


60 


40 
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IN WHICH AURORAS WERE DETECTED 


MAGNETIC DISTURBANCE INDEX (K) AT WATHEROO 


Fig. 1—The percentage of 3-hourly observing periods during which auroral 
forms were detected by the photometer against the magnetic disturbance 
index (K) at Watheroo. 


IV. AIRGLOW AND MAGNETIC DISTURBANCE 
It may be asked whether auroras make an appreciable contribution to the 
zenith airglow intensity. A study of the relation between airglow and the 
magnetic disturbance index K probably has some bearing on this problem. 
It will be seen from Figure 2 that, except for the sporadic effect of a few 
great auroras, the zenith green airglow intensity is independent of the magnetic 
disturbance index K. 
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Fig. 2.—The mean zenith airglow intensity against the magnetic 


disturbance index (K) at Watheroo. 


The zenith red airglow intensity, on the other hand, while independent of 
the magnetic disturbance index K so long as this is less than 4, increases rapidly 
with larger K. It would seem that the red airglow has a basic non-auroral level 
of about 175 rayleighs, plus a substantial auroral component on magnetically 
disturbed nights. 

Recently, Roach (1960) has shown that, although the mean intensity at 
5577 A is not significantly increased, there is a tendency for sporadic high 5577 A 
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TABLE 1 
LOCAL TIMES OF MAXIMUM 5577 A EMISSION 
Station Latitude tame of Observers 
Maximum 

Haute Provence, France .. 44°N. BYOH| Barbier, Dufay, and Williams (1951) 
Mt. Elbrus, Russia a 43°N. 01-6 Rodionoy, Pavlova, and Rdultoskava 

(1949) 
Cactus Peak, U.S.A. ie 36° N. 00-2 Roach, Williams, and Pettit (1953) 
Flagstaff, U.S.A. .. ¥ 35° N. 01-5 McLennon, McLeod, and Ireton (1928) 
Camden, Australia Bar 34° 8. 03 Duncan 
Sacramento Peak, U.S.A...° Sole. 23-01 Manring and Pettit (1958) 
Poona, India “t A 18° N. Minimum | Karandikar (1934) 

Ol 


intensities to occur at times of high K-index. Camden data support this con- 
clusion. McCaulley, Roach, and Matsushita (1960) have shown that the lack of 
correlation between 5577 intensities and magnetic disturbance index K is largely 
due to the imprecise nature of the K-index. They found a good correlation 
between 5577 intensities and the horizontal component of the geomagnetic field 
measured on a nearby magnetometer. 
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Fig. 3.—Mean nocturnal variation of the zenith green (5577 A) 
airglow. 


V. NocruRNAL VARIATION OF THE GREEN (5577 A) AIRGLOW 

The nocturnal variation of the zenith green airglow has been studied by a 
number of workers. Their results are summarized in Table 1. Observations 
at Camden confirm previous reports that the diurnal variation varies greatly 
from night to night. The mean nocturnal variations for each season, however, 
show a consistent pattern (Fig. 3). In winter a maximum is found ab 03 hr. 
In summer observations run only from 21 to 03 hr, but within this time interval 
the nocturnal behaviour seems similar to that found in winter. 

Observers at middle latitudes in the northern hemisphere (Table 1) have 


found maxima between 23 and 01 hr. 
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VI. NocTURNAL VARIATION OF THE RED (6300 A) ATRGLOW 

The zenith red (6300 A) airglow shows a far more pronounced and consistent 
nocturnal variation than the green (Fig. 4). Earlier workers (Elvey and Farns- 
worth 1942; Barbier 1957a) have described a dusk and dawn enhancement of 
the red airglow, but as a transient effect superimposed upon, and easily distin- 
guishable from, the “ true” nightglow. It would appear from Figure 4 that the 
dawn and dusk enhancements are simply parts of a smooth nocturnal variation. 
If this variation were to be explained entirely in terms of resonance scattering 
of sunlight, the terrestrial atmosphere would need to be effective at heights 
as great as 2000 km, for the winter dawn increase is already apparent at 03 hr. 
The nocturnal variation (Fig. 4) is not symmetrical about midnight. It is about 
twice as bright at dusk as at dawn, and the lowest intensity occurs not at midnight 
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Fig. 4.—The nocturnal variation of the zenith red (6300 A) 
airglow intensity. Mean for 20 nights during July and August 
1958. 


_butat 0lhr. This suggests that, in addition to an airglow component dependent 
on the solar zenith angle, there is a component which decreases steadily from dusk 
to dawn. Possibly this is due to the dissipation (as airglow) of solar energy 
absorbed by the upper atmosphere during the day. 
Bates and Massey (1946) have suggested that both the red airglow and the 
decay of the F’, region could result from the reaction 


Ose} 6+ O-L Ogun eee (1) 


In support of this Barbier (19576) found, for the first half of the night, an empirical 
relation 


I=(5-83 x 10°)(f,)? exp {(—h’ —200)/88} 


between the red airglow intensity J and the ionospheric height h’ and critical 
frequency fy. A comparison of Camden ionosonde and airglow records gives no 
support to this finding. Certainly there is support of a kind in the fact that 
both the red airglow intensity J and f,F’, decrease during the first half of the night. 
However, the airglow intensity at a given hour varies by a factor of 2 or 3 from 
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night to night and these fluctuations bear no direct relation to ionospheric 
parameters. Part of the variability is due to magnetic disturbance. As Figure 2 
shows, the average red airglow intensity I increases with magnetic disturbance 
while, as is well known, magnetic storms cause f,F’, to decrease and h'F 2 to increase. 
These changes are in opposite directions to those required by equation (2). 


Barbier also found that the second half of the night was characterized by 
an enhanced glow which appeared to the north and gradually spread over the 
entire sky. This phenomenon has never been seen to the south at Camden. 


VIL. ACKNOWLEDGMENTS 
This work formed part of the I.G.Y. programme of the Upper Atmosphere 
Section of the Commonwealth Scientific and Industrial Research Organization, 
and was done under the direction of Dr. D. F. Martyn. The photometer was 
lent by the American National Bureau of Standards. 


VIII. REFERENCES 
Barsier, D. (1957a).—C.R. Acad. Sci., Paris 244: 1809. 
Barsier, D. (1957b)—C.R. Acad. Sci., Paris 244: 2077. 
Bargsier, D., Duray, J., and Wiriams, D. R. (1951).—Ann. Astrophys. 14: 399. 
Barrs, D. R., and Massgy, H. S. W. (1946).—Proc. Roy. Soc. A187: 261. 
Duncan, R. A. (1959).—Aust. J. Phys. 12: 197. 
Exvey, C. T., and Farnsworru, A. H. (1942).—Astrophys. J. 96: 451. 
KaARANDIKAR, J. V. (1934).—Indian J. Phys. 8: 547. 
McCautiey, J. W., Roaou, F. E., and Matsusurra, 8. (1960).—J. Geophys. Res. 65: 1499. 
McLennan, J. C., McLzop, J. H., and Irreton, H. J. C. (1928).—Trans. Roy. Soc. Can. 22: 397. 
Mawnrine, E. R., and Perrirr, H. B. (1958).—J. Geophys. Res. 63: 39. 
Roacu, F. E. (1958).—U.S8. Nat. Bur. Standards Rep. 5591. 
Roacu, F. E. (1960).—J. Geophys. Res. 65: 1495. 
Roacu, F. E., Wiriiams, D. R., and Perrit, M. P. (1953).—J. Geophys. Res. 58: 73. 
Ropronov, S. F., Paviova, E. N., and Rountosxcava, HE. V. (1949).—C.R. Acad. Set. U.RS.S 
6613) 55; 
Sr. Amann, P. (1955).—Ann. Géophys. 11: 435. 
Seaton, M. J.{ 1956).—‘‘ The Aurorae and the Airglow.”’ Pergamon Press : London.) 


A NOVEL TYPE OF HIGH POWER PULSE TRANSMITTER 
By K. LANDECKER* and K. 8S. IwrrE* 
[Manuscript received June 14, 1960] 


Summary 


A system of generation of radio waves is described which makes use of a symmetrical 
circular array of condensers charged through resistors and discharged through spark 
gaps in the manner of the Marx impulse generator. It is shown that exponential wave- 
trains of very high peak powers, of the order of 10,000 MW, may be radiated. The 
radiation resistance and radiation field of the structure are given and modification of 
the field pattern by parasitic elements is considered. Formulae and graphical aids 
are given which facilitate the design of such transmitters and experiments with various 
model transmitters are described. Consideration is given to circuit losses, particularly 
spark losses, and means are described to minimize these losses. 


I. INTRODUCTION 

In recent years a demand for high power pulses of radio-frequency waves 
has arisen in various branches of physics, as for example in ionospheric and 
cosmic research as well as in radio direction-finding, radar, and communication. 
However, the generation of large radio-frequency pulses by means of conventional 
transmitters is at present limited to peak powers of the order of 1 MW. Trans- 
mitters for powers much in excess of this limit would become prohibitively 
costly. The rapidly increasing difficulty in the construction of large pulse 
transmitters is not only governed by the increase in the size and cost of the 
constructional elements, that is, the transmitting valves and their associated 
circuits, the aerial system, and the coupling elements between the transmitter 
and the aerial system, but also by the fact that new phenomena come into play 
which are not present in transmitters designed for moderate powers. We here 
merely point out one particular difficulty, namely, the power loss from corona 
discharges in the atmosphere surrounding the aerial system. 

The purpose of the present paper is to report on the results of experiments 
with various models of a novel type of pulse transmitter, which indicate that it 
is quite feasible to generate pulses of a peak power of the order of 10 000 MW 
at a small fraction of the cost of a corresponding conventional transmitter. In 
addition we shall refer below (in Section VI) to certain aspects of this principle 
of wave generation which appear to indicate that there exists a fundamental 
upper limit to the maximum power that, with materials available at present, 
may be radiated from a point in space. 


II. DESCRIPTION OF THE PRINCIPLE OF TRANSMISSION 


The principle of transmission may be conveniently explained with the aid 
of the schematic diagram Figure 1. A number of sets of electrical components 


L, 6, R, and 8 are shown arranged in circular form. The eight capacitors 0 
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are representative of an arbitrarily large number of condensers whose series 
combination forms the tuning capacity, while the corresponding tuning 
inductances L are located symmetrically with respect to the condensers and 
spark gaps S and may, at higher frequencies, be formed solely by the connecting 
leads. In a limiting case they may even degenerate into the internal inductances 
of the condensers. The total tuning inductance is then the inductance of the 
circular current path through the condensers and Spark gaps, and it will presently 
become clear that this is the desirable condition for maximum energy storage. 
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Fig. 1.—Schematic circuit diagram of the transmitter. 


The condensers © are charged in parallel through the resistors k and discharged 
~ in series through the spark gaps S in the manner of the well-known Marx impulse 
generator (Edwards, Husbands, and Perry 1951; reee and Meek 1954). 
However, the structure described forms not only the tank prcaie of the an 
mitter but also is extended spatially so as to form at the same time a HNO ae 
dipole or loop aerial which itself radiates the eeu encisy A Vatiation 
of this principle using parasitic elements in addition to the main oscillatory 
circuit will be referred to below. 
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Anticipating the results of calculations and experiments quoted below, 
which show that this system of transmission may be made to function and that 
very large power pulses may in fact be radiated, we enumerate some advantageous 
properties of this scheme : 

1. The energy stored in the tuning condensers is transformed into radiation 
in a most direct manner. This outweighs the fact that higher electrical energies 
may possibly be stored in unit volume by other means, for example by homopolar 
generators and the like, than it is possible to store in condensers with present-day 
dielectric materials. 

2. In the presence of an atmosphere a magnetic dipole is inherently superior 
as a radiator to an electric dipole. A magnetic dipole fails at high power because 
it ultimately gives rise to an electrodeless discharge in the surrounding air. An 
electric dipole causes a corona discharge at the field boundaries (the ‘ends? 
of the aerial wire) long before an electrodeless discharge can be initiated. The 
circular array of condensers of Figure 1 approaches a perfect magnetic dipole as 
the number of condensers is increased and their individual capacities increased 
simultaneously beyond all limits in such a way that the LC product of the entire 
circuit remains at a constant desired value. In practice it is sufficient to sub- 
divide the condenser bank to such an extent that the individual condensers are 
able to withstand their charging potentials. The latter are the only scalar 
potentials that arise in the circuit. This statement is related to a well-known 
electrodynamic theorem which has been frequently mentioned in the literature 
(Howe 1945) and according to which no scalar potential difference exists between 
any two points of a conducting annulus situated in a varying magnetic field of 
axial symmetry. Experience indicates that with slight departures from a 
perfectly circular form and from perfect symmetry there is still little tendency 
for a corona discharge or flashover to occur. 

3. Since no point of the structure is necessarily at ground potential, it is 
possible to ground the electrical centre of the voltage supply. The latter therefore 
need only be insulated for one-half of the condenser charging voltage. 


III. CALCULATION OF THE POWER OF THE EMITTED WAVE-TRAIN 
The stored energy will be radiated in the form of exponentially decaying 
wave-trains. To obtain an estimate of the power that such a system is capable 
of radiating we begin by calculating three characteristic quantities of the emitted 
wave-train : the average power during the first half cycle, the peak power of the 
entire wave-train which occurs instantaneously at the current maximum of 


the first half cycle, and the approximate average power of the total wave-train. ° 


Assuming that N condensers each of capacity O are distributed around the 
circumference of the circuit, and that they are designed to withstand a maximum 
charging voltage V, the total energy Ho; stored in the condenser bank is related 
to the current in the circuit through 


Ext= iNOV?2= ib TG sin? wt Oty a eee Co (1) 
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where k=h/2L, L is the circuit inductance, w=(N/LC —h?/41L?)3, R is a resistance 
representing radiation and circuit losses, and J, is the current in the circuit in 
the absence of damping (R-+0). 


From equation (1) we obtain 


IpR 0 
e—2kl (ce) fee) 
|| | all e—2k cos 2at dé 
—9] 
0 0 
a co 
og =| age COS 7 Cat Ol meee eer eel (2) 
0 


and, integrating the second term by parts, 


NOVI 
TR. 2h! le? 


co J-2 
+- I elt COB 2 ON 0h ae es ce (3) 
9 o 


On multiplying (2) by k?/w? and adding to (3) there results 


Nove \_ 1 
ER (at) say 7 ee oC ier (4) 
from which follows 
UN VV Ghd Gees ee ates «ae See eee (5) 


and therefore the average power during the first half cycle is closely 
Po tary =FIRR=$FR(NVJOL)® oo eee eee (6) 


In all cases of practical interest w~(N/LO)!, that is, oL»N/o0. Introducing 
the Q-factor of the circuit Q=wL/R=N/wCR, the average power during the first 
half cycle may also be expressed as 


Po (avy =4NC V20/Q =Hiotw/Q watts, cra Dad Gao Dose (12) 


when C is given in farads, V in volts, and H.4 in joules. The instantaneous 
peak power at the first current maximum is double this value, that is, 


Po (peak) =2 Bio4o/Q Wich US Sameera elaieps a optics (8) 


Lastly, in order to determine the approximate average power during one wave- 
train, we assume that the wave-train may be considered as terminated when the 
current amplitude has decreased to 1 per cent. of its initial value (Fleming 1919). 
Under these conditions the duration z of the wave-train is given by 


e-F=0-01, ] 
or 
bln, 201 —4-6)| 


* The same expression for P, (ay) (except for a factor (l1—e—§/?) not very much different from 
unity) is obtained by considering the wave-train to consist of a succession of half sinoids of 
amplitudes I), I,, I,,. . ., where successive amplitudes are related by the logarithmic decrement 
S=n/Q through J,/I,=I,/Is=Jo/Ta=. . .=0°/?. 
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but, since 
bk=Rj/QDswl2Q, soe. ees ie Se ee eee (10) 


we obtain for the duration of the train the value 
7=4-6/k=4-6 X2Q/m=9-2Q/@~10Q/w seconds. .. (11) 
The average power of the whole train is therefore given with sufficient accuracy by 
Per (ayy = 4 NOV20/10Q =Etot/10Q watts. ........ (12) 


Finally we note that from (11) there follows the useful expression for the number 
of cycles in the train 


n=or/2r=9°2Q0/27r=1-460~1-bQ.” awa. - (13) 


It is necessary to show that it is sufficiently accurate to include the radiation 
resistance in the total circuit resistance R of equation (1) as a part that 
is independent of frequency. 

So far in this discussion we have assumed that the wave-trains are generated 
in the manner described in Section II, without going into the process of trans- 
mission in detail. This process depends fundamentally on the Fourier com- 
position of the train. Fortunately, the frequency composition of an exponentially 
decaying wave-train may be stated in very simple terms. It is well known that 
the Fourier transform of such a train is identical with the frequency response of 
a tuned DRC circuit.* However, the radiation resistance of a transmitting 
aerial is a function of frequency and also a function, unfortunately not at all 
simple, of the spatial configuration of the aerial. The general problem of trans- 
mission and reception of arbitrary wave shapes is in fact exceedingly difficult. 
However, it will be shown in the following section that, though the radiation 
resistance of a circular loop exhibits a general increase with diameter but oscillates 
for large diameters, it nevertheless shows a simple monotonic increase within the 
range of interest (see Fig. 2). For all loop diameters which are practical the 
radiation resistance increases approximately with the fourth power of the diameter 
(see Fig. 3). From the foregoing it follows that it is permissible to assume 
that the frequency components of the radiation field have appreciable values only 
in the neighbourhood of the frequency of oscillation of the wave-train, provided 
that the train is not too short. It has in fact always been found that under these 
conditions the radiation field of a spark transmitter is a time function closely 
similar to the aerial current. The circumstance that in reality, in the process 
of radiation from our transmitting loop, the higher frequency components are 
always favoured compared to the lower frequency components, makes the above 


assumption of minor importance in practice, but this simplification greatly reduces 
the complexity of the problem. 


* See, for example, the Fourier pairs Nos. 448 and 449 in Campbell and Foster (1931). 
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Finally, it should be noted from equations (7), (8), and (12) that the power 
of the wave-train is simply inversely proportional to the Q-factor of the circuit, 
which is solely determined by radiation and circuit losses. It will be shown in 
the following that the circuit may always be adjusted so that the radiation 
losses may have any desired value, and that the circuit losses may be minimized 
by various means. The question of which particular value of the Q-factor 
is desirable will largely depend on the specific application of the emitted energy. 
For example, the wave shape most suitable for producing changes in an ionized 
medium will not be the same as that of a signal to be received against a strong 
background of noise. In transmitting practice it is accepted that the loaded Q 
of a transmitter should be at least 10 in order to ensure that the major part of the 
wave energy is radiated near the fundamental frequency. It can only be estimated 
that in general the choice of the factor will range between values of 10 and 50, 
and in the following we shall consider numerical examples in this range only. 


IV. THE RADIATION RESISTANCE OF THE CIRCUIT 

In the previous section the significance of the Q-factor was emphasized. 
It is therefore of importance to determine what fraction of the total damping 
resistance is made up of undesirable circuit losses. As will be shown below, the 
major portion of these losses is accounted for by losses in the sparks and means 
will be mentioned to minimize these losses. In order to render the problem 
manageable we therefore initially disregard the circuit losses and deal first with 
the radiation resistance. 


It is not immediately obvious that when a desired frequency of transmission 
and a maximum condenser voltage are specified the circuit constants may be 
made to converge towards values which are practicable and which ensure a large 
output of radiated power. The frequency fixes in the first instance the LC 
product in which the factor C should ideally be obtained from the largest number 
of series capacities it is possible to accommodate along the circumference of the 
circuit, allowing sufficient space for the spark gaps. The factor L, on the other 
hand, is a function of the diameter of the circuit and depends only weakly on the 
average cross section of the current path. The diameter of the circuit, however, 
also determines the radiation resistance but, as mentioned in Section III, the 
latter unfortunately is not a simple function of the former. If a suitable radiation 
resistance (a suitable Q-factor) is chosen, then the loop diameter, the tuning 
inductance, and the tuning capacity are determined. The number (NV) of 
condensers is limited by space requirements and it might well have been that the 
~ energy stored in the condenser bank and therefore also the radiated power were 
small. That the power is in fact large can only be shown by numerical examples 
or by the dimensional reasoning adopted in Section VI. 

The radiation resistance of circular current loops whose diameter is an 
appreciable fraction of a wavelength has been calculated by various authors 
(Foster 1944; Moullin 1946) and, in addition, the field distribution around such 
loops has also been worked out. 
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The radiation resistance is given by 
2k 
R*— bye | J,(~)da ohms, ...--------- (14) 

0 
where p)=1:2610-* henry/m, c=velocity of light—3 %108 m/sec, k=2na/a, 
a—radius of loop, A=wavelength, and J, («) is the Bessel function of the first 
kind and second order (Fig. 2). 
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Fig. 2.—The radiation resistance of a circular loop as a function of its 
diameter (in wavelengths). o indicates the region of interest (see 
text). 


This may also be written in the form of a power series in k as 
k2 k4 
ee 20mm (1— sta. . ); cbt vee (14a) 


which shows that the radiation resistance of a loop which is small compared to 
the wavelength approximates to 


A tld Vepsed et PA Utd (HEN Gee Se As aaa. 4 (14b) 
Equation (14b) still holds within 20 per cent. up to a value of k=1. Using the 
three terms of the power series in (14a) gives R* within a few per cent. at k=1. 


From equation (14) the radiation resistance R* of the loop has been plotted 
as function of the diameter in wavelengths in Figure 3 over the range which will 
be shown to be of interest to us. 
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The inductance L of the circular loop which together with the resistance 
A* of equations (14), (14a), or (14b) enters the Q-factor is given by 


Teen (ead 21-73 ee (15) 


where a=loop radius, as before, and )=(average) radius of the cross section 
of the current path. 
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Fig. 3.—The radiation resistance of the loop within the range of 
interest (see text). 


Since the ratio a/b and therefore also the term In (8a/b) are restricted by 
practical considerations to a range with rather definite limits the inductance may 
be written in an approximate form 


where D is a constant. 

Since the loop radius a enters equations (14b) and (15a) through the fourth 
and first power respectively it follows that a desired Q-factor may always be 
realized. Using (14) and (15) the Q-factor may be written 


o=! (m7 = —1- v5) \ {Je (a fama acetl (16) 
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Fig. 4.—Nomogram relating the circuit diameter as a fraction of the loop diameter, 
the Q-factor, and the loop diameter in wavelengths. 


This equation may be put into the form 


log Q=log f(b/a)+log g(2a/A), «2... eee ee ee (16a) 


where f and g are functions following from (16), and therefore a nomogram of Q, 
b/a, and 2a/\ may be constructed as shown in Figure 4. The nomogram shows 


at a glance that practical loop diameters will mainly be restricted to the range 
from 0-2 to 0°32. 


NOVEL HIGH POWER PULSE TRANSMITTER 647 


V. THE RADIATION FIELD OF THE LOOP AND THE EFFECT OF PARASITIC 
ELEMENTS 
The electric and magnetic vectors of field strength EH and B, at a distance r 
from the centre of a circular loop of radius a, have the following magnitudes 
(Foster 1944; Moullin 1946) 


1 
B= Spee bd,(k cos 9) = sin (o =) ero (17) 


and 


‘ov 
B= —}uol kJ, (k cos @) 7 sin (0 aa) Gore es (18) 


where I, is the amplitude of the current in the loop, rs~a, ¢ is the angle between 
the plane of the loop and the radius vector to the reference point, J,(#) is the 
Bessel function of the first kind and first order, and the other symbols have the 
same meaning as before. 

With the plane of the loop considered as an equatorial plane, the electric 
vector oscillates tangentially to a circle of latitude and the magnetic vector 
tangentially to a meridian. This should be noted when comparing the radiation 
field of the loop with that of an electric dipole. 

When the diameter of the loop is much smaller than a wavelength, that is, 
when k<1, these expressions approximate to 


E= tucel hk? cos 9 sin (o— ae ere ee se (19) 
j ope, 1 2 iar Wr 
B=—}y, Ik? cos ?- > sin ot — up sep creen (20) 


These approximations are correct to within about 10 per cent. even up to 
k=1, so that even when the loop is one wavelength in circumference the radiation 
field is still very like that of a small magnetic dipole. For reasons stated in 
Section IV this will nearly always be so for the transmitting loop we are 
considering. 

It should be noted that, whenever the Bessel function J,(x) passes through 
zero, extinction angles occur near the equatorial plane and progress towards the 
zenith with increasing radius of the loop until the entire quadrant is divided into 
lobes. The first zero of the Bessel function occurs for the argument 7=3°-83, 
that is, an extinction angle occurs for the first time when #=2ra/A=—3-83 or 
2a/A=3 -83/x=1-22. But since we are here always restricted to the condition 
0-2 <2a/A<0-3 this means that we have always only a single lobe. 

There are various possibilities of arranging parasitic elements (reflectors 
and directors) in the neighbourhood of the main loop for the purpose of modifying 
the radiation pattern. Some of these are shown schematically in Figure 5.* 
© represents a biconical reflector suggested by Moullin (1946) and D represents 


* Various configurations like the closed coplanar rings A and coaxial rings B call to mind 
various electric counterparts, such as broadside, endfire, Yagi array, etc., but it is doubtful 
whether any useful inferenc2s may be made from these analogies. 
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a parabolic reflector. It is to be noted that all axially symmetrical elements 
such as A and B share the property with the main circuit that no scalar potential 
arises between any two points of the circumference. 

The question whether it is possible to induce sufficient current into any of 
these parasitic elements and still effect a worth-while modification of the polar 
pattern by placing the element at a suitable distance is in general difficult to 
analyse and is best decided by experiment. The tuned secondary circuit 
mentioned in Section VII will also affect the radiation pattern to some extent 
and a choice may have to be made regarding its position with respect to the main 
circuit and whether more than one secondary circuit should be used. 


MAIN CIRCUIT 


Fig. 5.—Various forms of parasitic elements used in conjunction 
with the main driving circuit. 


VI. NUMERICAL EXAMPLES. SOME CONSIDERATIONS REGARDING THE 

MAxiImMuM AMOUNT OF WAVE POWER THAT MAY BE PRODUCED IN SPACE 

We now turn to some numerical examples and assume to begin with that a 
wave of a frequency of 1 Mc/s is to be produced and that the condensers will be 
designed to withstand a charging voltage of V=200kV. Condensers for such a 
voltage may today be considered to be standard items in engineering practice. 
The @-factor will be taken to be 10, and initially we shall ignore circuit losses. 
Selecting a value of (diameter of loop)/(wavelength) of 0-28, the nomogram 
(Fig. 4) gives a diameter of the current path of 0-06 of the loop diameter and we 
see that a circular structure of 80 m diameter and a cross section of the current 
path of 4:8m diameter is required. This would call for large tubular 
condensers with air or possibly polythene ribbon as a dielectric. The 
loop will have an inductance of L=170 yH and therefore requires a tuning 
capacity of O/N=1-5x10-4 uF. We now suppose that N50 condensers are 
distributed around the circumference, that is, that ~5 m of circumference are 
available for each condenser and its associated spark gap, though a design as 
shown schematically in Figures 6 (a) and 6 (b) would allow us to accommodate 


_—— 
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many more units around the circumference, even when the Spark gaps are operated 
at atmospheric pressure. The capacity of each of the 50 units is then 
C=50 x1-5x10-4=7-5 x 10-3 uF, and the capacity of the whole condenser 
bank NC=502x1-5 x10-4=3-75 x 10-1=0 -375 uF. The energy stored in the 
total capacity is given by 


Bop =4N OV? =0 -5ex 0-375 x 10-8 (2 x 10°)?=7 -5 x103 J. 
CONDENSER PLATES 
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Fig. 6.—Schematic diagram of space-saving arrangement of condensers. 


We then obtain from equations (7), (8), and (12) of Section III : 
Average power during the first half cycle 
Po ayy=7°5 X10? X 6-28 X 109/10 =4-7 x10 J/sec 
=4700 MW. 
Peak power at first current maximum 
Po (peak) =2P 0 (avy =9400 MW. 
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Average power during the wave train 
Ptr (avy =P (av) /10=470 MW, 


and equation (13) shows that the train lasts for n=1-5Q—=15 cycles. 

This example shows that it is quite feasible to obtain with this arrangement 
and a condenser voltage slightly larger than 200&V peak powers of the order of 
10 000 MW. The structure is, however, very large, even if it should be possible 
the reduce the dimensions by improvements in design and judicious use of parasitic 
elements, as discussed in Section V. Nevertheless, it covers no more ground 
than a conventional aerial array at this frequency. If used for ionospheric work 
it might be advantageous to erect the transmitter at or near the bottom of a 
natural valley in order to direct the radiation upwards. 

We next consider a more readily realizable example and assume a desired 
frequency of 30 Me/s, a maximum condenser voltage of 200 kV, and a (Q-factor 
of 15. Proceeding as before, we select from the nomogram a loop diameter of 
28m, which fixes the diameter of the current path at about 4 cm and its 
inductance at L=16 wH. This requires a tuning capacity of C/N=1-75 uyF. 
Accommodating N=40 condensers around the circumference allows 22 cm for 
the space to be occupied by one condenser and its associated spark gap, which is 
ample. Proceeding as in the previous example we obtain : 


C=40 <i-75 > 10-70 woe, 
NC=2800 pykF, 
Eo = 4+ X 2800 K10- 2 x 10°)7=56 J, 
Po (avy =56 X 6-28 X30 X 108/15 =700 MW, 
Po (peak) =2 X Po (ayy =1400 MW, 
Pe Gre Oey) WORN 
n=1-5 X15 22 cycles. 


We again note that very large powers are radiated and yet the circuit is 
very easy to construct. This example was selected because a transmitter 
approximating to these specifications is at present being designed and constructed 
in this laboratory. 

Finally, we consider the maximum peak power that it is possible to radiate 
with this type of transmitter. For this purpose it is convenient to rewrite 
equation (8) of Section III as follows 


Po wea == UN V7] 07 Wallies oe ee ee ee (8a) 


We now assume that for the construction of the condensers a low-loss 
dielectric material with a breakdown strength of A=5 x10? V/m is available 
(this is about the breakdown strength of polystyrene), that one-quarter of the 
circumference of the transmitting loop is taken up by dielectric, and that the 
other three-quarters is occupied by the spark gaps and the condenser plates. 
Tt will be shown in Section VII that it is quite feasible to control the spark length 
by pressurization of the spark gaps and that the above allocation of space for 
the spark gapsisample. The charging voltage of the condensers may be expressed 
as V=4naA/N and, since it follows from the foregoing that in any practical 
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design it is not possible to deviate much from the values 2a~0- 292, Y~10, 
Rk~100 Q, equation (8a) becomes 


Po (peak) = ($77aA)?/Q2R = (47 X 0-292A)2/Q2R 
=(1-2 x 107A)2/104=1012 watts. ........ (8b) 


From (8b) it follows that it should be possible to produce peak powers of 
the order of 104, 106, and 108 MW at wavelengths of 1, 10, and 100 m respectively. 
Even when circuit losses are taken into account ae calculation demonstrates 
that very large peak powers may be produced. Further, this result seems to us 
to have also a fundamental significance. Unless entirely new methods of 
radiating electrical energy from a point into space are discovered, equation (8b) 
may represent a limiting value for such radiation with present-day dielectric 
materials and insulating techniques. 


VII. EXPERIMENTS WITH VARIOUS SCALE MODELS 

In order to verify the feasibility of this principle of transmission and the 
calculations of the previous sections, experiments have been carried out with 
two transmitters. The first operated at a frequency of 10 Mc/s and used 6 
tubular air-insulated condensers. The average loop diameter was about 2m, 
that is, only 1/15 of a wavelength, and so the radiation resistance was very low. 
This transmitter was mounted in the open on a low wooden tower. It proved 
very difficult to protect it from the weather. A second transmitter was con- 
structed working at a frequency of 70 Mc/s. It uses 18 parallel-plate condensers 
with a dielectric of polystyrene and 18 associated spark gaps. The mean loop 
diameter is approximately 1m, which corresponds to a radiation resistance of 
about 100 Q. This transmitter is small enough to be operated in the laboratory. 
Both transmitters must be considered to be scale models because the condensers 
in the 10 Mc/s transmitter are not able to withstand a charging voltage of more 
than 5 kV and those in the 70 Me/s transmitter more than 15 kV. The charging 
resistances each had a value of 10 MQ in both models. In addition to experiments 
with these transmitters, many observations were made on simple tuned circuits. 


The results of these experiments may be summarized as follows : 


(i) It was found that the spark gaps of the transmitters could be adjusted 
to fire quite regularly and strong signals were received at distances of many 
wavelengths from both transmitters. It was immediately established that the 
transmitting loops exhibited the radiation pattern of a dipole and that the field 
vectors had the correct orientations. The signal strength was of the expected 
order of magnitude, but experience with these model transmitters has shown 
that accurate measurement of the field strength of short, strongly damped, high 
frequency pulses presents a special problem, and various detectors are being 
designed for this purpose. 

(ii) Near the transmitters the signals are easily displayed on the screen 
of a cathode-ray tube. Observations of the pulse length indicate that by far 
the largest part of the circuit losses is accounted for by spark losses, and various 
attempts were made to minimize these losses. The most effective means of 
controlling the pulse length proved to be the coupling of a tuned secondary 
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circuit with the main driving circuit. This of course used to be standard practice 
in spark telegraphy, but with the arrangement we are considering the problems 
are somewhat different from those encountered in spark telegraphy. The tuned 
secondary takes the place of one of the untuned parasitic elements of type A or B 
in Figure 5 (Section V). The radiation resistance of such a passive element 
by itself will be either exactly or nearly exactly the same as that of the main 
circuit. If an untuned parasitic reflector is coupled to the main circuit it modifies 
the radiation pattern and also increases the radiation resistance of the whole 
combination.* A tuned secondary, on the other hand, in addition reduces the 
energy losses in the whole system because while being driven by the main circuit 
it abstracts oscillatory energy from this circuit and continues to radiate after 
the sparks are extinguished. It was found that such a tuned secondary controlled 
the pulse length very effectively. Experiments with single tuned circuits at a 
frequency of 10 Mc/s showed that the pulse length could be extended by this 
means by at least a factor of 10 without significantly affecting the peak amplitude. 
Maximum improvement occurs at slightly below optimum coupling. Finally, 
it should be noted that the tuned secondary circuit shares with the main circuit 
the property that no scalar potential differences appear between any two points 
of its circumference provided the tuning condensers are able to withstand the 
oscillating voltage. 

The spark discharge was also studied at pressures higher than atmospheric 
and an adjustable spark gap was constructed which is capable of being pressurized 
up to 100 atmospheres. The most useful gas for this purpose proved to be 
nitrogen, in agreement with various statements in the literature. Replacing a 
spark gap working at atmospheric pressure in an oscillating circuit by a pressurized 
gap with the same breakdown voltage always increased the pulse length; that 
is, no combination of circuit constants was found in which pressurization reduced 
the pulse length. The actual increase in pulse length depends very much on the 
circuit constants and may be as much as by a factor of 2. Such a factor is very 
desirable but the main advantage of pressurization for us seems to be the 
possibility of reducing the space requirements of the spark gaps and, incidentally, 
also the noise from the sparks. Moreover, a design shown schematically in 
Figure 7, which is quite feasible at frequencies above 10 Mc/s, would protect the 
circuit from humidity and dust and from the effects of high altitude. We are 
at present constructing a 35 Mc/s transmitter according to this design. 


(iii) According to the calculations of Section IV the optimum diameter of 
the loop will practically always be within the limits of 0-2-0-3 of one wavelength 
so that the firing impulse propagates from one spark gap to the next along ie 
circumference in a finite time which is an appreciable fraction of a period. It 
was therefore anticipated that a separate firing impulse might have to be applied 
by some means from the axis of the loop, so as to arrive simultaneously at all 
the spark gaps. However, it was noted from the beginning that the firing of 


* Here and in the following the radiation resistance is referred to the driving current in the 
main cireuit. Detuning of the main circuit due to the coupled reactance is of no practical 
importance. ‘ a 
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all the gaps seemed to take place Surprisingly decisively, though one could not 
be sure that the circuit did not support some unwanted mode along its cireum- 
ference. This was tested, using the 70 Mc/s transmitter, by firing all gaps 
simultaneously by means of a flash of ultraviolet light from an auxiliary bright 
spark between aluminium electrodes situated at the centre of the loop. No 
Significant difference in operation could be detected with and without the 
auxiliary spark. It therefore appears that the desired principal current mode 
established itself within a sufficiently short time by the mechanism that is 
responsible for the operation of the Marx impulse generator. 


ORIVING ALUMINIUM 
CIRCUIT Vea 


Fig. 7.—Schematic drawing of cast aluminium envelope for pressurization 
of the circuit. 


VIII. Conclusion 

Summarizing, we can state that theoretical and experimental work on this 
problem has at the time of writing reached a stage where the mode of operation 
and the limitations of this transmitter are well understood and the properties 
of a practical design may be predicted with sufficient accuracy. We are at 
present constructing a 60 ft diameter parabolic reflector for a 35 Mc/s transmitter 
with the aim of radiating pulses with a peak power of about 1000 MW for a 
programme of ionospheric and cosmic research. 

In this work we have been led to numerous other problems whose investi- 
gation has had to be postponed. ‘These include the possibility of using lengths 
of transmission line in place of the condensers in the circular array and the 
possible use of electronic or ionic valves with negative current-voltage 
characteristics, instead of the simple spark gaps, with the aim of generating 
continuous or modulated waves. Though oscillating circuits in which the energy 
is fed into the tuned circuit in parallel are highly developed, there seems to be 
no reason why series-fed circuits should not find equally useful applications in 
certain cases. 

The principle of transmission described here has been the subject of a patent 
application by the University of New England (1958). It is hoped that the 
transmitter described will find use as a research tool. 
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VERTICAL CHARACTERISTICS OF TRAVELLING IONOSPHERIC 
DISTURBANCES 


By L. H. HEISLER* 
[Manuscript received December 17, 1959] 


Summary 

The vertical dimensions of travelling disturbance phenomena have been investigated 
by derivation of the associated true height distributions of ion density. The height at 
which they appear often has an upper limit which may fluctuate in height from day to day. 

It appears that direction of travel is related mainly to season, with the possibility 
of some additional form of solar control only evident at times of sunspot maximum. 

During a particular season, there is no obvious change of direction of travel with 
height in the height range under observation. 

No definite variation of speed with height is evident in the ionospheric region 
considered, which extends from 160 to 230 km. 


I. INTRODUCTION 

Munro (1950, 1953) has described disturbances on 5-8 Me/s fixed-frequency 
records, using a three-station triangulation to determine speed and direction, 
and has recently published (Munro 1958) comprehensive statistical data of 
eight years’ results. Munro and Heisler (1956a, 1956b) have described large 
disturbances which occur in variable-frequency ionosonde records, and Heisler 
(1958), comparing records from Australian ionosonde stations, has shown that 
these disturbances travel distances of at least 3000 km, with fronts possibly 
broader than 1000km. The existence of similar disturbances in the northern 
hemisphere, first mentioned by Munro (1957), has recently been confirmed by 
Valverde (1958), using backscatter techniques. 

Recently there has been an increasing appreciation by many workers of the 
importance of N(h) electron density profiles in ionospheric research. The usual 
ionosonde presentation, which relates virtual height to frequency, has many 
shortcomings when used for investigating ionospheric phenomena and, in fact, 
can give rise to entirely misleading results. For the purpose of many investi- 
gations, it is the knowledge of the true height at which a certain ionospheric 
event occurs, and of the change in electron density N at this point, which is of 
importance. In particular, recent probings of the ionosphere by rockets and 
artificial satellites have necessitated a knowledge of true heights of ionospheric 
layers for comparison with experimental measurements. Furthermore, the 
derived N(h) curve may be used to estimate n, the total ion content of the layer 
up to maximum ion density, and variations in this quantity give valuable informa- 
tion on ionization and recombination processes, as shown in applications by 
Ratcliffe (1951). 

Many methods have been developed to obtain the necessary N(h) curves 
from the available h'f presentation, and as early as 1938 Booker and Seaton 
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suggested a simple technique based on a parabolic distribution of electron density 
with height, neglecting the effect of the Earth’s magnetic field (Booker and Seaton 
1938). This has been superseded by the more accurate integral methods of 
Manning (1947) and of Kelso (1954), the latter method including the magnetic 
field correction. These analyses and others are discussed in an excellent survey 
by King (1957). 

The complicated and tedious calculations necessary, particularly in 
accounting for the magnetic field, have always been a deterrent in employing 
these methods for ionospheric research. However, Jackson (1956), King (1957), 
and more recently Duncan (1958), have developed techniques which are especially 
suitable for programming automatic electronic computers, and a program 
prepared by Duncan has been used by this laboratory in SILLIAC, the electronic 
computer situated at the Physics Department of Sydney University. The data 
consisting of virtual heights and corresponding frequencies from the ionosonde 
record, enough readings to sufficiently describe the curve, are supplied to the 
machine on punched paper tape. The results are given in true heights and 
electron density, which are then plotted to give the requisite N(h) curve. The 
complete analysis of one h'f curve is performed by the machine in approximately 
20 seconds. 

Approximately 200 of these calculations have been made during an investi- 
gation of travelling ionospheric disturbances. 

Published data refer mainly to horizontal dimension of movements. The 
present investigation considers vertical dimensions and variation of characteristics 
with height using improved methods for the deduction of true heights. It is 
to be understood that the term height as used in this paper always refers to true 
height in the ionosphere, not virtual height. 


II. OBSERVATIONAL DATA 

In any study of heights of F-region travelling ionospheric disturbance 
phenomena it is important that the observational data do in fact pertain to that 
ionospheric layer, and that the effects observed are not due to disturbances in 
lower parts of the medium through which the probing radio ray has to pass. 
Briggs and Spencer (1954) have indicated that this uncertainty about precise 
height of measured movement is an inherent disadvantage of the fading method 
due to Mitra (1949), as applied to F-region measurements. There is a similar 
ambiguity in observations using radio star scintillations, and the height at which 
contributing irregularities occur does not seem to be known with certainty. 


Speeds and directions of disturbances used in this paper are based on time 
differences of anomalies in F' traces recorded at three spaced stations, a technique 
described by Munro (1950). A feature of these records is that the manifestations 
at a particular station usually occur at different times on the o- and a-ray traces 
(Munro and Heisler 1956a, 1956b). Moreover, the time differences so obtained 
are consistent with expected separations of o- and «-reflection points in the F 
region. The observed effects, therefore, are almost certainly due to phenomena 


occurring in this region at a height corresponding to the reflection height of the 
probing frequency. : 


—— 
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Consideration must also be given to the validity of analysis as applied to 
the observational data. If 6 is the angle which the probing ray makes to the 
vertical and ¢ is the sweep time occupied by an ionosonde record, then in the 
integral method of actual height analysis d0/dt must be zero or a very small 
quantity. This ensures that, for each frequency, the integration process is always 
along the same path through the ionized medium. 


A travelling disturbance manifests itself on an ionosonde record as a com- 
plexity due in part to non-vertical reflection during portion of the frequency 
Sweep. According to the above hypothesis, analysis of such a record would not 
be valid. However, in the results which follow below, ionosonde examples for 
analysis have been carefully chosen, either immediately before the advent of a 
disturbance or during that phase of the disturbance where non-vertical complexity 
is absent. 


III. HEIGHT OF OCCURRENCE OF TRAVELLING DISTURBANCES 
(a) Piaed-frequency Observations 

Most of these observations at this laboratory have been made on a fixed 
frequency of 5-8 Mc/s. This immediately predetermines the height of observa- 
tion as that height in the ionosphere at which an electron density of 4:17 x105 
electrons/cm?® occurs, and will consequently vary with season, sunspot number, 
and time of day. The observations considered in this paper were made during 
the years 1952 to 1954 inclusive and a determination of heights from carefully 
selected disturbance-free ionosonde records shows that the median seasonal 
midday values are as follows: winter, 189 km; equinox, 197 km; and summer, 
205 km. In each case there is a scatter of heights 10 km each side of these 
values. 

Disturbances are identified on the film record as complexities in the PF, 
trace, and are due to off-vertical reflections from curved isoionic surfaces which 
possess an apparent velocity relative to the observing station. Munro (1953) 
and Price (1959) have shown that the appearance of the disturbance as a com- 
plexity depends on the relationship between the radius of curvature r of the 
distorted isoionic surface and its height h above the observing point. If r>h 
the disturbance passage will cause height and intensity changes in the record only, 
and these are not normally recorded as disturbance effects. Also, for a given 
value of r, the number of complexities and therefore recorded disturbances should 
increase with h. Absence of complexities therefore does not imply absence of 
disturbances and the film records must be carefully studied for height changes. 
Unless otherwise stated the statistics considered here concern only those 
disturbances which have produced complexities on film records. 


Munro (1958) has shown that the number of disturbances observed varies 
considerably from day to day and indicates that, while ease of observation 
depends to some extent on the ionization gradient at the time, this does not 
account entirely for the fluctuations. A variation in occurrence of disturbances 
is found to be often associated with change in height of observation. This is 
illustrated in Figures 1 (a) and 1 (b). These are N(h) curves for approximately 
the same time on consecutive days. On March 24, 1954 there were a large 
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number of 5-8 Mc/s disturbances, whereas on March 25, 1954 there were no 
5-8 Mc/s disturbances and very few height rises. It will be noticed that on 
the disturbed day the 5-8 Mc/s reflection height was approximately 23 km lower 
than on the quiet day. Munro (1953) has indicated that travelling disturbances 
are more readily recognized on ionosonde records in a region where electron 
density changes slowly with height. The resultant manifestations on the film 
record are then more likely to be complex and hence more easily recognizable. 
It is possible therefore that differences in electron density gradient may influence 
the number of observed disturbances. On this basis, since the gradient is steeper 
on March 25, it would be expected that more disturbances would be observed on 
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Fig. 1.—Electron density distributions illustrating relationship between height and 

occurrence of disturbances. (a) Camden 1243 hr, March 24, 1954, a day of numerous 

5:8 Mc/s disturbances. (b) Camden 1230 hr, March 25, 1954, a day of no 5-8 Me/s 
disturbances. 


this day, whereas the converse is true. In this particular case therefore it would 
appear that gradient is not a contributory factor and height of observation 
governs the number of disturbances observed. Several similar cases seem 
to%indicate that disturbances are more prevalent at lower reflection heights. 


More definite information on this point was provided by a series of observa- 
tions taken over a period of several months, using three frequencies simul- 
taneously. A typical record on October 2, 1956 showed disturbed conditions 
on 4°5 and 5:8 Me/s, while simultaneous observation on 9-8 Mc/s showed very 
few disturbances even when height rises as well as complexities were carefully 
checked on the film record to avoid electron density gradient influence on obser- 
vational results as discussed above. The N(h) curve of a typical ionosonde 
record for this period shown in Figure 2 places the 4-5 and 5-8 Mc/s reflection 
heights at 175 and 202 km respectively, whereas the 9-8 Mc/s reflection point 
is at the much greater height of 270 km. 


The conclusion from these fixed-frequency observations is therefore that 
in the height range under observation there is a definite tendency for disturbances 
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to be more prevalent at the lower heights, with the probability of an upper 
boundary to the height at which they are apparent. This has been previously 
Suggested by Heisler (1958). 


(b) Variable-frequency Observations 

The fixed-frequency technique described above is particularly valuable for 
studying small ionospheric disturbances, as it provides a very sensitive means 
of detection. It is limited, however, to observation of one isoionic contour and 
therefore cannot record the disturbance effect on a cross section of ionospheric 
heights. Large disturbances are more readily observed and their nature studied 
from sequences of ionosonde records made with a panoramic type recorder. 
Anomalies in such records due to travelling ionosonde disturbances have already 
been described (Heisler 1958). One particular type described therein as a C type 
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Fig. 2.—Hlectron density distribution at Camden for 1303 hr, October 2, 
1956, showing true height of 4-5, 5-8, and 9-8 Me/s reflection points. 


anomaly occurs as a cusp-shaped trace at the top of /, forming a double peak, 
and gradually travels down the Ff, trace. An N(h) curve of the ion density 
distribution giving rise to this type of disturbance is shown in Figure 3 (a). The 
F, peak in the original virtual height curve corresponds to an electron density 
of 3-1 105 electrons/em? and occurs at a beight of 160 km. On this day there 
were no 5-8 Mc/s disturbances or height rises on fixed-frequency records but 
frequent type C anomalies on ionosonde records, so, while disturbances occurred 
at a height of 160 km and below, there were no disturbances at the 5-8 Me/s 
height of 185 km. 

In this case the upper boundary must be close to the F, peak and below the 
5-8 Mc/s reflection height, i.e. at about 170km. By contrast, records several 
days later, of which Figure 3 (b) is an N(h) analysis, showed both Py and 5-8 Me/s 
disturbances to be present. In this case the F, peak reflection point is at 165 km, 
approximately the same height as before, but the 5-8 Me/s reflection point is at 
210 km; on this day, therefore, the upper boundary is much higher. It will be 
seen that the gradient is noticeably different on the two days and this may have 
some connexion with the change in boundary height. 
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IV. HEIGHT AND DIRECTION OF DISTURBANCES 
The variability of direction of travel of fixed-frequency disturbances has been 
previously described by Munro (1958). During the summer the mean direction 
is 120° east of north and during winter 30° east of north with an abrupt change 
from winter to summer conditions during September-October and a more gradual 
change from summer to winter conditions during March-April. 
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Fig. 3.—Electron density distributions illustrating change in upper bounding 
conditions for propagation of disturbances during equinoxial months. 
(a) Camden 1040 hr, March 18, 1952. (6b) Camden 1327 hr, March 22, 1952. 
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Fig. 4.—Variation of true height with direction of travel (-+-, summer 
values; @, winter values). 


Figure 4 shows a plot of direction against height of observation for 
disturbance cases in which summer values are plotted with crosses and winter 
values with dots. There are two obvious groupings; summer disturbances 
with heights above 190 km and directions ranging from 100 to 180°; and winter 
disturbances with heights below 190 km and directions ranging from 0 to 100°. 
This is a typical distribution, as summer heights of reflection for a particular 
frequency in general are higher than winter heights. However, the majority 
of winter disturbances observed at heights above 190 km also have directions 
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in the range 0-100°, and the few summer disturbances observed at heights below 
190 km still have directions in the 100-180° range. It would appear therefore 
that change in direction is mainly seasonal and it is not valid to associate it 
directly with change in height of observation. 


This apparent change of direction with height of observation can be further 
studied during equinoxial transitional periods when directions are variable. 
The N(h) curves for such a case are shown in Figure 5. Figure 5 (a) is an N(h) 
curve for 1151 hr, September 8, 1952, a day of predominant summer directions, 
and shows that a disturbance travelling 152° east of north at the time of this 
ionosonde recording was observed at a height of 195 km. However, September 11, 
1952 was a day of predominant winter directions and on a similar curve shown in 
Figure 5 (b) for 1301 hr a disturbance travelling 33° east of north was observed 
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Fig. 5.—Electron density distributions illustrating relationship between 
height and direction of travel of disturbances. (a) Camden 1151 hr, 
September 8, 1952, a day of predominant summer directions. (b) Camden 
1301 hr, September 11, 1952, a day of predominant winter directions. 


at a height of 177 km. It will be noticed that this height is 18 km lower than 
previously. This difference in observational height for disturbances travelling 
in different directions has been observed in several equinoxial examples. In 
all cases low heights of observation are associated with winter directions while 
high heights of observation correspond to those disturbances with summer 
directions. 

Since Figure 4 would suggest that during any particular season there is no 
direction gradient with height in the region of the ionosphere considered, it 
would appear that during the equinoxes there is some form of unstable control 
which not only causes the ion density distribution to fluctuate between normal 
summer and normal winter conditions but also governs the direction of movement 
of disturbances in these distributions. 

A similar fluctuation in disturbance directions has been observed during 
summer, but only during recent sunspot maximum years when on occasional 
days predominant winter directions occur. A case in particular occurred on 
January 19, 1957 and Figure 6 (a) shows a typical electron density distribution 
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just prior to a disturbance on this day. Figure 6 (b) is an NV (h) distribution 
prior to a summer direction disturbance on the next day. It will be noticed 
that the 5-8 Mc/s reflection heights are almost identical. Moreover, in contrast 
to the equinoxial case of Figure 5 there is no very marked change in overall 
electron density distribution. Several such examples have been analysed and 
show the same features. It would appear, therefore, that in these cases the 
factors controlling direction of travel mentioned in connexion with Figure 5 
are not significant and, since the particular phenomenon is observed only during 
sunspot maximum, some form of solar control is indicated. 
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Fig. 6.—Electron density distributions illustrating relationship between 

height and direction of travel of disturbances during summer at sunspot 

maximum. (a) Camden 1134 hr, January 19, 1957, 5-8 Mc/s disturbance, 

direction 40° E. of N.  (b) Camden 1115 hr, January 20, 1957, 5-8 Me/s 
disturbance, direction 144° EK. of N. 


V. HEIGHT AND SPEED OF DISTURBANCES 

As previously described (Heisler 1958), large ionospheric disturbances 
extend over a range of ionospheric heights and travel long distances with no 
apparent change in form or amplitude. Moreover, observations (Munro and 
Heisler 1956a, 1956b) indicate that a travelling disturbance always has an 
apparent vertical component of progression which is assumed to be the result 
of a forward tilt in the wavefront of the disturbance. If this were due to a 
height-speed gradient it would be difficult to understand how the disturbance 
could propagate without considerable change in form, particularly over large 
distances where tilt of the wavefront would become almost horizontai and hence 
vertical progression of the anomaly on the ionosonde record would be extremely 
slow. This is supported by the evidence of Figure 7, which shows the relationship 
between speeds of disturbances and height of observation. The random scatter 
of points suggests that in the region 160-220 km no correlation exists. 


Distance of travel and duration of small fixed-frequency disturbances have 
not been fully investigated, but initial examination of several cases indicates 
that they are much less than those of the large ionospheric disturbances. It is 
possible, therefore, that an individual fixed-frequency disturbance could possess 
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a speed gradient with height. This would not be obvious on the statistic plot 
of Figure 7 and would restrict the scale of propagation of the disturbance. 

Briggs and Spencer (1954) suggest that a height-speed gradient may exist 
at times of high magnetic activity, but, to arrive at this opinion, use Mitra 
method # and F region and radio star scintillation speed data. 

As emphasized previously, there is uncertainty about precise heights of 
measured F-region movements using the Mitra method and a similar ambiguity 
exists in radio star scintillation measurements. Moreover, there is doubt as to 
whether the same phenomena are being observed in each case and whether indeed 
these phenomena are the same as those measured at this laboratory by different 
experimental techniques. In our observations no evidence of such an effect 
has been found. 
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Fig. 7.—Variation of true height of occurrence with speed of disturbances. 


Ratcliffe (1954) refers to a height-speed gradient in the ionosphere of 
1msec-!tkm-}. This apparently originates in a private communication from 
Beynon and Thomas, quoted by Briggs and Spencer (1954). There is no indica- 
tion whether this refers to all disturbance phenomena or to a particular case. 

It is of little value to postulate a general gradient of speed with height in . 
the ionosphere when there is a high variability of observed speeds in the region 
160-220 km. It is probable that the inconsistency exists because of different 
phenomena being measured by different experimental techniques or is due to 
different geomagnetic observational locations. The use of virtual heights, in 
particular, can give very misleading results. 


VI. CONCLUSIONS 

The heights of observation of travelling disturbances and the mechanics 
of the phenomenon have been investigated by the study of associated N (h) 
electron density profiles. Results apply particularly to the sunspot minimum 
years 1952-1954 and are as follows: 

(1) Travelling disturbance phenomena are often contained by an upper 
boundary when anomalies are observed below this height only. The boundary 
is not always evident and may fluctuate in height from day to day. 
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(2) During a particular season there is no definite correlation between 
direction of travel of disturbances and height. Direction of travel seems to be 
mainly related to season, change in season being accompanied by a change in 
ion density distribution and a change in direction of travel of disturbances. 
Further, variability in direction of travel occurs during summer, but only in 
sunspot maximum years. There is little accompanying change in ion density 
distribution or height and the effect may be related to solar-geomagnetic 
phenomena. 

(3) There is no significant correlation between speeds of disturbances and 
height of observation in the ionospheric region considered, which includes heights 
between 160 and 230km. The observed long-distance propagation of large 
travelling disturbances also suggests that a height-speed gradient is improbable. 
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AN INTERFEROMETER FOR THE MEASUREMENT OF RADIO 
SOURCE SIZES 


By B. R. Gopparp,*+ A. WATKINSON,*{ and B. Y. Mitis*t 
[Manuscript received May 11, 1960] 


Summary 
Modifications have been made to the 85-5 Mc/s cross-type radio telescope at Sydney 
to permit the measurement of radio source sizes in the range of 10” to 1’ arc. The basic 
modification involves the addition of another aerial at a distant site connected by radio 
lnk. A new form of automatic gain control ensures very good stability for the system. 
The modified instrument is described in general terms and calibration techniques 
are discussed. 


I. INTRODUCTION 

Knowledge of the angular sizes of discrete radio sources is of the greatest 
importance as an aid to the identification of the sources, for an understanding 
of the basic physical processes causing the emission, and in the application of 
radio-astronomical data to cosmology. However, the measurement of size is 
usually difficult, and little reliable information is available. Some of the closer 
and brighter sources, particularly those of galactic origin, are quite large, of the 
order of a degree or more, and the measurements in this case are relatively easy. 
Several of these were detected during 1950 and 1951 in an interferometer survey 
of the sky (Mills 1952a). However, it was clear then that the majority of detectable 
sources were much smaller than this and that special equipment would be needed 
to resolve them. Such equipment, utilizing very widely spaced interferometer 
aerials, was constructed in both Australia and England and almost simultaneous 
measurements made of the strongest unresolved sources (Brown, Jennison, and 
Das Gupta 1952 ; Mills 1952b ; Smith 1952a). It was shown that the five sources 
measured, Cassiopea A and Cygnus A in England and Cygnus A, Taurus A, 
Virgo A, and Centaurus A in Australia, all had angular sizes of the order of a 
minute to several minutes of are, consistent with the identifications which had 
been made with visible nebulae. These initial measurements were extended and 
improved (Smith 1952b; Mills 1953; Jennison and Das Gupta 1956) and for 
several years represented the only results available on the small discrete sources. 
Subsequent measurements extended the number of known sizes to a couple of 
dozen, all in the range 1’ arc or greater (Carter, personal communication ; Edge 
et al. 1959) and, more important, showed that a few sources appeared to have 
sizes less than 12” are (Morris, Palmer, and Thompson 1957). 


* Division of Radiophysics, C.S.I.R.O., University Grounds, Chippendale, N.S.W. 
+ Present address: University of New South Wales, Sydney. 
+ Present address: School of Physics, University of Sydney. 
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On completion of the survey of the southern sky with the Sydney cross-type 
radio telescope, it became apparent that to progress in the identification of the 
radio sources and to apply their statistics in obtaining cosmological information 
it was necessary to obtain estimates of the sizes of a substantial proportion of the 
catalogued sources. When the catalogues are complete the total will exceed 
2000. 

Present results are based on the analysis of the initial catalogue between 
declinations +10° and —20°, containing nearly 1200 sources ; it appears that the 
majority are not resolved by the radio telescope, which is a pencil-beam instru- 
ment of beamwidth 50’ are. The instrument has now been modified, principally 
by the addition of another aerial at a great distance, to study the sizes of these 
unresolved sources. In the present paper we describe the modified instrument 
and its operation ; the principal observational results will be given elsewhere. 

The original instrument has been described in detail by Mills et al. (1958) ; 
briefly, it consists of a cruciform arrangement of two arrays of dipoles, each 
approximately 1500 ft long by effectively 15 ft wide. It operates at a wavelength 
of 3-5 m (85-5 Me/s) and by combining the outputs of the arrays in a switching 
arrangement, produces a pencil-beam response. Analysis of some unpublished 
results of A. W. L. Carter, who had performed a relatively low sensitivity survey 
using an interferometer with spacing between arrays of 1000A, suggested that 
the most useful arrangement for the present investigation would be initially an 
interferometer of about 3000 spacing. With this spacing, only sources with 
sizes of the order of 1’ arc or less would be visible, which would enable the 
recognition of the very distant and small sources. The relative visibility of the 
sources at this spacing would also allow crude estimates of angular sizes. 

However, estimation of the angular size of a radio source from measurements 
of its relative visibility using widely spaced interferometer aerials requires 
some assumption about the distribution of emission across the source. In so far 
as a good estimate can be made of this distribution the results of a measurement 
may be quite accurate; if, however, the distribution is complex or merely 
elongated, quite erroneous values may be obtained. In fact, some strong radio 
sources which have been investigated do show complex structure and are not 
circularly symmetric, so that for accurate results it is necessary to make measure- 
ments with base lines of different spacing and azimuths. It is expected that the 


present instrument will be used in this way at the completion of a sky survey now 
in progress. 


II. DESCRIPTION OF THE INSTRUMENT 
The modified system comprises, in effect, two interferometers using the 
existing east-west arm of the Cross as a common aerial. The other aerials 
consist of (a) portion of the existing north-south arm of the cross 300 ft in length 
centred 30A south of the centre of the east-west array and (b) a new array 
identical with (a) and situated 10 km distant in an east-west direction. The 
signal from aerial (b) is transmitted to the central site by means of a radio link 
similar to that used in an earlier investigation (Mills 1953). 
The general arrangement is shown in Figure 1, together with an indication 
of the response patterns of the interferometers. In each case the primary 
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pattern is formed from the product of the response of the two arrays in the 
interferometer, and consists of an ellipse 0°-8 by 3°-7 to the half response points 
in east-west and north-south directions respectively. This response is broken up 
into “lobes”, 2° apart in a north-south direction for the close-spaced ‘‘ local ”’ 
interferometer, and 1-2’ arc apart in an east-west direction for the wide-spaced 
“link” interferometer. The overall gains of these interferometers are made 
equal so that information about the angular size of a source may be derived 
from the ratio of their responses. 
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Fig. 1.—TIllustrating the layout of the arrays and the form of the resulting 
response patterns. 


The natural period of the interference pattern with the east-west interfero- 

meter varies according to the declination of the source, and is given by 
ie sec 5 hours, 
™ a 

where d is the spacing between aerials and $6 is the declination. With the 
operating wavelength 3-5m and the actual spacing of 10-2 km (2920) the 
fastest period is about 5 sec, occurring when a source is at zero declination. If 
this interference pattern were recorded directly, it would not be possible to 
determine the relative visibilities by direct visual observation, since a very short 
output time constant would be needed, with consequent reduction in visibility. 
To overcome this, the period is increased by continually sweeping the phase of 
one half of the interferometer, causing the lobe pattern to follow partially the 
source’s progression (Brown, Palmer, and Thompson 1955). In this way an 
interference pattern of any desired period can be obtained ; It was decided to 
use a standard 1 min period for all declinations. To enable direct comparison 
between the outputs of the long- and short-spaced interferometers, the period 
of the short-spaced pattern is also made 1 min by continuous phase sweeping. 
Some such arrangement is necessary, since a north-south interferometer has an 
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infinite period at transit and, with a fixed pattern, the response would depend 
on the position of the source relative to that of the pattern. 

Since the accurate measurement of relative visibilities at the two spacings 
is required, the amplitude of the interference pattern is of prime importance. 
It is therefore necessary that the gain should be stable in each part of the system. 
The degree of stability required necessitated continuous automatic gain control 
on all receivers. This is achieved by maintaining constant the amplification of a 
‘ control” signal injected at the aerial inputs. The gain is then independent of 
variations in receiver noise factor and interfering signals. 

Some of the essential parts of the system will now be discussed in more 
detail. 

(a) The Aerials 

The east-west arm of the Cross is used without alteration as the common 
element of the two interferometers. The constructional elements and per- 
formance of this array have been described elsewhere (Mills et al. 1958). Briefly, 
it consists of two collinear rows of half-wave dipoles separated by approximately 
half a wavelength and backed by a horizontal wire-mesh reflector. The overall 
length is approximately 1500 ft. 

The remaining element in each interferometer is a north-south array of 
full-wave dipoles similar in arrangement to the original north-south arm of the 
Cross. The spacing between dipoles is 4, and the phasing is arranged so as to 
produce a fan-beam response at right angles to the meridian plane. The position 
of this fan beam (i.e. its angle of elevation) is altered by a rearrangement of the 
phasing of each dipole by changing the point at which the coupling transformer 
is connected into the feed line, as in the original north-south arm of the Cross. 
However, the north-south arrays in this case have only 51, instead of 251, dipoles. 
With coupling and current distribution set as in the original north-south array, 
the aerial efficiency would be greatly reduced. To offset this, the coupling has 
been increased and the current distribution adjusted for rather less taper, resulting 
once again in an efficiency approaching 50 per cent. These alterations have 
caused an increase in the side-lobe responses, but they are only of the order of a 
few per cent. in the meridian plane and much less elsewhere: in view of the 
reduced sensitivity of the instrument compared with the original Cross, these 
are not serious. 

(b) The Recewers 

The overall receiving system consists of three separate receivers inter- 
connected as shown in Figure 2. The three receivers are identical except for the 
inclusion of a radio link in the remote receiver and a compensating delay line in 
one of the local receivers. Two signals are transmitted over the radio link, 
one is the original 85-5 Mc/s cosmic noise signal converted to 159-5 Me/s and 
re-radiated, the other is a 245 Mc/s c.w. signal derived from the local oscillator. 
When these are combined at the local site the original 85-5 Me/s signal is recon- 
structed with phase information unaffected by the conversion processes. The 
245 Mc/s signal also carries a pulse modulation for synchronizing the see 
gain control system as described later. The level of transmitted power is abou 


4 W in each case. 
D 
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The 85:5 Me/s receivers use double frequency conversion, first to 25-5 Me/s 
and then to 1-5 Me/s, the delay line being inserted at the latter frequency. The 
system bandwidth of 250 ke/s is determined at 1-5 Mc/s after combining the 
signals from each aerial; the bandwidth up to that point is greater than 1 Mc/s 
in all receivers. Phase switching and phase-sensitive detection is used in both 
interferometer systems, the phase-reversing switch being inserted in the local 
oscillator feed to the common east-west converter (85-5 Mce/s>25-5 Me/s). 
This corresponds exactly to the original Cross receivers (Mills et al. 1958). 


The lobe sweeping of each interferometer is accomplished by continuous 
phase rotation of the local oscillator signal fed to each north-south receiver. 
The phase changer consists of a pick-up loop rotating in crossed magnetic fields 
excited in quadrature. The loop is driven by a small synchronous motor fed 
from a variable frequency oscillator, to enable the period of the interference 
pattern to be kept at 1 min for all declinations. 


Square-law detection is achieved by an electronic multiplier at 1-5 Me/s. 
This multiplier is a balanced type to improve the square law and eliminate 
unwanted components (Coates 1957). 


(c) Automatic Gain Control 

The automatic gain controls operate by keeping constant at the outputs 
of the receivers a standard modulated noise signal injected at the input to each 
preamplifier. The ‘“ control” signal is derived from a square-wave-modulated 
noise generator. The average noise diode current is maintained constant by a 
feedback control system operating on the output of the oscillator supplying power 
to the diode filament. The fundamental component of the control signal in the 
receiver output is extracted from the remaining noise by a filter and phase- 
sensitive detector, and a feedback control system maintains the amplitude of 
this component at a preset level. This action, in conjunction with the control 
of the injected noise, maintains a constant overall gain which can be set to any 
desired value by altering the comparator voltage. Since it is the fundamental 
Fourier component which is measured at the output and not the average value of 
the “‘ square ’ wave, which is controlled at the input, it is important that the 
on-off ratio of the square wave be accurately maintained. For best performance 
this ratio should be unity and to achieve this a bistable multivibrator is danas 
triggered by a free-running pulse generator at double the frequency. 


In order that the input to the phase-sensitive detector be proportional to 
the gain of the system, regardless of the total level of the noise, it is necessary 
to extract the control signal after square-law detection. 


Referring to Figure 2, it can be seen that the output of each detector contains 
components of : 
(i) Control signal of the east-west receiver 
(ii) Control signal of the north-south receiver 
(iii) Coherent noise due to 429 c/s phase switching. 


To avoid interaction, a modulation frequency of 1080 ¢/s was chosen for the 
east-west receiver and 632 c/s for the north-south receiver. 
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The phase-sensitive detector used is a gating type switched by a square 
wave which must be identical to the one which modulates the noise diode. In 
the systems located wholly at the local site this is achieved by connecting the 
phase-sensitive detector and the noise modulator to the same Square-wave 
generator. In the case of the remote system it is necessary to generate another 
Square wave at the local site, which must be coincident with the Square wave 
modulating the noise diode at the remote site. 


The double frequency triggering pulses previously referred to are used to 
modulate the 245 Mc/s local oscillator of the link, and these pulses are then 
extracted at the local site and used to trigger a similar bistable multivibrator. 
A sensing device in the output of the phase-sensitive detector is used to ensure 
that the second square wave is generated in phase with the one at the remote site. 
This is necessary because the bistable multivibrator has an inherent 180° 
ambiguity of phase. 

In normal operation the automatic gain controls appear to have a stability 
of better than 2 or 3 per cent. over periods of 24 hours, with a long period stability 
over weeks of better than 10 per cent. With frequent calibration, as described 
in the next section, the relative sensitivities of the two interferometer receiving 
systems have an uncertainty of about 2 per cent. 


III. CALIBRATION 

Since the measurement of angular size involves the ratio of the responses 
of the two interferometers, only a knowledge of relative sensitivities is required. 
It is convenient to consider the two main parts of the system, the aerials and the 
receivers, separately. 

The two north-south aerials are made as nearly as possible identical. The 
methods of setting up have been described in detail elsewhere (Mills e¢ al. 1958) ; 
briefly, the current distribution is adjusted by feeding power from a CaM 
controlled oscillator into the aerial via its normal output connection, which is 
situated at the centre of the array. The current in each dipole is sampled in 
turn by a pick-up loop and fed to a detector. A Pei ceen oe signal is obtained by 
inserting a directional coupler into the feed line at a point Close to the centre 
of the array, i.e. near the feed point. This reference signal is passed through a 
variable attenuator and a variable phase shift to the same detector. The 
magnitude and phase of the current in any dipole relative to the current at the 
reference point can then be determined by adjusting the attenuator and phase 
shift for a null. Using this technique, the individual currents are adjusved to 
conform with the required law and, if the corresponding reference point is used 
when setting up each array, the dipole currents can be adjusted for approximate 
equality in the two arrays; any remaining plutprences care included i the 
overall calibration. Similar measurements made at different times have indicated 
an uncertainty of about 5 per cent. in the relative gains of the arrays. 

To check the receivers, a standard signal of a special type is injected into 
each preamplifier in turn and the output is measured on the recorder. The 
ratios of the deflections on the recorders then give directly the ratios of the gains. 
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In order to include correctly the effects of the output time constant on the 
responses, the calibration signal is arranged to simulate the standard 1 min 
interferometer pattern. The standard signal is an accurately controlled noise 
signal, square-wave-modulated at a frequency differing from the phase-switching 
frequency by 1 ¢/min, obtained by the use of a continuously rotating goniometer. 
Thus, the rectified receiver output has a component at a frequency which differs 
from the gating frequency of the phase-sensitive detector by 1 ¢/min, producing 
the required sine-wave output on the recorder. This calibrating signal, which 
except for the frequency chosen is similar to the A.G.C. control signal, is injected 
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Fig. 3.—Some calibration records—on the left is shown the two east-west calibrations per- 

formed by injection of the signal into the common preamplifiers and beside them the cali- 

bration of the local north-south system. On the right is shown the corresponding calibration 
of the remote preamplifier transmitted over the radio link. 


into the preamplifiers of each of the three systems in turn, i.e. the local and remote 
north-south systems and the east-west system. Typical calibrations are shown 
in Figure 3. When injecting into the north-south systems a pattern is produced 
only on the corresponding side of the record, but the common east-west system 
gives a response on both sides. 


IV. SizE MEASUREMENT 
As indicated above, the chief use of the instrument in its present form is to 
recognize radio sources of small angular size. A further aim is to obtain estimates 
of the actual sizes of the sources, but, because of the expected diversity of shapes 
and orientations, estimates obtained with two spacings only must necessarily 
be very crude. However, they have proved of some use in identification work 
and a study of their statistics should yield information of cosmological interest. 
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In order to represent the measured visibility ratios by the size of an equivalent 
model, we assume that all sources have circularly symmetric Gaussian brightness 
distributions, i.e. possibilities of fine structure, elongation, and asymmetry are 
ignored. 

This distribution is defined by 

Bie) by Bx (LCi Ua), hese tee ere cee (1) 


where r is the angular distance to the centre and 0, is the angle between half 
brightness points. 

It is well known (e.g. Mills 1953) that the relative visibility of a circularly 
symmetric distribution B(r), is given by 


Aa b 7(0) cos (Qnn0)d0, ....0-.....2.. (2) 
0 
1(0)=2 | G2—62)10" aie. 8° (66 (6 6 (4, 10; 8 «18 1610). 6 18.16 eo (3) 


where A, is the relative visibility at a spacing of n wavelengths. 
Inserting (1) in (2) and (3), it is readily shown that the equivalent size,of the 
model source is given by 
G,=31 -4.40e, X)* seconds of arc, .:....: 42.4; (4) 
where X is the ratio of visibilities at the two aerial spacings, i.e. XY = A49/A g909- 
This expression is plotted in Figure 4. 
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Fig. 4.—A plot of the equivalent source size against the ratio of 
responses on the two systems. 


For large values of X the model source becomes meaningless, as fine structure 
in the distribution will dominate the picture: the indicated size then will be 
usually smaller than the actual. Below values of XY ~3, however, the value of 
9, represents approximately the half brightness angle in an east-west direction of 
quite a variety of distributions. As Y approaches unity, the errors in measure- 
ment become important and, if the accuracy is about 10 per cent., as expected, 
a lower limit to the resolution is set at about 10” are. For the weakest sources, 
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statistical fluctuations reduce the accuracy considerably and, hence, also the 
effective resolution : for these sources the resolution limit is of the order of 30”. 
A substantial number of relatively strong sources have already been observed 
with apparently equal amplitudes at both spacings, i.e. X~1 +10 per cent. : 
for these we have 6)<10’. 

In Figure 5 a record of two sources is reproduced, showing, in one case, a 
strong source Hydra A (09-14, Mills, Slee, and Hill 1958), which has an appreciable 
angular size (XY¥=4-0 and §,=45”) and a weaker source (11—18), which is 
unresolved by the equipment (X ~1, 0,<10”). 
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Fig. 5.—Some sample records showing, on the left, a strong source resolved by the 
instrument and, on the right, a weaker unresolved source. 


Finally, one should mention the effects of the ionosphere on the measurement 
of size, since irregularities of electron density which cause the scintillation of 
radio sources might be expected to cause difficulties. If the ionosphere were 
sufficiently irregular to scatter coherent radiation into each of the widely separated 
aerials simultaneously, the whole basis of the method would be undermined : 
it is easily shown, however, from the known size of the irregula: ities that the 
scattering angle is too small by several orders of magnitude for this to occur. 
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Another possible cause of drastic error is the differential Faraday rotation 
at each aerial, but in this case, too, the known amounts of ionospheric refraction 
suggest that the effects would normally be negligible. One is therefore left with 
the usual effects of amplitude and phase variation, which, because of the great 
separation, are largely uncorrelated at the remote and local aerials. The latter 
variation is the more serious of the two: it can cause substantial changes in 
the period of the observed interference pattern and, occasionally, can be 
sufficiently rapid to destroy the pattern entirely, leaving it indistinguishable 
from random noise. To guard against the latter possibility it has been found 
desirable to take several observations at each declination and this has led to a 
considerable lengthening of the observational programme. Some of the results 
of this programme have already been described (Mills 1960) and others will be 
reported in subsequent papers. 
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A CATALOGUE OF RADIO SOURCES BETWEEN DECLINATIONS 
—20° AND —50° 


By B. Y. Mrts,**.O. B. SLEE,* and E. R. Hit1* 
[Manuscript received June 15, 1960] 


Summary 

A catalogue has been prepared of the radio sources observed between declinations 
—20° and —50°, using the Sydney cross-type radio telescope at a wavelength of 3-5 m ; 
a total of 892 sources is listed. This supplements an earlier catalogue in the declination 
zone +10° to —20°. In addition to the positions and intensities of the sources, angular 
sizes of 50 of the strongest are given ; several are found to have a size less than 15” are. 
As before, identifications with bright optical objects have been sought, and a number of 
galaxies of apparently abnormal radio emission listed. Statistical analyses of the 
distribution of the radio sources give results very similar to those obtained using the 
earlier catalogue. Within the uncertainty in the data, the distribution appears uniform 
in depth and there is a significantly greater number of sources of large apparent size than 
expected from chance blending effects. 


I. INTRODUCTION 

Analysis has been continuing of the records taken with the Sydney 3-5 m 
cross-type radio telescope. A catalogue of radio sources has now been prepared 
between declinations —20° and —50° to add to the previously published catalogue 
between -+-10° and —20° (Mills, Slee, and Hill 1958); the zone between —50° 
and —80°, which will complete the series, is in preparation. In the present 
catalogue a total of 892 radio sources is listed in the area of 2-66 steradians. 
The source density, 336 per steradian, is somewhat less than the previous figure 
of 356 per steradian, but not significantly so in view of the inclusion 
in the catalogue of the brightest regions of the Milky Way and the brightest 
southern extragalactic source: these cause a marked reduction in the listings 
of faint sources in their vicinity because of increased background temperatures 
and side-lobe effects. 

The catalogue has been examined in the same manner as the earlier one for 
identifications with optical objects and for the statistics of the radio source 
distribution. Since the present data add little to the earlier conclusions and 
since more data will be available soon from the final catalogue zone, these results 
are discussed only briefly. 

The present catalogue differs from the earlier one in that angular sizes are 
given for about 50 sources unresolvable with the pencil-beam aerial. These 
have been obtained with the newly developed angular size interferometer 
(Goddard, Watkinson, and Mills 1960). Since the sizes of many more sources 
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in the catalogue will be available shortly, the analysis and discussion of the 
results is postponed. It is interesting, however, that a number of radio sources 
are listed which could not be resolved by the instrument. These have been 
estimated to be smaller than 10-15” are. 


IT. THE CATALOGUE 

Preparation of the catalogue follows the methods described in two earlier 
papers (Mills and Slee 1957; Mills, Slee, and Hill 1958); subsequently these 
papers will be referred to as paper I and paper II respectively. The catalogue 
is divided into three zones covering declinations —20° to —30°, —30° to —40°, 
and —40° to —50°; they are given in Tables 1, 2, and 3. The previously used 
scheme of reference numbers has been adopted in which the first two digits of 
the number denote the hour of the Right Ascension ; these are followed by the 
sign of the declination and the tens digit measured in degrees and, finally, an 
italicized serial number arranged in order of increasing Right Ascension within 
the 1-hr period. Only the latter italicized numbers are given in the tables, as 
the others are evident; for example in Table 1, the second source would be 
referred to in the text as 00—22. 

As before, the probable error in the final digit of a position is indicated by 
a superscript. The positions given are directly as measured, with allowances 
only for precession. In the catalogue of paper II the sources of highest accuracy 
for which the estimated probable error in R.A. was --0-:1™ were corrected by 
+48, This correction was based on the differences between the apparent radio 
and optical positions of six bright identified galaxies. However, two of these 
galaxies are of large angular size and complex structure so that differences are 
probably not significant, and in two other cases our uncorrected positions were 
confirmed by other observers. It therefore appears probable that the systematic 
error in our positions is negligible and no such corrections have been applied in 
the present catalogue. Radio sources resolved by the aerial beam have been 
treated as before, both their peak flux density and their integrated flux density 
being given, the former in parentheses. In general, only sources with apparent 
angular size less than 2° have been included in the catalogue, two exceptions 
being the well-established sources at the galactic centre, 17—2/3, and 
Centaurus A, 13—42. When the size of a source has been measured using the 
angular size interferometer it is given in the footnotes in units of seconds of are. 
These are ‘“ equivalent” sizes and refer to a circularly symmetric source of 
Gaussian brightness distribution which has the same ratio of visibilities at spacings 
of 30A and 2920 (Goddard, Watkinson, and Mills 1960). The measurements 
are preliminary and are likely to be improved later. Possible identifications 
with bright nebulae are also given in the footnotes ; when the NGC number 
is in parentheses there appears to be no reason for thinking that this is other 
than an accidental coincidence in position. These and other possible identi- 
fications are discussed briefly in the next section. 

A small portion of the catalogue is included in an area examined by Rishbeth 
(1958), who gave contour diagrams and a short catalogue of radio sources in the 
Vela-Puppis region near 089 R.A. In general our catalogues agree quite closely, 
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Sources which may be “‘ extended ”, that is, resolvable, are indicated by a dagger. 
been placed beside uncertain flux densities. 
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TABLE 1 


A colon has 


Angular sizes, where available, are indicated in the 


footnotes. Details in Section II 
Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. Density 
No. R.A. Dec. (10-26 No. R.A. Dec (10 ee 
S. W m-? (e/s)—) Ss. W m-? (e/s)-) 
hm bas h m a 
00 01 
i 01-28 22 408 6 10 24-28 22 284 8 
2 02-18 23 227 10 oh] 28-7? 260215 18 
3 07-18 28 55§ 8 IZ 34-14 24 058 6 
4 10-08 23 43° 8 13 35-94 23 05° 5 
5 16-08 Dowote 13 14 39-24 Zh a 7 
6 ie 20 408 13 15 39-63 27/217 16 
if 20-7? 25 225 21 16 43-34 24 166 12 
8 21-24 PAU N55 14} Ly 48-7? 29 444 63 (43) 
9 21-98 29 434 33) 18 50-04 Dinos 11 
10 22-02 26 168 rice a 55-05 21 097 10 
hl 24-34 20 477 7 
12 24-34 27 518 5 02 
13 25-37 28 217 9 I 02-84 23 50° 10 
14 29-98 24 325 if 2 05-23 22 327 16 
15 30-54 22 097 8 3 08-54 28 Ol? 8t 
16 33-2? 20 184 19 4 09-34 23 488 9 
17 35-25 26 418 5 5 16-48 24 595 15 
18 35-5? 23 125 10 6 21-58 28 31° 11 
UD 36-94 29 247 8 rf 22-28 ras ee dS 19 
20 38-18 27 478 6 8 24-15 23 58? 10) 
21 42-0? 22 215 9 9 26-85 24 037 gt) 
22 45-1? 25 384 29°) 10 29-94 20 437 ll 
23 54-44 26 237 9 11 31-44 23 438 17 
24 56-53 29 218 8 12 33°5* 29 118 8 
13 39°34 28 128 Sf 
Ol 14 42-68 21 568 9 
w| 00-12 22 114 358) 15 44-08 23 57? 8 
2 00-63 27 438 16 16 45-42 29 287 10+ 
3 06-2 29 06’ 9 Ly 46-7? 20 427 14 
4 10-68 22 198 9+ 18 53-44 22 097 9 
5 13°78 23 O16 8 L9 54-22 23 328 28 
6 13-98 21 075 24 (18) 20 54-78 20 46° 13 
7 13-95 28 416 12 rail 54-78 26 187 a 
8 17-94 27 258 a 22 59-84 29 185 8 
a) 22-44 PAsy PASI 9 
@) 20". 
(2) NGC 253. 
(3) 30". 


‘9 Sources 02-28, 02-29 form perhaps one extended source. 
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TABLE 1 (Continued) 
| 
Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. Density 
No. R.A. Dec (CLORAS No. R.A. Dee (oss 
8. W m—(e/s)-) S. W m-~? (e/s)-1) 
hem pawn” hm oes 
03 05 
1 04-14 | 24 O08 7 1 01-95 | 25 518 8 
2 05-24 | 22 318 17 @ 03-62 | 28 414 60 (30) 
3 13-82 | 27 148 23 3 08-43 | 22 068 19+ 
4 25-53 | 23 077 15(5) 4 19-03 | 20 318 19:(8) 
5 27-23 | 26 087 14 5 25-54 | 23 217 9 
6 30-78 te 27 237 8 6 28-64 | 20 55? 8 
7 37-23 | 28 465 9 7 40-94 | 24 168 13 
8 37-65 | 24 148 8 8 41-14 | 27 286 fi 
9 38-14 | 21 407 14 (9) 9 41-93 | 21 086 10 
10 45-42 | 29 265 15 10 43-23 | 26 275 12 
11 AGB 7£a) eet asi? 9 iil 49-73 | 21 208 8 
12 49-72 | 27 575 536) 12 52-94 | 22 478 9 
13 52-14 | 25 457 11; 
14 58-74 | 24 388 9 06 
1 02-84 | 28 568 7 
04 2 04-62 | 20 135 23 
1 06-33 | 24 155 14 3 15-93 | 28 167 9 
2 06-73 | 22 526 10 4 19-14 | 24 088 9 
3 12-92 | 29 405 17 5 19-43 | 27 216 g(9) 
4 13-73 | 21 04? 26 6 20-13 | 26 057 7 
5 14-73 | 22 078 8 7 21-03 | 25 187 11+ 
6 20-34 | 26 334 11 8 24-54 | 20 357 11 
7 22-54 | 20 048 Fi 9 30-44 | 27 168 12 
8 26-74 | 27 038 15+ 10 34-71 | 20 374 85 (67) 
9 27-93 | 28 075 12 11 37-43 | 27 367 17 
10 30-45 | 23 347 7 12 37-73 | 28 408 17 
11 30-73 | 29 498 13 13 51-03 | 21 208 9 
12 30-94 | 24 447 8 14 51-24 | 22 327 8 
13 32-73 | 25 338 8 15 55:53 | 26 407 16} 
14 34:65 | 20 458 5 16 56:91 | 24 134 5940) 
15 34-83 | 22 328 10 
16 36-75 | 27 218 8 
17 37-15 | 24 398 8 | 
18 42-71 | 28 188 827) | 
19 45-93 | 20 405 19 | 
20 50-43 | 28 367 12 | 
21 52-93 | 22 045 38 (21) / 
22 53-3? | 20 364 18 | 
23 56-54 93 43" 7 | 
{ 


(5) Possible confusion with IAU 03S3A side lobe. 
(6) 40”, 

(2) 35”. 

(8) Possible confusion with TAU 05S4A side lobe. 
(®) (NGC 2217). 

(10) 45”, 
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TaBLe | (Continued) 


Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. Density 
No. R.A. Dec. (Ose! No. R.A. Dec. (10528 
Sp W m-? (e/s)-1) s. W m-? (e/s)~7) 
h m em hema: ee 4 
07 09 
al 04-28 22 588 6 | 05-74 29 218 9 
2 06-98 29 14? 45 (22)0)) 2 20-48 27 24? 8 
3 09-0? 20 344 71 (33)00 3 215 21 058 ri 
4 15-48 24 578 17 4 24-54 25 bi 10 
5 16-18 26 288 18+ 5 26-24 29 48° 13 
6 2018 S25 oe in 6 30-78 20 077 12 
7 21-98 23 (287 15+ 7 35-34 28 508 16:7 
8 24-78 24 516 10 8 36-58 25 408 12 
9 25-58 20 398 10 9 42-75 27 DOE 8 
10 25-84 28 457 15+ 10 47-4? 24 598 19 
a 27-28 22 06° 167 TH 512105 23 48° 9 
12 39-24 24 188 38 (24) 12 56-08 28 558 30 (18)8) 
13 40-78 28 158 13 13 56-68 20 258 6 
14 50-54 21 478 8 14 59-68 23 388 15 
15 50-94 26 207 13 
16 54-44 23 288 7 10 
1 02-8? Zl y3® 4804 
08 2 06-54 30 00° 10 
i 12-38 23 338 7 3 08-54 28 308 gee) 
2 12-98 24 487 10 4 10-98 27 508 12 
3 18-33 29 46° 16 5 12-58 23 388 10 
4 25°12 20 155 26 6 24-68 29°46 Bilis 349 
5 25-14 23 137 11 7 30-88 23 37° 18} 
6 26-75 21 308 8T 8 35-54 24 508 8 
7 33-14 24 498 Ta 9 35-94 28 587 10 
8 34-58 26 007 9 10 35-98 26 105 15 
9 34-88 27 O18 16 11 40-14 23 20° 12 
10 35°74 22 197 it) eve 42-78 28 467 147 
11 42-75 26 308 14+ 13 43-04 22 107 OG) 
12 43-04 23 30° 15 14 48-94 20 197 10 
13 43-48 29535" 1202) 15 53-23 26 007 ll 
14 43-94 22 S8 9 16 53-54 27 438 13+ 
15 50-84 22 527 8 
16 50-98 20 366 19 
7 55-33 21 187 a 
18 56-48 27 188 13 
19 59-6? 25 495 5418) 


{1) Perhaps several sources. 

“2) Perhaps one extended source with 08-37. 
(13) 32”. 

(14) —1Q”, 

5) A doubtful source. 

46) Perhaps a background irregularity. 
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TaBLE 1 (Continued) 


Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. Density 
No. R.A. Dec. (Gat No. R.A. Dec. ons 

S. W imn-2 (e/s)) Ss. W m-? (c/s)-) 
h m oe hm at wilt 
Uy 13 
1 03-2 24 337 | g16) v4 45-94 95 147 17 
2 03-35 20 525 16 8 50-74 BY tek 9 
3 08-18 22 51s | 16 9 55-74 23 358 9 
4 T5298 25 557 9 10 56-73 21 55® 8T 
5 27-15 28 53° 19+ 
6 29-65 26 567 8 14 
7 atefon le 26 184 28 i 03-8? 27 18° 20 
8 39-0? 28 30° | 2 08 - 08 22 27" 8 
9 49-8? 30 00° 7 3 08-85 27 29% 13 
10 Sli 28 088 9 4 09-8? 20 188 12 
TT || 52-138 Zo) 11? 13 5 10-92 23 38? 19 
12 SO 22 O78 15 6 14-18 21 245 il 
7 19-5? 29 188 15 
12 to} 20-28 27 145 40 (20) 
1 08-48 25-028 14} 9 20-6% 24 52° 14 
2 09-8> 26° 5610 7 10 PAT fe: 29 426 26 
3 18-58 27 407 10 itu 22-04 23 487 iyi 
4 19-24 26 467 10+ 12 PAS fe 21 4710 18 (11)0) 
a) 2214 25 515 15t 16 39-08 26 28° 24 
6 Hagin 20 528 12 14 39" 0* 24 38° us 
Uf 32-9? 24 574 28 15 54-78 28 21° i 
8 SH CH 20) 107 9 
9 38°3° lo 94 15 
10 40-93 PATE athe Vit il 03°88 22 207 10 
algh 45-04 28 588 18 2 05-17 26 278 9 
12 51-7? 29 034 4737) 3 05:84 29 23° 35 (20) 
13 55-84 25 257 8 4 07-08 20 138 9 
14 56-6 23 08° 27: 5 07-9° Zo Eb" 7(16) 
15 po 02 20 09° 16 6 iors 22 as Oren) 
L305 28 au 
13 } 25-24 26 42° il 
1 02-73 | 20 538 11 9 26-74 | 27 38° 18+ 
2 03-5 25 008 15 10 26+ 87 20 108 8 
3 09-0? 22 004 6108) hl 44-34 22 42° 3 
4 12-24 24 587 9 12 50-0 25 10% 12 
5 34-73 29 326 3619) 13 56-28 21 35° 36 (2) 
6 36-93 sual O}s%" 14 14 56-48 24 058 18 
(a BTR, 
(18) 15”. 
(9) NGC 5236. 


(20) <12”, 
(21) Superimposed on background irregularity. 
(22) > 30°. 


(25) Kepler’s Supernova—size > 45”. 


®°) TAU 17S52A, complex distribution—only the integrated flux is given. 


SSSI 
(28) ~30". 
(29) ~ 60”. 
(30) ~ 40”. 
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TABLE 1 (Continued) 
Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. Density 
No. R.A. Dec. (10=-© No. R.A Dec (10-45 
8. W m-? (e/s)-1) Ss. W m- (e/s)-*) 
h m eS o h m s 3 
16 18 
il 01-98 2Srpos 35 a 02:71 a Pere 32 2539) 
2 05-04 20) 2EE 20 2 02-84 26 037 1945) 
3 18-25 Domoon 14 3 04-0? 20 104 70 
4 19-28 28 008 18 4 05-8? 27 00° 34 
5 28-78 26 508 27 5 12-88 24 095 PIF (718s 28) 
6 35-98 23 388 10@5) 6 14-95 22 40° 24 
vf 36-64 21 33? 8 7 16-2° 25 458 307 
8 39-64 25 251° 19 8 22°58 20 057 70 (27) 
9 43-48 22 225 22 if) 23-88 24 037 31 
10 53-3? 25 408 14 10 24-83 29 10° 26 
11 53-63 20ets* 11 Il 33-68 23 597 39 
12 593% 28 437 25>) 12 48-58 25 576 2316) 
13 50-95 26 527 22,46) 
17 14 53-04 22 27° 19@6) 
1 OL: 77 24 24? 15) 
2 08-28 20 568 19 19 
3 09-5? 23 1p> 4223) 1 00-18 Zo O° 35T 
4 09-98 28 00? 24 Zz, 04-08 25 095 17 
5 10-8? 25 005 47 (24) 3 12-8? 26 585 537 
6 13-24 29 50? 38 4 13° 6® 24 498 35 
7 Irons} 28 568 32(25) 5 16-48 21 417 18 
§ 22-16 20 198 17 6 20-48 PAE wis 8 
9 22-44 26 50 18 7 25-63 25 445 21; 
10 23-05 28 05° 1479) 8 29-24 26 408 42 
iil 21-81 21 298 110!25) 9 33-48 24 08° 9 
12 37-04 21 337 21 10 34-98 22 228 11 
13 43-79 28 4122 | 450026) 11 42-6* | 27 05° 12 
14 52-24 2) 30? 194 (87) 12 46-44 28 O1§ 18 
15 52-44 23 05° 82+ 13 50-44 21 476 11 
16 57-44 23 268 900°) 
17 59-98 28 008 35 (28) 
SoU 
(24) > 40", 
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Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. Density 
No. R.A. Dec (10-25 No. R.A Dee ((UO=e6 
S. W m2 (ce/s)) 8. W m~ (e/s)-1) 
hm say h m ae Te 
20 22 
t 02-83 25 03° 7185) il 14-74 28 238 23 (14) 
Zz 15=1¢ 29 407 1] 2 16-72 20 507 ll 
3 1834 22 23" 8 3 22-84 27 328 9 
4 21-14 21 115 ll 4 26-58 29 307 9 
5 24-34 20 346 9 3 32-74 23 028 9 
6 24-54 21 445 18 6 45-94 24 238 9 
7 27-48 29 407 gas) 7 46-54 26 537 8 
8 2eos 23 005 22 8 47-38 23 195 13 
9 34-83 24 188 6 9 48-14 22 176 7 
10 35-78 20 16? 9 10 55-33 23 01° 9 
11 38-83 27 578 12 11 59-48 26 556 8 
12 40-63 26 43° 18+ 
13 45-28 24 167 12 23 
14 53-08 20 07® 25 iD 03-88 25 268 13 
15 58-7? 28 134 59; (31) 2 07-14 28 20° 7 
16 59-94 21057 8116) 3 10-94 Zilaone 8 
4 nL 27 388 23 
21 5 17-58 22 138 13 
1 04-61 25 398 100'32) 6 18-18 24 255 12+ 
2 06 - 23 23 55° 21f if 22-98 23 208 12 
3 1362 mae aS 24 8 26-48 21 255 20 
4 15-33 25 078 15 g 25% O45) Wed 27 (16) 
5 18-93 26 456 17 10 42-58 24 105 10 
6 23414 23° 58? 9 ate 45-83 28 208 15) 
7 24-98 29 247 10 12 46-78 26 257 6 
8 26-63 22 557 9 13 47-98 25 164 16 
g 28-38 20 54% 16 14 49-88 23 198 7 
10 34-78 20757 11 15 51-98 22 258 7 
11 36-33 25 598 15 16 55-44 21 307 1] 
12 36-78 27 58* 8 
13 50-04 20 037 15 
14 50-78 28 108 12 
15 59-9? 28 36° 10 
(31) <20”. 


82) > 50". 
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TABLE 2 
SOURCES BETWEEN DECLINATIONS —30° AND —40° 
Sources which may be ‘‘ extended ’’, that is, resolvable, are indicated by a dagger. A colon has 
been placed beside uncertain flux densities. Angular sizes, where available, are indicated in the 
footnotes. Details in Section II 


Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. Density 
No. R.A. Dec. (te No: | Reza. 7] Dec. (10-26 
S. W m= (e/s)—) 8. W ma-*(e/s)—) 
h m ee h m pes / 
00 01 
vi 00-2 Syl oS} 17 10 30-45 38 31® | 10 
2 02-74 | 36 508 8 1 31-62 | 36 443 Ets) 
3 06-98 on Le 14+ 2 41-28 39 485 | 15'6) : 
4 09-84 Sooe 12 13 44-38 oy fies 2 ee 9 a 
S 12-98 38 275 19 Id 54-28 36 235 16 A 
6 14-18 | 31 284 17 15 57-93 | 31 108 260) ; 
7 14-5% | 36 418 i : 
8 Zoot Sey OK 20 02 
9 25°3* 39 21° 8 7 04-53 | 37 528 10 
10 32-95 | 33 417 9+ 2 10-08 | 34315 | 13 
Il 33-0" 37 478 9 3 16-0? | 36 455 297(8) 
IZ 34-38 30 417 9 4 18-88 / 38 216 19 
13 36-2? 39 245 30) 5 24-5? 30 51° 14 
14 Spode 33 067 ll 6 25-08 34 05° 9 
Le 41-9? 35 384 13 7 38-54 oi 13° 8 
16 51-68 36 38> 15 8 45-38 3b 34? 10 
17 63-55 38 035 12@) 9 59-52 34 368 137 
OL 03 
i 00-68 31 365 12 v4 20-6? 37 23* 950 (825) 
2 03-38 38 25° 9 2 32-64 39 095 | 10 
3 07-28 | 37 028 8 3 36-22 | 35 355. | 9ga0 
4 07-38 | 39 208 10 4 36-73 | 32 O14 14 
5 07-93 | 35 O18 30 (15) 5 42-809 “S) Ver Tay 
6 12-24 32 307 17(4) 6 44-92 34 354 330) 
Hy 12-75 ol 218 10 ‘4 45-74 35 187 8 
8 19-98 37 495 19 8 65-18 30 215 12 
i) 24-78 36 417 10 9 57-14 31 025 14 
qd) i Baye 
@) NGC 300. 


3) Perhaps several sources. 

‘) Sources 01-36, 01-37 perhaps form one extended source. 
(5) S45”. 

‘6) A doubtful source. 

(ORNS 

(8) >40", 

‘) Fornax A; IAU 0383A. NGC 1316 (NGC 1317). 

0) NGC 1399 (NGC 1404), 

(11) S452 
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TABLE 2 (Continued) 


Position (1950) Flux Position (1950) Flux 
Ref. | _ —_______ Density ——_— Density 
No. R.A. Dec. (10-28 R.A. Dec. (10-26 
Ss. W m-? (e/s)-*) 8. W m-? (e/s)-1) 
hm eae hm o 4 
04 06 
ii 05-33 32 067 11 1 01-1? 34 274 22 
2 06-3? Ble ite 127 2 02-23 32 208 LEB 
3 10-6? 34 334 16 3B 02-43 35 248 9 
4 19-85 31 018 6 4 07-84 30 467 8 
5 20-5? 33 265 12t 5 14-34 390327 8 
6 27-12 36 384 35122) 6 14-34 35 068 7 
7 31-13 39-57¢ 14 7 17-9? anf ike 18 
8 37-83 31 098 10 8 25-52 35 264 26 
9 38-54 36-11 8 9 30-74 36 388 137 
10 41-34 34 458 6 10 31-14 39 16° 8 
os 43-8? 33 00° 8 11 40-73 36 058 7 
12 45-05 39 128 18 hy 46-7? 39 564 26 
13 46-48 35° 53° 10 13 50-23 34 198 9 
14 54-5? 30 16° 78 (43)0%) 14 50-23 32 408 9 
15 55-48 31 587 9 15 58-84 34 348 5 
05 07 
1 03-84 38 557 8 1 00-08 ai 337 14} 
2 04-64 37 51? 13 2 00-94 33 3410 11 
3 04-64 32 218 6114) 3 06-04 39 588 12 
4 04-9? BOL 7 (6) 4 07-63 35 538 15 
5 11-5? 30 34° 29 5 15-3? 36 275 187 
6 21-48 36 246 66:15) 6 17-24 31 408 9 
7 22-98 32 488 18 7 18-9? 34 157 1807) 
8 24-98 30) 55? 8 8 19-03 35 018 ga) 
9 34-23 31.237 9 g 25-88 35 407 9 
10 38-13 33 401° 30 (16) 10 27-68 30 108 8 
11 41-9% | 38 428 8 11 29-74 | 32 037 14 
12 45-84 | 30 018 7 12 35-88 37 595 2118) 
13 46-33 33 04° 1416) 13 36-0? 30) 21% 19 
14 50-03 35 06% 13 14 37-93 36 207 lly 
15 52-04 36a312 10‘) 15 49-44 38 458 24 (13) 
16 54-33 32 147 14 16 50-34 35 068 8y 
17 58-25 AE UE 9 Sug 51-8 31 068 gi19) 
9-3? 38 407 16 18 53-84 35 408 7 
E e -04 30 218 15119) 
19 59-8? 39 495 20 19 55-0 
(12) < 10”, 


43) Perhaps two sources. 

(14) (NGC 1800). 

(15) 20”, 

46) Sources 05-310, 05-313 perhaps form one extended source. 

(47) Sources 07-37, 07-38 perhaps form one extended source. 

: 48) Perhaps extended or superimposed on galactic emission feature. 
49) Sources 07-317, 07-319 perhaps form one extended source. 
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TABLE 2 (Continued) 


Position (1950) Flux Position (1950) Flux 
Ref. | ,—___ Density Ref. |————_,-—___ Density 
No. R.A. Dec. (Gee No. R.A. Dec. (10,28 
8. W m-? (c/s)-1) Ss. W m-? (c/s)-*) 
h m ° y h m ° , 
08 10 
1 07-6? 38 505 38 (19) 12 51-78 34 038 127 
2 25-54 31 427 19 (20) 13 52-38 37 265 18 (12) 
3 27-28 30 314 20 14 56-3? 35 525 12 
4 28-55 32 517 30 (20)20) 15 59-58 31 007 10 
5 32-3? 34 216 10 poses see 
6 40-88 39 428 17 1l 
7 43-08 30 398 12(22) 1 00-38 32 108 8 
8 43-08 33 39° 7, 2 23-04 37 008 7 
9 45-44 35 437 8+ 3 23-58 35 128 17 
10 48-18 34 388 9+ 4 25-88 32185 107 
11 50-5 33 407 ll 5 27-13 34 207 11 
12 51-84 31 508 9 6 35-14 39 02° 1l 
13 59-84 34 308 6 7 35-54 36 558 8 
————_|—____|_———_|__ 8 36-08 31 588 28 
09 9 43-44 35 20° 9 (28) 
1 01-98 36 487 15 10 43-5? 31 416 27 
2 02-68 38 346 25+ 11 44-18 33 14§ 16 
3 21-08 32 11% 9: 2 46-58 37 587 1l 
4 21-12 39 565 9 13 48-74 35 408 8 (24) 
5 34:74 32 258 14 14 50-8? 34 498 104) 
6 38-88 39 168 8 15 55-58 31 308 13 
7 41-58 30125? 8 ss 
8 42-38 35 407 8 12 
9 44-04 33 008 8 1 01-28 35 376 13+ 
10 45-3? 36 618 ll 2 05-48 33 437 17 
th 50-08 38 3510 25 (12) 3 18-48 36 567 14+ 
12 51-48 39 267 17 4 32°38 33 287 227 
5 40-23 38 408 lly 
10 6 41-85 36 001° 6 
1 02-28 36 58? O(2s) 7 43-28 35 001° 7 
2 03-54 32 308 14 8 54-28 30 095 28 
3 11-54 31 458 10 9 56-28 33 148 18 
4 14-38 30 088 11 10 59-68 36 546 14 
5 16-54 331? 14 
6 Irie 32 488 19} 13 
7 29)°3? 35 588 25 (15) 1 05-35 30 037 12 
8 30-8? 34 098 20 2 09-04 36 617 9 
9 35-24 39 50? 6 3 33-9? 33 425 70(25) 
10 43-98 38 008 7 4 46-94 39 007 29 
11 45-18 35 545 9 5 59-14 32 458 8 (6) 


0) Perhaps part of 08-34. 

@1) Hxtended parallel to galactic plane. 

(22) Perhaps one extended source with 08-273, 
®3) Perhaps a background irregularity. 


4) Sources 11-39, 11-313, 11-314 perhaps form one extended source. 
(25) +50”, 14296. 


(26) +45” (NGC 5419). 


COSMIC RADIO SOURCES BETWEEN DECLINATIONS —20° AND —50° 687 
TABLE 2 (Continued) 
Position (1950) Flux Position (1950) Flux 
Ref. Density Density 
No. R.A Dec. (10-26 R.A. Dee. (10-26 
Ss. W m-? (e/s)-) Ss. W m-? (c/s)-) 
hm ae hm Oe 1 ¥ 
14 17 
1 00-6? 36 168 20 26-23 34 455 60 
2 01-53 33 53* 5726) 26-7? 34 024 120 
3 18-83 30 557 17 28°37 32 435 95 
4 21-73 38 107 20 36-6? 30 50° 165 
5 28-13 31 53° 16 39-93 36 018 31t 
6 29-63 34 546 12(6) 
vf 40-33 30 288 147 18 
8 5i-37 36 224 41 1 03-73 36 208 9 
9 51-43 34 107 19 2 06-93 39 525 9 
3 17:8? 39 115 4] (28) 
15 4 26-13 38 115 18 
1 04-23 a9 226 14 5 oD? 2 Sill alg 19 
2 07-48 30 307 gs) 6 aM lial 36 178 12(29) 
3 17-68 39 36° 10 if 48-63 35020" 1b 
4 22-48 at 20° 14 8 49-48 34 414 17(6) 
5 30-83 ol 245 12(23) 9 58-64 38 248 12 
6 40-98 OOo sie 21:03) 
if 53-84 35 567 10 19 
12-88 37 037 18 
16 16-9? 35 448 16 
1 04-94 32 467 19(6) 27-03 36 06° 16 
2 10-93 39 246 56 (28) 37-18 36 115 13 
3 17-53 36 408 13 56:4 35 494 45 (30) 
4 20-0? 34 308 13 
5 22-48 32 058 24 20 
6 22 -5* 30 367 297 1 04-54 39 008 177) 
if 26-65 35 097 ll 2 13-5? 36 057 be 
8 57-88 33 578 19 3 13-94 31 44° 10 
4 14-08 31 007 157 
17 5 24-44 | 32 O15 18 
il 02-13 31 548 32 6 29-83 37 408 7 
Zz 06-64 37 10° 1007 ff 32-3? 35 045 41 (32) 
8 11-12 38 234 500 (300) 8 48-03 37 O78 197 
4 19-83 35 405 75 (27) i) 56-23 33 007 9 
5 21-85 38 147 TAS 
i 


(27) > 45”. 

(28) —90”, 

(29) Superimposed on background irregularity. 

(30) Flux measurement uncertain due to Cygnus A side lobe ~20’. 
(31) A doubtful source, perhaps extended or a background irregularity. 


82) <10", 
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TABLE 2 (Continued) 


Position (1950) Flux Position (1950) Flux 
Ref. |_———_—_—,—_—__ Density Ref. |————__—— Density 
No. R.A. Dec. (10-26 No. R.A. Dec. (ites 
8. W m-? (c/s)-*) 8. W m-? (c/s)-*) 

h m ° i ° ¢ 

21 
1 00-33 B9naT" oo Le 9 
2 O79 34 047 37 335 21 
3 10-44 SomlOs 
4 14-68 30 36° 
5 16-04 33 35? 32 267 13 
6 23°18 31 225 34 317 9 
if 45-58 30) 238 32 566 16 
8 57-34 38 04° 34 586 24 

37 38° FS) 

22 38 396 11 
al Iifsyoa les 30 45? 34 594 39(83) 
2 24-88 30 465 36 085 18 
8 43-48 30 425 
(33) > 40”. 


EE es 
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TABLE 3 


SOURCES BETWEEN DECLINATIONS —40° anpD —50° 


Sources which may be “‘ extended ”’, that is, resolvable, are indicated by a dagger. A colon has 


been 


placed beside uncertain flux densities. Angular sizes, where available, are indicated in the 
footnotes. Details in Section II 


Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. |——__ —_____ Density 
No. R.A. Dec. (10-26 No. R.A. Dec. (10-76 
S. W m-? (e/s)-1) 8. W m-~? (e/s)-1) 
hm oe a hm OF 
00 02 
1 01-68 48 238 10 1 01-98 44 025 9 
Z 03-83 42 508 17 2 07-24 42 005 16 
3B 08-2? 44 404 60 (31) S 14-08 48 035 65 (31) 
4 28 41 195 13 4 19-72 45 335 a; 
5 17-83 49 278 9 5 22-63 42 037 17 
6 19-18 47 175 13 6 24-24 43 158 7 
7 32-74 Ab 338 7 27-48 40 147 24 (14)) 
8 34-28 49 128 8 30:5? 41 225 i 
9 36-64 41 167 9 34-63 43 498 8 
10 39-8? 44 374 10 40-28 42 03? 12 
11 43-7? 42 254 itl 43-58 45 075 10 
12 44-04 49 578 iy 55-08 44 057 7 
13 48-84 44 467 13 55-34 48 495 16; 
14 50-14 43 235 
15 52-98 49 278 03 
1 03-5? 46 427 8 
Ol y 15:68 44 008 10 
1 03-2? 45 225 A41(5) 3 18-83 45 207 19°) 
2 03-6 41 294 25 (17) 4 35°5* 41 285 14 
3 07-44 48 128 5 7) 38-48 40 378 10°) 
4 12-94 44 468 8 6 38-53 47 465 13 
5 14-08 47 345 346) yi 42-0? 44 208 11 
6 15-73 | 41 066 g(2) 8 43-64 | 45 467 9 
i 24-34 43 257 8 9 51-28 4] 235 9 
8 Done AL 7s 33T 10 53°08 43 048 8 
9 31-48 44 568 18 i 54-18 48 244 14 
10 39-53 44 407 9 
11 39-84 45 518 10 
12 40-54 43 587 14 
13 52-38 43 46° 10 


(4) Perhaps two sources. 

(2) A doubtful source. 

<5", 

(4) +40", 

(5) 93”, 

(6) > 40”, 

(7) > 40”. 

(8) haps several sources. 

(9) neolnce subject to side-lobe influence of IAU 03S3A. 
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TABLE 3 (Continued) 


Position (1950) Flux 
Ref. |——__, —_—__- Density 
No R.A. Dec (10-8 
S. W m-? (e/s)-1) 
h m ° / 
04 
1 08-24 47 068 
2 13-98 42 446 
3 18-18 40 086 
4 18-28 41 358 
5 25-18 49 06° 
6 34-28 49 507 
7 36-5 47 368 
8 37-88 45 008 
9 38-28 43 315 
10 55-48 40 30° 
I1 55-68 43 068 
12 55-63 46 207 
05 
4 02-05 41 50? 
2 11-94 48 325 4] an 
3 18-31 45 488 57002) 
4 25-83 40 567 14 
i) 31-58 45 408 19°) 
6 34-88 49 468 16 
7 45-28 48 175 15 
8 46-28 44 358 13 
9 46-24 45 576 17 
10 47-8? 40 524 317.48) 
if] 47-84 42 038 7 
06 
1 06-13 | 49 387 15+ 
2 07-34 42 278 9 
3 12-28 47 188 27 
4 16-63 48 445 9 
5 19-88 45 078 12 
6 20-14 47 007 7 


CNR DAR & WMH 


R.A. 


h 


06 


08 


Position (1950) 


m 


57: 


0) Perhaps extended or several sources. 


an <920”, 
{2) Pictor A; IAU 0584A 55”. 
(13) ~ 39”, 


{4 Sources 07-412, 07-414 perhaps form one extended source. 


5) Puppis A; TAU O8S4A. 
(16) 55”, 


385 
257 
257 
268 


Flux 
Density 
( 1 Q-26 
W m-? (e/s)~*) 


12 
97 (14) 
12 
184) 


27 
10 
10®@ 
690 (514)"5) 
1100 (276)%) 
10+ 
26+ 
17 
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TABLE 3 (Continued) 


Position (1950) Flux Position (1950) Flux 
Ref, | _—________ Density Ref. |———__—______ Density 
No. R.A. Dec. (10>** No. R.A. Dec. (10-26 
8. W m- (e/s)-) 8. W m-? (c/s)-1) 
h m ad hm 955 2 
09 12 
1 06-94 46 378 7 3 16-24 46 O17 197 
2 09-72 45 078 ll 4 33-74 41 257 11 
3 12-6? 43 347 gti) 5 46-88 40 565 4518) 
4 27-04 46 428 8 6 51-9? 47 435 14: 
5 27-85 45 147 7 7 55°78 41 338 16 
6 43-24 49 537 9 
7 46-48 41 227 127 13 
8 49-3° 46 508 15) 1 02-65 49 077 2009) 
9 51-84 45 167 15 2 22-4? 42 414 8700 (2170) 2) 
10 51-98 43 278 12 3 38-33 40 327 9 
4 47-14 45 367 100 (32) 0) 
10 5 56-0? 41 384 35: 
1 02-94 49 446 OF 
2 03-3? 46 275 8 14 
3 14-94 44 257 9 1 08-98 41 33° 14 
d 17-93 42 268 5140 2 12-68 43 238 387) 
5 19-14 45 487 10 3 14-44 47 237 14 
6 19-24 49 227 10 4 17-64 49 447 167 
vi 24-74 40 57° 13 5 25-28 43 075 90: (30)? 
8 25-84 47 528 6 25-78 47 488 29 
g 26-43 45 007 7 32-94 45 05° 10 
10 32-25 41 178 8 42-48 42 138 15 
iit 39°34 46 05° 9 45-03 46 57° 33 
12 46-44 49 157 10 45-14 48 098 10 
13 52-24 40 007 11 47-08 40 18° 14 
12 50-08 41 007 9 
11 13 61-3? 42 178 19) 
1 04-43 41 56° 14 57-5? 48 04° 17+ 
2 06-7? 48 237 15 59-9? 41 425 5b: 
3 16-13 43 247 rl 
d 23°07 48 13° 15 
5 33°9° 43 157 a 04-48 42 547 14} 
6 43-33 48 208 2 07-64 43 568 14 
2 49-18 44 567 3 27-4? 42 21 100 
8 55-58 42 167 4 2-45 49 448 21+ 
5 38-78 40 427 18 
12 6 41-23 48 047 47 (25) 
rf 10-74 41 437 y, 45-94 45 04° 24°) 
2 11-53 44 568 8 | 55°3° 46 127 507 (23) 
| 


7) 30” TAU 1084A. 

(18) <15” (NGC 4696). 

(9) NGC 4945. 

(20) Centaurus A; IAU 13S4A. NGC 5128. 

(21) Possibly a galactic feature. 

(22) Perhaps several sources, the interpretation is doubtful. 


23) + 30”. 
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TABLE 3 (Continued) 
Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. |——————_ Density 
No. R.A. Dec. (10-26 No. R.A. Dec. (10-26 
S. W m- (e/s)-) S. W m-? (c/s)-) 
h m as hm cs 
16 19 
il 03-08 44 258 19@) 3 21-78 43 01° 17 
2 17-44 45 51? 26; (24) 4 25-28 41 237 17 
3 19)59% 49 287 76 (25) 5 31-88 49 535 9 
4 24-94 48 407 64 6 32-6? 46 25° 141 (28) 
By) 25-14 42 108 307 7 32-98 48 31’ 8 
6 33-2 40 327 25 (2) 8 35-48 42 347 8 
di 36-35 46 314 3306) 9 36-74 48 01° u 
8 51-18 44 00° 557 10 40-48 40 41° 387) 
9 52-14 42 258 32 td 45-68 47 085 11 
10 57-65 46 228 30: 12 shou) 49 596 14 
11 59-04 41 185 90 (70) 13 53-08 42 385 31t 
17 20 
1 03-88 44 246 25 1 03-58 47 40° 9 
2 10-48 45 367 35 2 05-88 42 378 157 
3 14-48 44 007 18+ 3 65 42 238 9 
4 14-88 48 207 25 f 20-54 49 537 8 
By} 16-6° 49 337 29 i) 27-38 41 30° 22 
6 folk 46 478 18 6 28-98 45 234 10 
7 31-33 48 318 127 7 37-38 47 528 1l 
8 54-48 41 337 30f 8 43-94 43 378 16 
9 47-05 44 485 13 
18 10 51-48 48 228 12 
il 04-65 45 28° 77 (45) 
2 12-85 48 207 16: 21 
3 39-8? 48 365 4127) 1 09-44 43 128 15 
4 40-93 | 40 277 22 2 10:24 | 45 268 7 
5 43-93 | 49 497 12 3 16-94 | 47 508 7 
6 44-65 | 43 067 8 4 22-34 | 44 56? 7 
7 53-24 | 43 008 7 5 27-78 | 48 258 10 
8 59-63 | 41 255 13 6 37-15 | 44 264 10 
7 40-32 | 43 214 37(30) 
19 8 50-53 | 46 39° 14 
i 12-08 46 28° ll 
2 Higjot 45 395 13 


4) Possibly a background irregularity. 


(25) 45”, 
(26) + 45”, 
i 30m 
(28) 19”, 

(29) ~ 30”, 
(30) ~ 40”, 
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TABLE 3 (Continued) 


Position (1950) Flux Position (1950) Flux 
Ref. Density Ref. Density 
No. R.A. Dec. (IE: No. R.A. Dec. (ons 
8. Witms(C/s)is2) 8. W m-? (e/s)-) 
hm ae hm ae 
22 23 
1 06-8* 49 576 8 1 07-88 48 427 iil 
2 06- 8° 43 598 13 z 10-33 44 437 10 
3 25-38 40 56° 28 3 24-28 40 345 15 
4 35-05 43 448 6 4 31-9? 41 424 5032) 
5 44-14 44 267 6: 5 39-08 44 507 6 
6 50-1? 41 105 42(31) 6 42-98 47 258 © 
vi 56-88 41 09° 14 7 49-75 43 248 8 
8 57-08 46 09% 10 8 | 51-78 45 267 fe) 
| 
(31) 30”, 
{330267 


although we had available a number of ‘‘ non-scanning ” records of high sensitivity 
not available to Rishbeth. His catalogue is valuable in that it is based on the 
contour diagram, gives information about the shapes of some extended sources, 
and includes several large concentrations not included in our cataloguing scheme. 
There are, however, often large discordances in the flux densities of weaker 
sources. 


ITI. IDENTIFICATIONS 

A search for identifications with bright optical objects has been carried out 
as in papers I and II; the atlas of Becvar (1951) has again been used in this 
search. Some of the area, north of declination —33°, is included in the National 
Geographic Society-Palomar Observatory Sky Atlas. The previous catalogue 
and a small portion of the present one has been compared with this Atlas in a 
systematic way (Mills 1960) and no doubt a similar comparison of the remaining 
common area will prove interesting, but it is outside the scope of the present 
paper. 

In contrast to the earlier work, no Hm regions can be unambiguously 
identified in the present catalogue ; this arises partly as a result of their uneven 
distribution and partly because the zone of the catalogue crosses the galactic 
plane, in one of the two places, close to the galactic centre. Here the background 
temperature is of the same order or higher than the electron temperatures in 
the nebulae so that, if they are observed at all, they appear as “ holes” in the 
background emission: some of these are listed elsewhere (Mills 1959). Co- 
incidences with planetary nebulae and globular clusters are again not significant : 
there is a moderately close coincidence between the radio source 17 —2/ and the 
globular cluster NGC 6284, but this would be expected by chance. Two other 
identifications with galactic objects are well known. These are Kepler’s 
supernova, identified with 17—21/ (Mills, Little, and Sheridan 1956), and 
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08 —44 (Puppis A) identified by Baade and Minkowski (1954) with a galactic 
nebulosity of a peculiar filamentary type. 

There are numerous possible identifications with bright galaxies in the 
catalogue. Those listed include all coincidences within 1™ in Right Ascension 
and 20’ in declination ; there are 14 such coincidences and the chance expectation, 
based on the numbers of catalogued galaxies and radio sources in the area, is 
about 6. To reduce effects of chance we restrict attention to these coincidences 
within 0™-7 in Right Ascension and 13’ in declination, as in paper II (these 
numbers represent approximately twice the mean probable error in each 
coordinate). There are 10 such coincidences with 3 expected by chance, but 
since 5 of these were previously known the new information is not great. They 
are listed in Table 4. 


TABLE 4 
POSSIBLE IDENTIFICATIONS WITH BRIGHT GALAXIES 
: Galax 

Beale —— My.9-M » Notes 

Source NaC Type 
00—221 253 Se 1-2 
03—31 1316 Sap —4:6 Fornax A 
03—31 1317 SBb 8a) Probably a chance coincidence 
03—33 1399 Eo —1:6 
06—25 2217 SBb —0-5 
12—46 4696 El —3-1 The angular size of the radio source is 

small. Probably a chance coincidence 

13—41 4945 SBe 1-3 
13—42 5128 Ep —3-9 Centaurus A, a well-known identification 
13—33 14296 E —4-1 The radio size supports the identification 
13—25 * 5236 Se 0-9 


Of the previously known systems we have three “normal” galaxies 
NGC 253, NGC 4945, and NGC 5236, and two ‘“ radio’ galaxies NGC 1316 and 
NGC 5128. Two of the coincidences can almost certainly be ascribed to chance, 
NGC 1317 and NGC 4696. In the former, the radio source is, by virtue of its 
position and size, almost certainly associated with the neighbouring galaxy 
NGC 1316 ; it would be a remarkable coincidence if NGC 1317 also were con- 
tributing to the radio emission. In the case of NGC 4696, which has been noted 
also by Basinki, Bok, and Gottlieb (1959), the angular size, <15” arc, appears too 
small. There are numerous fainter galaxies in the vicinity; but it appears 
most likely that the radio source is considerably more distant than any of those 
visible on the plate of Basinki, Bok, and Gottlieb. The galaxies NGC 1316 and 
NGC 5128 have been discussed extensively in the literature and will not be 
considered further in this paper (for the most recent results see Sheridan (1958) 
and Wade (1959)). 

It is interesting to see whether the above coincidences might be expected 
assuming a radio “ luminosity function ” of the form Suggested by Mills (1960). 
In this, the probability of a galaxy, chosen at random, having a ratio of 
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radio to optical emission, defined by My.9—M,, between m—4 and m+} 
(where m=m,.,—m,) is given by 


log P=1:0—0-42m. 


Applying this result to the galaxies listed in the catalogue we have used, we 
would expect 6 galaxies with m,.,—m p<—2 and 1 galaxy with m,..—m,< —4. 
In fact there are 5 coincidences in the former category and 3 in the latter, a 
satisfactory but not particularly significant agreement. 

Two galaxies, NGC 300 and NGC 55, which are within the catalogue area, 
had previously been listed as probably observed (Mills 1955) but are not in 
Table 4. The former is listed in the footnotes to Table 2 as a possible coincidence 
with the source 00—317, that is, it is within 1™ in Right Ascension and 20’ in 
declination. However, the identification is not certain because of a complex 
brightness distribution in the region ; this difficulty was noted previously. The 
radio source associated with the latter galaxy is too faint for detection by the 
methods used in preparing the catalogue and is not included ; originally it was 
necessary to average several records to observe it. 

Finally, one should mention the two possible identifications found in the 
previously mentioned comparison of part of the catalogue area with the Palomar 
Sky Atlas. These are the sources 17—23 and 21—2/; they are discussed 
elsewhere (Mills 1960). 


IV. STATISTICS 

The question of the distribution of the radio sources has been discussed at 
some length in papers I and II. The present results do not alter substantially 
any of the earlier conclusions, so that a very brief discussion should be adequate 
at this stage. More useful statistical investigations may be carried out when 
the angular size data are complete. 

Previously we had divided the analysis into two parts, the first concerned 
with the two-dimensional distribution over the celestial sphere in an investigation 
of possible clustering effects, the second with the distribution in depth of the 
sources. These are not independent, but the division is convenient and will be 
adhered to here. 


(a) Clustering 

Earlier analysis had provided some slight, but not conclusive, indications 
of large-scale clustering of the radio sources. However, addition of the present 
data weakens the case and, on confining attention to regions of the order of 
10° to 30° across, no convincing evidence exists that the distribution is non- 
random. 

However, in the present catalogue a curious distribution arises, namely, 
the density of radio sources over the last few hours of the catalogue is very much 
less than over the first few. For example, the number of sources between 00 
and 02% igs 103 and between 22" and 002 only 56. The probability of this 
difference arising by chance in a random distribution is less than 0-001, Since 
the catalogue was prepared in a number of stages each commencing at 00" and 
working through to 24, it seems possible that the result is an artefact. However, 
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the catalogue of paper II was prepared in essentially the same manner and shows 
no such effect and, in addition, a careful re-examination of the original records 
supports the idea that the difference is real. The region of high source con. 
centration includes the south galactic pole, but no reason is known why this 
should provide a surplus of radio sources. 


(0} 
(0) 
NENG 
Se 
NS 
L NX 
100 \ A, 
XN 
Ye 
z 
< 
ra) 
< 
ira 
Ww 
EF 
oO 
Ho 
ae al 
WwW 
U 
© 
=) 
ie) 
oD 
SLOPE = 


1 i oe | 
10 20 40 80 160 320 640 


FLUX DENSITY (107 2©Wm~-2(c/s~') 


Fig. 1—Counts of Class II radio sources. 


In paper II the number of “ extended” sources was compared with the 
number of blends of radio sources expected on the basis of a random distribution, 
It was concluded that there was a very significant difference, indicating that 
either the sources were extended regions of comparatively low brightness, or 
close blends of physically associated sources of small size. The present results 
support this conclusion although, surprisingly, the proportion of extended 
sources listed is much less than before. For example, between Right Ascensions 
of 092 to 15» and 21" to 154 (i.e. at high galactic latitudes) the numbers of extended 
sources with flux density greater than 4010-26 W m-? (c/s)-1 is 9, and with 
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flux density >20x10-?* W m-?(e/s)-! is 19: the corresponding numbers 
expected on the basis of the previous catalogue are 16 and 29 respectively. 
Similarly, the number of sources with flux densities + 2010-26 W m-2 (c/s)-1 
listed as “ perhaps extended” is now 7 compared with an expectation of 14: 
however, the availability of data from the angular size interferometer has bern 
to eliminate some doubtful examples. Since the chance blends above the same 
flux density limits may, from the data of paper I, be expected to be 2 and 11 
respectively, the previous conclusion, that these are inadequate to explain the 
observations, is upheld. 


(b) Source Counts 

One of the most interesting questions arising from a study of the statistics 
of radio sources is their distribution in depth, which is investigated by means of 
counts of sources with fluxes above defined levels. Since there is a class of 
galactic source fairly closely confined to the plane it is usual to divide the analysis 
into two parts, dealing with the low latitude and high latitude sources separately. 
We will follow this procedure here, but, since the information about the galactic 
sources adds little or nothing to that presented in paper II and by Mills (1959), 
we will discuss only the high latitude or Class II sources (sources for which 
| b |>124°). 

In Figure 1 is shown the count resulting from a combination of the present 
catalogue and that in paper IJ. The actual numbers are listed in Table 5. 


TABLE 5 
THE NUMBERS OF CLASS II RADIO SOURCES ABOVE DEFINED FLUX DENSITIES 


Flux density (10-2 Wm-? (c/s)—*) bi a 10 20 40 80 160 320 640 


Number of sources, NN... Be ~ 2 L658 12205332 94 24 a 6 4 


The total area used in making these counts is nearly 5 steradians (2-72 
steradians in paper IT and 2-07 steradians in the present catalogue) and altogether 
there are 1658 sources with flux density >7x10-? W m-?(c/s)-? (the total 
number of Class II sources in the areas is 1700). Previously, the counts had been 
restricted at high intensities because of the fortuitous lack of such sources in the 
area of the catalogue; now, however, it may be extended much further, the 
range of flux density being nearly 100: 1. The standard error in each point in 
Figure 1 is indicated by vertical wings ; it is equal to \/N, where N is the number 
of sources counted. 

It has been shown in paper I that the counts need correction at low flux 
densities because of the effects of random noise and finite aerial resolution. 
Another effect was discovered when compiling the present catalogue, namely, 
a systematic tendency to over-estimate the flux density of a source when the 
signal-to-noise ratio is low. This appears quite distinct from the effects of noise 
previously considered ; it was detected by comparison of flux densities measured 
on “ scanning” and ‘“non-scanning”’ records. The former have a noise level 
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approximately double that of the latter and it was found that flux densities 
around 1010-26 W m-? (¢/s)-! were, on the average, estimated higher on the 
scanning records by a factor of about 30 per cent. For sources of double the 
flux density the effect was still noticeable but much reduced. These differences 
appear to arise because a source position is always selected for which the noise 
fluctuations combine with the deflection due to the source to produce the greatest 
apparent flux density. This position is, in general, displaced from the true 
position, which, if known and adopted, would result in a lower estimate for the 
flux density. The effect of this type of error on the source counts is not large ; 
this is fortunate, since it is difficult to assess an accurate correction without 
re-examination of all the original records. The most likely correction for the 
combined counts is a reduction of 30 per cent. for sources having fluxes near 
10 x10-26 W m-? (c/s)-!, and 10 per cent. for sources having fluxes near 
20 x 10-26 W m-? (c/s)-!.. In Figure 1 the crosses mark the corrected points, 
using both the data of paper I and the above corrections. 


A straight line of slope —1-5 is shown, drawn to fit the points as closely as 
possible; the fit is very good. Since this is the slope to be expected with a 
uniform distribution of source in a static Euclidean universe, the source counts 
indicate no divergence from uniformity and no obvious cosmological effects. 
It has been suggested by Mills (1960) that, among those sources which are 
associated with external galaxies, a source count having a mean slope of about 
—1-3 might be expected, if account is taken of red-shift effects. A line of this 
slope is shown dotted ; although the agreement is not particularly good, it is not 
inconsistent in view of the uncertain corrections to be applied at low flux densities 
and the statistical uncertainties at higher fluxes. It is known, however, that 
not all the high latitude radio sources in the catalogues could belong to this 
class (unless there is physical clustering of such sources) because of the inclusion 
of apparently ‘‘ extended ”’ sources in the catalogue, as discussed earlier. 


Finally, one should mention a slight difference in the counts if the two 
catalogues are analysed separately, the numbers in the present catalogue being 
systematically less at a given flux density. Part of the discrepancy at low flux 
densities may be due to the availability of more non-scanning records in the 
area of the present catalogue and consequently less systematic error in the weaker 
fluxes, whence, for a given number of sources, the flux density will be lower. 
The remainder may be the result of a small calibration difference ; this is likely, 
since the first catalogue includes zenith angles between 44° and 14°, the second 
includes zenith angles between 16° and 0°. A difference of 10 per cent. in the 
average calibrations in these areas would be adequate to explain the effect : 
this is about the accuracy of the calibrations, as discussed by Little (1958). 
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THE RADIO BRIGHTNESS OF THE QUIET SUN AT 21 CM 
WAVELENGTH NEAR SUNSPOT MAXIMUM 


By N. BR. LABRUM* 
[Manuscript received June 27, 1960] 


Summary 


An investigation has been made of the radio emission from the quiet Sun at 21-2 cm 
wavelength in 1958 (near sunspot maximum). Two different methods have been used, 
both involving observations with very high angular resolution, to distinguish between 
the quiet-Sun component and the radiation from localized active regions. In one 
method, the Sun was scanned with a narrow pencil-beam ; in the other, a fan-shaped 
aerial beam was used to give one-dimensional strip scans. In both cases it was necessary, 
when analysing the data, to take into account the residual effects of the very intense 
radiation from the localized sources. The two independent measurements gave results 
which agree within the limits of error. The apparent disk temperature was found 
to be approximately 140 000 °K, or twice the value for the same wavelength at sunspot 
minimum. 

The fan-beam observations also provide some evidence on the distribution of 
quiet-Sun brightness with respect to heliographic latitude. There is limb darkening 
at the poles, and the distribution does not appear to have changed in shape between 
sunspot minimum (1953) and the time of the present series of observations. 


I, INTRODUCTION 

The Sun’s radio emission at wavelengths of the order of 20 cm is made up 
of two main components. One of these varies slowly in intensity from day to day 
and originates in localized active regions in the solar atmosphere. These emitting 
regions may conveniently be termed “radio plages”’, by analogy with the 
optical plage areas with which they are closely associated. The second component 
consists of thermal radiation from the undisturbed or ‘“ quiet ” Sun, and remains 
at a constant level for periods of months or years. Radio measurements of the 
quiet Sun are of considerable importance, since they can provide information 
on temperatures and electron densities in layers of the solar atmosphere which 
are difficult to observe optically. 

The existence of a quiet-Sun component is inferred from the fact that, at a 
particular wavelength, a constant base level can be recognized in the varying 
intensity of the solar radiation. For an experimental investigation of the quiet 
Sun, it is necessary to find a method of separating the fixed component from the 
slowly varying plage radiation. This may be done by statistical analysis of a 
large number of daily observations of the total power received from the whole 
Sun at a chosen wavelength. An alternative is to use an aerial system with 
sufficient angular resolution to distinguish the active regions on the solar disk 
from each other and from the quiet-Sun background. Pawsey and Yabsley 
(1949) used a statistical method. A series of daily values of the apparent disk 
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temperature (at 50cm wavelength) was plotted against the corresponding 
projected sunspot areas. The distribution of the points on this diagram showed 
that the two quantities were closely correlated, and that the radiation consisted 
of a constant base level (quiet Sun) and a component proportional to sunspot area. 
The quiet-Sun temperature was estimated by drawing the line of best fit and 
extrapolating back to zero sunspot area. 

Many quiet-Sun determinations of this kind have since been made. Taken 
over a period of years, they show that the temperature at a given wavelength 
appears to increase towards the time of sunspot maximum (Christiansen and 
Hindman 1951). The technique is, however, open to objections because of the 
underlying assumption of a linear relationship between sunspot area and radio 
emission. There is evidence that individual radio plages generally have a 
considerably longer life than the corresponding sunspot groups. The method 
therefore tends, especially at sunspot maximum, to ascribe to the quiet Sun some 
radiation from active regions which no longer contain visible sunspots. 


Piddington and Davies (1953) and Dodson (1957) used alternative methods 
of analysis in attempts to overcome this difficulty. In each case it was found 
that, with the particular sets of data used, there was no clear indication of any 
change in quiet-Sun temperature over the sunspot cycle. On the other hand, 
Allen (1957) made a somewhat more elaborate analysis of the available observa- 
tions and concluded that a real variation does occur. His results indicate that 
the greatest change (about 2:1 during the cycle) takes place at wavelengths 
near 25cm. Christiansen, Warburton, and Davies (1957) independently con- 
firmed Allen’s conclusions. In view of these conflicting results, it is obviously 
desirable to make quiet-Sun observations by high resolution methods, and so 
avoid dependence on assumptions of correlation between radio emission and 
visible solar phenomena. 


II. High RESOLUTION TECHNIQUES 

This paper describes two separate determinations of the quiet-Sun 
temperature at a wavelength of 21-2 cm (frequency =1423 Me/s) in 1958, shortly 
after sunspot maximum. Both use the high resolution Christiansen crossed 
grating interferometer at Fleurs, N.S.W. (see Christiansen, Mathewson, and 
Pawsey 1957). 

This aerial system consists of two long arrays, each made up of 32 steerable 
paraboloid elements, 6 m in diameter and spaced 12-2 m apart on a 380 m base 
line. The base lines are aligned north-south and east-west respectively. The 
arrays can be combined in the same way as in the Mills Cross (Mills et al. 1958), 
to give a multiple response pattern of pencil beams about 3’ in half-power width 
and 1° apart. The separation between adjacent beams is sufficient to ensure that 
no more than one of them can lie on the Sun at any one time ; thus the multiple- 
beam response does not introduce any ambiguities in solar observations. <Alter- 
natively, either array can be used alone, to give a multiple fan-beam pattern 
with high resolution in only one direction. The system is thus suitable for either 
pencil-beam scanning over the solar disk (with the complete cross), or one- 
dimensional strip-scanning (with a single array). 

e 
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At sunspot maximum, a large part of the quiet Sun is at any time obscured 
by intense radio plages. Because of this, it has not been possible to observe 
details of the brightness distribution of the quiet-Sun component by either of 
these scanning methods. Both, however, were used to determine the quiet-San 
temperature ; the two independent results are in satisfactory agreement with 
each other. The fan-beam observations have also provided some information 
on the distribution of quiet-Sun radiation in heliographic latitude. 


OCT. 5.1958 


Fig. 1—A typical contour map of the 21 em solar brightness distribution, 
based on pencil-beam observations. Contour interval= 100 000 degk. 


III. PENcIL-BEAM MEASUREMENTS 

The solar data obtained with the crossed interferometer are summarized in 
daily contour maps of the radio brightness of the Sun; Figure 1 shows a typical 
example. It might be supposed that quiet-Sun temperatures could be read 
directly from each of these maps at points which on that day happened to be 
clear of plage areas, and that in this way a complete map of the quiet-Sun 
brightness distribution could be built up. The situation is, however, complicated 
by the presence of spurious responses to the radiation from plage areas. Like 
all aerials of the Mills Cross type, the crossed interferometer has rather large side 
lobes. The latter may be either positive or negative, and occur when a source 
is in a main lobe of the fan-beam pattern of one of the arrays. Their positions 
and magnitudes cannot be accurately predicted, as they are largely due to small 
accidental maladjustments in the aerial system. 
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Examination of the records showed that, on most days in 1958, the largest 
side lobes were about 4 per cent. of the main response (in terms of received power). 
From this, the effect of side lobes on the accuracy of quiet-Sun measurements 
was estimated for each of the records made between May and November 1958 
It was found that in most cases the range of uncertainty was at least 
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Fig. 2.—Estimation of quiet-Sun brightness temperature (at 
centre of visible disk) from selected pencil-beam records. The 
horizontal strips show the limits of error for individual records ; 
the vertical shaded band indicates the range of temperatures 
which is consistent with all the observations in the group. 


+100 000 degK, which is of the same order as the expected quiet-Sun temper- 
atures. Under these conditions, it was impossible to derive any significant 
information on the distribution of quiet-Sun brightness ; however, in a number 
of particularly favourable occasions there were points near the centre of the solar 
disk at which the quiet-Sun temperature could be determined within 
+50 000 degK. These selected observations were used in making an estimate 
of the central brightness temperature of the quiet Sun. 
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The temperatures were read from the original records, since the contour 
interval on the maps was inconveniently large for this purpose. The results are 
summarized in Figure 2. The mean value of the quiet-Sun temperature at the 
centre of the disk is 102 000 °K, and all the individual measurements are 
consistent with any value between 90 000 and 110 000 °K. 


IV. FAN-BEAM MEASUREMENTS 

For the one-dimensional strip scans, one of the interferometer arrays was 
used alone, giving a multiple fan-beam response. The method is a modification 
of that used by Christiansen and Warburton (1953a, 19536, 1955), who obtained 
sets of daily strip scans across the Sun in chosen directions and found the lower 
envelope of each set. These envelopes were the corresponding quiet-Sun profiles ; 
the area under any one of them was a measure of the flux density of the quiet-Sun 
component. 

The plage radiation is so intense at sunspot maximum that it becomes 
impossible to isolate the quiet-Sun component by this simple lower-envelope 
technique. The active regions are, however, mainly confined to narrow bands of 
heliographic latitude; it will be shown that, because of this feature of the 
distribution, a quiet-Sun profile can be derived even at sunspot maximum, 
provided that the scans are taken in a direction parallel to the solar axis. ) 


The use of strip scans is a less direct approach to the problem than is the 
pencil-beam method described above. However, it avoids appreciable errors due 
to side-lobe responses ; in the fan-beam pattern these are always very small 
(less than 0-2 per cent. in received power). 


(a) Observations 

The aerial used was the north-south array of the crossed interferometer. 
In directions normal to the array, where the angular resolution is greatest, the 
width of the fan beams (to half-power points) is 2’ are ; in the oblique directions 
in which the Sun was observed during this work, the beamwidth is about 3’. 
The directions of maximum response, when projected upon the celestial sphere, 
appear as small circles centred on the line of aerials. In Figure 3, a typical 
aerial beam is plotted in coordinates of declination and hour angle, together 
with the apparent diurnal path of the Sun across the sky. As each beam crosses 
the Sun, the receiver output records a one-dimensional strip scan. The direction 
of this scan, relative to the celestial meridian, depends on the solar hour angle 
and the declination, whilst the angle between this meridian and the Sun’s axis 
varies annually over a range of +26°. The scanning angle 0 measured from the 
central meridian of the Sun (see Fig. 3) therefore depends on the time of day and 
year; by a proper choice of observing time, scans parallel to the solar axis 
(8=0°) can be obtained. 

The records used for this work were made between September and 
November 1958. They consist of 20 strip scans taken on different days; the 
scanning directions are all within +-5° of the Sun’s central meridian. 
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(b) Method of Reduction 
(i) Positions of Scans.—The most precise method of locating a sean in 
relation to the Sun is to calculate the time at which the centre of the appropriate 
aerial beam crosses the centre of the solar disk. It was found to be sufficiently 
accurate, however, simply to determine this time by measurement on the record 
to the mid point of the scan. In all cases where comparisons were made, the 
times so estimated were within 5 sec of the calculated values. 

The rate at which the aerial beam moved across the Sun was calculated for 
each observation. The scale of time in the original record was then converted 
into one giving the perpendicular distance from the centre of the scanning strip 
to that of the solar disk. To compensate for the small annual variation in the 
angular diameter of the Sun, all positions were expressed in terms of the apparent 
diameter of the visible disk at the time of observation. 


DECLINATION 


APPARENT 
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Fig. 3.—Strip-scanning across the Sun with a fan beam. The shaded 
circles are successive apparent positions of the Sun, at approximately 
4-min intervals. 6 is the scanning angle, as defined in the text. 


(ii) Intensity Calibration.—The overall sensitivity of the equipment inevit- 
ably varied somewhat from day to day. A complete daily calibration was 
impracticable, owing to the difficulty in measuring losses in the complicated 
aerial and feeder system. Instead, a relative measure of the total energy received 
from the whole Sun was found for each day by integrating the strip-scan record ; 
this was then compared with an absolute value of the total flux density obtained 
from an independent direct measurement. In this way a multiplying factor was 
found, to bring each of the scans to a common intensity scale. 

The Sydney 1420 Mc/s radiometer was unfortunately not operating con- 
- tinuously over the required period ; however, daily flux records at 1500 Mc/s 
were available from the Heinrich Hertz Institute, Berlin, Germany. These were 
adjusted to 1423 Mc/s on the assumption that the flux density is proportional to 
frequency. An examination of published data at frequencies from 1000 to 
2000 Me/s indicates that, for the small frequency change involved in this case, 
this procedure is unlikely to lead to errors of more than 1 per cent. in the derived 


values. 
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(iii) The Quiet-Sun Profile—A typical scan with 6=0° is shown in Figure e! 
while 20 such scans have been superimposed in Figure 5. Thee Ring TAnis aps 
clearly the characteristics of the distribution of plage Bear with heliographic 
latitude. At sunspot maximum, practically all the major radio plages are centred 
in the middle latitudes ; strong sources are unusual within +10° of the equator 
and very rare in latitudes higher than +40°. 
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Fig. 4.—A typical strip scan, with scanning direction parallel to 
the central meridian (9=0°). AR is the radius of the photosphere. 


With the scanning direction parallel to the central meridian, therefore, the 
quiet-Sun profile is seriously obscured only in positions corresponding to the 
middle latitudes. The complete profile can be estimated by determining the 
equatorial and polar sections and drawing a smooth curve through the intervening 
unobservable region. The lower envelope is clearly defined in the polar regions ; 
at the centre of the disk, however, there is a large scatter between individual 
scans, and a more elaborate method is necessary to determine the quiet-Sun value. 
Sources centred near the solar equator occur so infrequently that their contribu- 
tion is likely to be unimportant; any excess over the quiet-Sun flux at the 
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central minima in the scan must be due mainly to radiation from the outer parts 
of the intense active regions centred in middle latitudes. The influence of these 
is accentuated by the smoothing effect of the finite aerial beam. The problem is, 
therefore, to subtract the component due to these strong sources from the observed 
‘central’ values. 
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Fig. 5.—A set of 20 scans at 0=0° (Sept.-Nov. 1958). 


This has been done by correlating the observed minimum value, p,, of the 
flux near the centre of each record with a parameter p,, which is used as a measure 
of the combined effect of the middle latitude sources. The choice of p, 18 
somewhat arbitrary. An obvious possibility is to use the average of the two 
maximum flux values on the record. The positions of these maxima, however, 
vary somewhat from day to day, and it appeared preferable to take Py as the 
average of the two values at equal fixed distances on either side of the minimum 
position. A suitable spacing was found to be one-sixth of the solar radius, which 
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was in all cases rather less than the distance between the positions of maximum 
and minimum flux density, as indicated in Figure 4. 

In Figure 6, p, is plotted against p,. The coefficient of correlation between 
the two quantities is 0-80, and the points in the diagram have a clearly defined 
lower envelope—the straight line AB. The slope of this line is interpreted as 
representing the contribution of the middle latitude sources to p, ; the upward 
displacement from it of individual points corresponds to the occasional presence 
of weak sources close to the solar equator. 

It is now necessary to find what part of p, corresponds to quiet-Sun radiation. 
For this purpose, it is tentatively assumed that the central part of the quiet-Sun 
profile is flat, so that p,—p,, in the absence of any radiation from plage areas. 
On this assumption, the quiet-Sun value for p, is given by the intersection of AB 
with the line p.=p, (Fig. 6). 
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Fig. 6.—Graphical determination of central ordinate on quiet-Sun 
profile. The units for » and », correspond to the vertical scales in 
Figures 4 and 5. 


This value of p, was combined with the polar sections of the original lower 
envelope, to give an estimate of the north-south profile of the quiet Sun (Fig. 7). 
The middle latitude part of the curve is not fixed by this procedure, but is merely 
interpolated by hand. It appears probable, however, that the true profile lies 
within the shaded area in the diagram. For any such profile, p, does not differ 
from p, by more than 3 per cent., so that the initial assumption that p,—p, 
is not invalidated. 

A possible source of error in this method is the influence of sources centred 
near the equator on the measured value of p,. Pencil-beam records for the dates 
in question show that any such sources were always too small for an effect of this 
kind to be appreciable. 

(iv) The Quiet-Sun Flux Density—The flux density of radiation from the 
quiet Sun can now be calculated, since it is proportional to the area under the 
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profile in Figure 7. The constant of proportionality is determined by comparing 
the area under any one of the individual daily scans (on the same scale) with the 
corresponding observed flux density. The results range between 5-7 and 
6-6 x10-*! Wm-? (c/s)-1 for the profiles within the shaded area of Figure 7. 
The corresponding apparent disk temperatures (for a uniform disk equal in size 
to the photosphere) are 130 000 and 150 000 °K. 
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Fig. 7.—The quiet-Sun profile for 0=0°. For comparison, the 
corresponding sunspot-minimum profile (Christiansen and 
Warburton 1955) is drawn on the same scale. 


V. DISCUSSION 
(a) Apparent Disk Temperature 
The apparent disk temperature of the quiet Sun in 1958 has now been 
determined, and an independent estimate has also been made of the brightness 
temperature at the centre of the disk. In order to compare these two results, it 
will be assumed that the brightness distribution in 1958 was the same as that 
found by Christiansen and Warburton (1955) at sunspot minimum. It will be 
shown later that this is a reasonable assumption. The pencil-beam value of 
90 000-110 000 °K for the central temperature then gives a disk temperature of 
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130 000-160 000 °K. This agrees with the value derived from the fan-beam 
observations, within the limits of error of the two experiments. For 
the remainder of this discussion, the apparent disk temperature for the 
quiet Sun in 1958 will be taken as 140 000 °K. The corresponding flux density is 
622 107! Wm (ef): 

This result shows that the quiet-Sun temperature increased by a factor of 2 
between sunspot minimum and sunspot maximum ; Christiansen and Warburton 
(1955) give a value of 70000 °K for the apparent disk temperature in 1953. 
There is thus good agreement with Allen’s (1957) calculations, which indicated 
a variation of 2-0: 1 and were based on observations at 25 cm wavelength in 
the period 1947-1953. The value of the quiet-Sun flux at sunspot maximum 
also agrees well with an earlier estimate. Christiansen, Warburton, and Davies 
(1957) found the flux at 21 em by interpolating between published 10 and 25 em 
results for 1947, and obtained a value of 6-4 x 10-7! Wm > (c/s). 


(b) Brightness Distribution 

For sunspot minimum, Christiansen and Warburton (1955) constructed a 
complete map of the quiet-Sun brightness distribution from a series of profiles 
at various scanning angles. This was not practicable in 1958, owing to the high 
level of solar activity. Some inferences can, however, be drawn from a com- 
parison of the north-south profiles for 1953 and 1958. The 1953 curve given by 
Christiansen and Warburton is shown in Figure 7. The two profiles are similar 
in shape ; the ratio of corresponding ordinates in both the equatorial and polar 
sections is approximately 2:1. It therefore seems likely that the brightness 
distribution is much the same at sunspot maximum and minimum. In particular, 
the pronounced limb darkening at the poles, which was detected from the earlier 
observations, was also present in 1958. The same conclusion follows from recent 
eclipse observations at 21 cm wavelength (Krishnan and Labrum, in preparation). 


There is no detectable variation over the sunspot cycle in the polar diameter 
of the radio Sun at this wavelength. The 1953 profile is slightly the broader of 
the two, but this difference can be explained by the lower angular resolution of 
the earlier aerial system. 


(c) Interpretation of the Results 

The primary purpose of this investigation was to determine whether the 
observed radio emission from the quiet Sun changes during the sunspot cycle. 
This has been done by making measurements at a single suitable wavelength. 
In order to make any detailed deductions regarding the temperature and density 
in the solar atmosphere, similar data for a range of wavelengths would be 
necessary. It appears, however, that the present empirical results are consistent 
with other evidence. In particular, van de Hulst (1949) estimated the coronal 
electron density from optical measurements and found that the density at any 
height increases towards sunspot maximum. From this result and assuming 
that the kinetic temperature of the corona does not vary, he predicted that the 
brightness temperature of the Sun in the decimetre-wavelength region would 
show changes of the order of 2:1 during the sunspot cycle. 
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VI. CONCLUSION 
The flux density of radiation from the quiet Sun at 21-2 cm in 1958 was 
approximately 6-2 10-2! Wm? (c/s)-1;_ this corresponds to an apparent disk 
temperature of 140 000 °K. The brightness of the quiet Sun at this wavelength 


increased by a factor of 2:1 from 1953 (sunspot minimum) to 1958 (slightly 
after sunspot maximum). 


The data indicate that the distribution of quiet-Sun brightness with helio- 
graphic latitude does not change appreciably during the sunspot cycle. In 
1958, as in 1953, there was substantial limb darkening near the poles of the Sun. 
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DIRECTION-FINDING ON DIFFUSE SOURCES OF 
ELECTROMAGNETIC RADIATION 


By D. G. CARTWRIGHT* 
[Manuscript received June 27, 1960] 


Summary 


It is shown that, for sources of large angular size, the response of an Adcock type 
direction-finder is independent of the extent of the source in altitude. On the other 
hand, the response of a rotating loop does depend on altitude. By combining the 
characteristics of both types of direction-finder, the position and size of an extended 
source can be found, provided that a brightness profile can be assumed. 


I. INTRODUCTION 

Observations of the time variations of very low frequency emissions from 
the Earth’s upper atmosphere (Ellis 1959) have shown the need for a means of 
estimating the size of the apparent sources of the radiation. Investigations so 
far have been made at a frequency of about 5 ke/s, a long wavelength inconvenient 
for use with conventional high resolution antennae. However, information 
about the angular distribution of the recorded radiation can be obtained with 
relatively simple antennae such as loop and Adcock direction finders. For 
instance, using rotating loops, Ellis and Cartwright (1959) have found that the 
radiation sources have an apparent azimuthal size of between 60° and 90°. 

In a recent paper (Wait 1959) it has been demonstrated that the 
characteristics of a downcoming radio wave may conveniently be measured by 
making use of rotating loop and Adcock direction-finders. Expressions were 
derived for finding the angle of arrival, the azimuth, and the polarization, in 
the case of a single plane wave-train incident on the observation point. 

In the present paper the response of a loop and an Adcock system will be 
analysed for sources of large angular size. It will be shown that the direction- 
finding characteristics of an Adcock system are independent of the size of the 
source in elevation. 

If the source is extensive in altitude as well as in azimuth, then the directional 
properties of a rotating loop are modified by a term containing the altitude. 


II. COORDINATES 
The aerial systems will be considered relative to a spherical coordinate 
system, with the Earth a conducting surface in the vy plane and the antenna in 
the w-z plane. An incident wave is specified in terms of its spherical coordinates : 


angle of elevation above the w-y plane, 9, and its azimuth measured from the 
z-axis, 9. This is shown in Figure 1. 


* Upper Atmosphere Section, C.S.I.R.0.. Camden, N.S.W. 
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I. Apcock AERIAL 
A pair of suitably connected antennae rotating about a vertical axis midway 
between them is equivalent to a conventional 4-wire Adcock and goniometer 
arrangement. The response of this system to a small source will be considered 
first. Ideally, it responds only to the vertical component of # of the incident 
wave. If the antennae are on the w-axis and close to a good conductor (owe, 
where o is the conductivity and < the permittivity of the surface beneath the 


az 


ADCOCK 


ELEMENTS INCIDENT WAVE 


Fig. 1.—Coordinate system, showing position of loop and Adcock 
antennae. 


loop in appropriate units, and w is the angular frequency of the wave) the power 
transferred to the receiver is proportional to cos* 9 cos? 0, and the point source 


response can be written 
p4=cos' ¢ cos? 0. 


If the extended source is a random emitter with uniform brightness then the 
power at the receiver can be written 


P,=| [padds. 


Assuming the source has a uniform brightness over an azimuthal width of 
2@ and extends in elevation from zero to @, then the maximum power recorded 


as the aerial rotates is ert 
Pras. =2 | | p,dbdo. 
oJ0 


This gives : 
Pras. =4(O+9 sin 20), 


P 
a= | cos* edo 
0 


i/3 : TeV 
= (50 +80 20+ asin 0), 
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and at the minimum 
4m oD) 
ange { | p dvds, 
in—OJ 0 
Pin. = «(© —4 sin 20). 


Defining the modulation in power for the Adcock, M ,, resulting from the antenna 


rotation as 
Big ee —Prin. 


en dee (1) 
Prax. +P min. 
gives 
sin 2) 
M=5 STEMI a Pe (2) 


Since the spacing of the dipoles at very low frequency is very much smaller than 
a wavelength, the minimum power corresponds to the azimuthal centre of the 
source in the usual way. 

Thus an Adcock gives information about the position and size in azimuth. 
The relationship is independent of the brightness profile in ®. 


IV. Loop AERIAL 
The case of a small rotating loop will now be considered. The magnetic 
components of the incident wave are H, in the plane of incidence and H,, normal 
to this plane. For a loop in free space the Cartesian components due to H, are 


H,=H, sin ¢ cos 8, 
H,=H4H, sin 9 sin 9, 
H,=H, cos 9, 

and those due to H, are 


For random polarization, H,—H, and the powers due to each component can 
be added. The loop absorbs power only from H,, so the small source power 
response for a small loop in free space is proportional to 


cos? 0+ sin? sin? 0. 


This is unaltered if the loop is close (in terms of the wavelength) to a good 
conductor in the #-y plane, so we can write 


p,=Cos? 0+8in? ¢ sin? 0. 


For an extended source of uniform brightness, the total power at the receiver is 


P= | [rsddde. 
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If, as before, the source is of width 2 in azimuth and extending from 0 to ® 
in elevation, the maximum power recorded is 


Prax. =(1+8)0 +3(1 —8) sin 20, 


where 


and as before 

Prin. =(1+8)0 —4(1 —8) sin 20. 
The modulation in power for the loop then becomes 
1—B\sin 20 
M, feo) 510 (3) 


M, is indistinguishable from M, if 8<1, ie. ® is small. 
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Fig. 2._The azimuthal size 2@ of a diffuse source, shown as a function of the modulation in 
power, M, produced by a rotating loop direction-finder, for different source heights ®. The case 
“@=0 corresponds to an Adcock direction-finder. 


V. COMPARISON OF Loop AND ADCOCK RESULTS 
The modulation in power of an Adcock direction-finder is independent of 
the brightness profile in elevation, but, if the brightness distribution in azimuth 
is symmetrical, the minimum recorded power corresponds to the direction of 
the centre of the source in the usual way. 
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The power modulation of a rotating loop, on the other hand, is aftected by 
the distribution of brightness in elevation if the source extends more than about 
5° above the horizon. 

By comparing the records from an Adcock and a rotating loop the height of 
the source can easily be found. 

For instance, 


M, 1+8 
M, 1-8’ 
or 
M,—M, B—4(0—}F sin 20) 
M, +M, 5, 
ae 1-0 
o-9 iT. 
0-8 re 
O-7 ay 


ro) 
o 
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Fig. 3.—The relation between the modulation in power, 
M=(Pmax.—Prmin.)/(Pmax. +P min.) 
and the ratio Pmin/Pmax.» OT Vmin./ Vmax.= V(Pmin./Pmax.)- 


VI. Discussion 
Figure 2 shows the expressions (2) and (3) with 20 plotted as a function of 
the modulation in power M. The most convenient quantity to record is the 
amplitude or voltage produced by the direction-finder. From (1) we have 
M=(Prax, —Prin.)/(Pmax. +Prmin.)- 
If we denote Vinin,/ Vmax. by v, then v=+/(Pmin,/Pmax.), and so 
O=4/{ (LB) | (LS ae ee e e (4) 


The relationship between or Pyin,/Pmax, and M is shown in Figure 3. Using either 
(4) or Figure 3 we can obtain the source size in azimuth as a function of v for any 
given ®. This is shown in Figure 4. 
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This method is best suited to sources for which ® ig greater than about 30°. 
From Figure 2 it can be seen that the curves for ® less than 30° are very closely 
spaced in comparison with curves for ® greater than 30° and, in fact, the accuracy 
of measurement is not high enough to resolve curves for which © is less than 20°. 
An exception to this is when the source is of small size in azimuth, since the 
curves of Figure 4 are quite widely spaced below 20 =40°, say. 
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Fig. 4.—The azimuthal size 2© of a diffuse source shown as a function of the 

depth of the amplitude null, Vyyin./Vmax,, produced by a rotating loop direction- 

finder, for different source heights ®. The case D=0 corresponds to an Adcock 
direction-finder. 


VII. ACKNOWLEDGMENTS 
The author acknowledges many suggestions made during discussions with 
Dr. G. R. A. Ellis of this organization. 


VIII. REFERENCES 
Exus, G. R. A. (1959).—Low frequency electromagnetic radiation associated with geomagnetic 


disturbances. Planet. Space Sci. 1: 253-8. ' 
Exus, G. R. A., and Cartwricut, D. G. (1959).—Directional observations of radio noise from the 


ter atmosphere. Nature 184: 1307. + i 
o's ins sche vewatetine radio waves measurement of characteristics. Electronic & 


Radio Engr. 36: 106-7. 


THE GENERAL THEORY OF THE MOTIONS OF IONS AND 
ELECTRONS IN GASES 


By i G2 oH. Hoxtey’ 
[Manuscript received June 20, 1960] 


Summary 


In this paper the more important general formulae for the drift velocities and 
diffusion coefficients of ions and electrons in gases are derived by the application of 
dynamical principles. These formulae agree with those already established by a proper 
application of the method of free paths. Formulae for the distribution of speeds of 
agitation are also derived. 


I. INTRODUCTION 

There is an extensive literature on the theory of the motions of ions and 
electrons in gases but the derivations of the formulae for drift velocities, 
coefficients of diffusion of ions and electrons, and the conductivities of weakly 
ionized gases in direct or alternating electric fields accompanied or unaccompanied 
by a magnetic field, that find most frequent practical use, are not readily 
accessible. 

In what follows the chief formulae are established in a general form by the 
application of dynamical principles. Although the treatment is purposely 
elementary it does not lack rigour. 


The aim is to provide a compact and uniform summary of the theory of the 
subject in a useful form (general references: Allis 1956 ; Margenau 1958). 


IJ. ELECTRONIC MOTION IN GASES 
(a) General 

When electrons move freely among the molecules of a gas in the absence of 
an electric field they interchange energy and momentum with molecules in 
collisions and their steady state of random agitational motion is one in which at 
any instant the directions of their velocities c are distributed isotropically and 
their speeds ¢ are distributed according to Maxwell’s formula which states that 
the proportion of a group of n electrons whose speeds exceed ¢ but do not exceed 
c-+de is 


dn, 4 
ce Oxp (= oF] a2 ode. sui ee (1) 


where « is the most probable speed. In addition a condition of equipartition of 


energy prevails in which the mean kinetic energy of agitation 4me? of an electron 
is equal to that, }MC2, of a molecule. 


* Commonwealth Scientific and Industrial Research Organization, Melbourne ; present 
address: Australian National University, Canberra. 
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This motion becomes modified in important respects when electrons move 
in a steady state of motion in a gas in the presence of a uniform and constant 
electric field E. It now comprises a steady drift velocity W of the centroid of 
the group with a superimposed random motion of agitation such that the speed W 
is much smaller (a few per cent.) than the mean speed of agitation ¢é. 

The speeds ¢ do not, in general, conform to Maxwell’s distribution formula 
and the mean kinetic energy of agitation 4mc? of an electron exceeds that, 4MC?, 
of a molecule by a factor k (Townsend’s energy factor) which is a function (specific 
to each gas) of the ratio H/N of the electric field strength to the number of 
molecules V in unit volume of the gas. When the measurements are referred 
to a standard temperature (15 °C) the parameter H/p is more commonly employed 
than E/N, p being the pressure of the gas. 


These and other aspects of electronic motion in gases are considered in greater 
detail in what follows. 


Fig. 1 


(b) Nature of the Distribution Function in a Steady State of Motion 

Consider a general isotropic distribution of the speeds ¢ represented by a 
function f(c) with the meaning that the proportion of electrons of a group n whose 
speeds exceed ¢ but do not exceed c+de is dn,/n=4rf(c)c*de. 

This function is represented in velocity space by a spherical distribution of 
points, such that the number of points contained within an element du dv dw 
of velocity space at a distance c(u,v,w) from the origin is n f(c)du dv dw. ; Consider 
the unsymmetrical distribution F(c,w) obtained from f(c) by increasing the uw 
component of every velocity ¢ by an amount V(c) as shown in Figure i “Lhe 
centre of symmetry has been displaced from O to O’ through a distance V, 
thus leaving unaltered the number of points in the element of space dz at a 
and P’ respectively, but the velocity ¢ increases from c=OP on the left to c=OP 
on the right. Let 0 be the angle between OP’ and the w-axis and suppose that 
_ V<c, then cos §=u/e and O’P’=c—V cos 0=c—Vule. 

The number of representative points in dz on the right is 


nBlo,n)az—nf(e— oY Jae 


whence 
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It remains to determine the value of V and the form of the function f(c) associated 
with steady electronic motion in a gas under the influence of a steady and uniform 
electric field E, it being assumed throughout that the number of electrons in 
unit volume is much smaller than the number of molecules in unit volume so 
that mutual interactions of electrons are unimportant. 


(c) Specification of an Encounter 

In Figure 2 (a) an electron or ion e with mass m approaches a molecule O 
with mass VM at a relative velocity ¢ and is deflected to leave the vicinity of O 
at a velocity g’. It is supposed that when e is sufficiently distant from a molecule 
its trajectory is rectilinear. 

In Figure 2, b is the perpendicular distance of O from the direction of § 
and § the angle of deflection (the angle between g and g’). 9(b,g) is a function 
of b and g. 


Fig. 2 


It will be assumed that there is a limiting distance o such that if b exceeds o 
the deflections 0(b,g) are zero or so small that even in aggregate they are 
unimportant. 

The vector diagram of velocities relating to an encounter is (following 
Maxwell) depicted in Figure 2 (b). c, and c, are the velocities respectively of 
the ion (or electron) and the molecule before an encounter and C1 and Co those 
after the encounter. 

G is the velocity of the centroid X of m and M and §, and g, are those of 
‘mand M relative to X before the encounter and 2) and So those after it. It 
follows that $=$,+8, and $’=$) +5. 

In a system in which the gas as a whole has no mass motion, and to which 
the velocities c refer, the actual deflection of ¢ is the angle 0, between c, and c, 
and in general this angle depends not only upon b and g, but when c, is given, 
upon c,, However, when m/M<1 and also c,>¢,, as is the case if e is an electron, 
then ¢~c, and 0(b,g) =6(b,c,)=0,. Also ¢; differs little from c, in most instances. 


An important quantity is cos 0, the average value of cos 6, taken over all 
encounters in which the ion or electron is travelling with speed c, before an 
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encounter. Since the velocities c; are distributed with axial symmetry about 
the direction of c,, in the case of an electron where 0,=6(b,c,), 


o 


cos 0,= Ear COS 6(6,e)bab| [rot=2/0t| cos (0,¢,)bdb. 
0 


The mean value of cos§, averaged over the whole distribution of speeds c¢, 


would be written cos 6,, that is to say it is the mean value of 9 in an encounter 
of any kind. 


(d) Free Paths 
Consider a large number p of free paths 2, ,. . ., v, all traversed at the 
Same speed ¢ by the ion or electron, but not necessarily consecutively. 
If J), is the mean free path, then 


n 
Plo. =La,. 
1 


Let g be the mean speed of the ion (with speed ¢) relative to the molecules which 
move at random and let N be the number of molecules in unit volume. 

The sum of the times spent in traversing the paths x, is t=pl,,/c, that is to 
say, p=ct/l,,. Butp=gt.Nxo?, consequently, if 1,=1/N ro, ly, =c/9.N mo? =(¢/g)ly. 
For electrons g=c; 1|,,=1/Nzxo?=l,. o is the limiting value of the impact 
parameter beyond which deflection of the ions are unimportant. 


Fig. 3 


(e) Calculation of V 
According to Section 2 (b), the velocity of an electron is the vector sum 
V-tc of a fixed velocity V and velocity of agitation c. 
Consider the vector diagram of velocities of a collision of such an electron 
with a molecule with momentum MC (Fig. 3). 
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In Figure 3, c.=V+c; G=(mc,+MC)/(m+W) is the velocity of the 
centroid XY of m and M, and r is the velocity of the electron (or ion) relative to X. 
These are the velocities before the encounter. After the encounter the corres- 
ponding velocities are c’, C’, G, and r’. 

If no change results in the internal energy of the molecule then r’=r. 

Consider first the case in which all directions of r’ are equally probable 
when C and c, are given, as would occur in encounters between rigid smooth 
spheres (Maxwell 1890). 

The mean velocity of electrons and molecules (c, and C as specified) after 
encounters is therefore G. Consider next, the mean residual velocity in all 
encounters of electrons for which c, is given but C is allowed to range over all 
possible directions. Draw XP parallel to C to meet c, at P. Then G=OP+PX. 
But OP={m/(M+m)\c, and XP={M/(M-+m)}C, consequently P is a fixed point 
and XP has a fixed length when C is constant. The possible end points of C 
lie on the surface of a sphere about O and those of G on a sphere about P. Thus 
the mean projections of both C and PX on the direction of c, are zero. The mean 
residual velocity after all such encounters is OP. Finally c may assume all 
directions with respect to V and therefore QP does likewise. The mean residual 
velocity is therefore OQ={m/(m-+M)}V for encounters in which the agitational 
speed has a fixed value c. 

Suppose that all directions of r’ are not equally probable but that r’ is 
distributed with axial symmetry about r. Since the scattering occurs with axial 
symmetry about the direction XQ, the mean values of the components of c’ and r’ 
normal to XQ are zero, so that the mean value of the vector c’ when the scattering 
about YX is not isotropic is a vector parallel to r. It may be written ar. 

But, from the triangle PXQ, «r=a{M/(M-+m)}(c,—C), so that when C 
is distributed over all possible directions in space the residual mean value of c’ 
is a vector parallel to c, and equal to {«M/(M-+m)}c,. But c,=V-+c, so that 
when c ranges over all possible directions the mean value of c’ reduces to 
{a M/(M-+m)}V, where «cos 0, where 0 is an angle of deflection of the velocity r 
relative to X in an encounter when the speed of the ion or electron relative to a 
molecule is g={(M-+m)/M}r. 

The mean residual momentum after encounters is therefore 
{(m+aM)/(M-+m)}mV, which is the same as the mean momentum that the ion 
possesses when it enters the speed group c. The mean momentum lost in an 


encounter is 
m+aoM mM 
mv(1— MU tm eee Me 


and the mean rate at which an ion or electron loses momentum in encounters is 


mM 
M-+m 


(1—a)'V, 
where J,=1/N7o?. 


In a steady state of motion the mean momentum mV is the sum of the mean 
momentum {(m-+«M)/(m+M)}mV at entry and that acquired from the field, 
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namely, Eel,/g. The mean momentum imparted by the field is therefore equal to 
the mean momentum lost in an encounter. Consequently, 


Eel) mM 
g m+uo*Ys 
or 
Vee JECT ee Ol ab) tae, Se ens (3) 
where 
lag ao “x=cos 0 


Lis an equivalent mean free path and g is the mean speed of the ion or electron 
with agitational speed c, relative to the molecules. For electrons, g~c, m/M<l, 
1=I,/(1—a). The equivalent cross section is here 


A(c)=1/NI=(1 ~0550)/t4=27 | {1 —cos 0(b)}bab. 
0 


When electrons and ions interact as point centres of repulsive force P=k/r’, 
then, if v>2, the integral 


a 


on | ” {1 —cos 0(b)}bdb 
0 
is convergent and defines an equivalent cross section and mean free path 
1=1/2nN i ” (1 —cos 0(b))bdb=1/NA(C), 
0 


where A(c)=27BA,(v)/c/@-D, in which B=(k/m)2/-) and A,(v) is a function 
of v only (Chapman and Cowling 1952, p. 171; Huxley 1957a, p. 125). For 
instance, when v=5; A(c)o1/c. 
If the scattering is isotropic, 
A=2n[’ bdb=no®= Ay. 


0 
In this event, «=—0, 


ya ala) _ Bet = 


mg M g \m' My 


In general, J=1,(M-+m)/M(1—«) may be written J=1,+S (Huxley 1957, 1960), 


where 
S=(n aM )l,jM(1l—a).  neccserscesseses (5) 


(f) Formula for Drift Velocity W in a Steady and Uniform Field E 

The mean speed of drift in the direction of E (parallel to Ox) of the group of 
electrons with speeds c is the mean value of the component w of the velocities c 
averaged over all directions in space. From equation (2) it follows that the 
mean value of w=c cos 0 is 


fee es oe 


W (c)=ef(e) cos 0 ag : Be do! 
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so that, from equation (3), 


The mean value of W(c) taken over all speeds ¢ is the drift speed W of the centroid 
of the whole group. Thus, 


He (” 4 ele dj 


—- = ae 
3M J 5 g de 


Ee( Azle? .| so gl vel “4 ; ‘ 

=s—i— 4 — ~|—])f.c?de 

3m | g uf 7 an ce aay f 

> He = a fle® 7 
3m . (6) icq) a6 2 a ee 6. 6 6 Sk Oe se... * & 0 oe Bee © Cle ( ) 


where 1=1,-+S ; (equation (5)); g is the mean speed of an ion or electron relative 
to the molecules. For electrons, m/M<1, g=c, equation (7) becomes 


_— 


the bar denoting an average with respect to c. 


This formula was previously derived directly by a correct application of the 
method of free paths (Huxley 1957a, 1960). 


Throughout, the atomic charge e is regarded as algebraically positive, that 
is to say, a negative value should be substituted in the case of electrons. 


(g) Drift Speed in an Alternating Electric Field 
Let the electric field be considered to be a rotating vector in the XOY plane, 
E=E, expipt. It is assumed that the frequency p/27 is sufficiently large that 
fluctuations in the mean energy 4mc? are unimportant. This assumption is also 
correct for the special case p=0. 
An expression for the velocity V is first obtained. The momentum mV 
fluctuates in a time-dependent field because its rates of loss by encounters and 
of gain from the field are not equal. Thus, from Section II (e), 


dv : M 
ms =H, exp (ipt) — (Wan —a) inv, 


or 
dVy G55 hee : 
ai sit core CXP (1b), 5 a sire em other aie «Gua wines ee (9) 
whence 
_ ee 5 
ear eXp (pt), . we bes etcetera (10) 


where v=g/l. 
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When the vector V is not directed along an axis of coordinates the appropriate 
form of equation (2) is 


F(6,u,0,w)=f(e) — — =. (11) 


So that, when H=H, exp (ipt)=H,+i#, and V=V, +iV,, the appropriate form 
of equation (11) is 


F(ejuy0) =f(e) — eR ot 


It follows, as in Section II (f), that, with H=E, exp (ipt), 


ee (oe anc df, 
3m}, (+ip) de® 


Ee (® Oly Ge 
=>=— as 2 
3m} le a ( =| 4rc*fde 


w= 


(h) Drift Velocity in the Presence of a Magnetic Field 


Let a magnetic field B act along the direction Oz and H along Ox and write 
@=—Be/m. Then equation (9) is to be replaced by 


dV _ Kye F : 
ape om °xP (ipt)+iwV, 


or 
se (9405) V = —= OX IT) eee dane ete. eee (13) 
Whence 
= ee co. dete’ Segoe (14) 
and 
Wa eect Pore: (15) 


Equations (12) and (15) can be derived directly by use of the method of free 
paths (Huxley 1957b). 


With electrons v=g/l becomes v=ce/l. 


(i) Magnetic Deflection of an Electron Stream in a Gas 


In equation (15) let p=0, then 


ge A ee We Ee _,4[(v+io)e | 
W=W, exp 0) =F) Sm “Ab aseste awe (16) 


«= —Be/m. 
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Choose the axes of coordinates so that +Oy is parallel to B and +0z to E, then 
it follows from equation (16) that 


W=W,+iw,=W, (cos 6+ sin 9). 


he kf we 
We =f “alte =) 


3 
WF O56 £ ) Joha: 


t ga Vs , 5a e ek ve 
ea a As ata 7 Coser 


In laboratory experiments in gases at pressures of a few millimetres of mercury 
and with fields B less than a few tens of gauss the condition w?<v? holds. 


In this event, with electrons, 


—— eee 
Ww Ee pel d (; \=5n os d (2)=W, 


23m de\y 3m de 
Ee | 2d 2 
Wi3,, ac ao! Cc), Pr) 
cd ae: Batre ne : 
ee ae Te =e PAP) 
tan 0=a-c del! e) fe que”) 
d d d 3 WB 
Sa aE cp =) 2242 \i| oe 
@:e aol! c)-¢ Tree if |< ag? | OE? 
where 
wl Sa 2 7 2d 2 
0=5|¢ rad | Je aol! ec); 
a dimensionless factor. Thus 
CE 
| W |=-, | tan 6 | eae 0 segs Bias SEE ony Sa aeeeen (19) 


Since tan § can be measured directly, W can be calculated if C is known. The 


value of G depends upon the distribution function f(c) and the dependence 
l=Il(c), of | upon e. 


(j) High Frequency Conductivity of a Weakly Ionized Gas 


An alternating electric field H,=X, cos (pt-+«), in the directions +O# can 
be resolved into two oppositely rotating vector fields as follows : 


E,=4X exp (ipt)+4X* exp (—ipt), 
where YX =X, exp (ix) and X*=X, exp (—ia). 
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It follows from equation (15) that the components of W are given by 
W,t+iW,=(We +iW,)+(Wz +iW, ) 
Sees d (03 ees exp (ipt) _X* exp aaa 


6m de/ v—i(w—p)  v—i(w+p) is 


The associated current densities are J,=new, and J, =neW ,. 


When the electric field also possesses components L,=Y cos ( aaa ) and 
H,=Z cos (pt+y) it can be seen that the complex eoerent density J(J iJ J.) 
whose real parts give the physical current density is 


Re x xX* ] 
=r eh Hy sad Racal (ipt)+| o- 1 +] ese—ton, 


J 


in which 
Eels Senha \oxz Oxy 0 | 
pot ayy-oyy > 0) cand) |.or |= opowaj 0+ |, 
O00 sag OpaOMi<a,7 
with X—X, exp (ix), Y=Y, exp (i$), Z=Z, exp (iy), and wa(21) 
Pee eee + gd: e 
Ox = Oyy =10 zy = —1 dy, = (ne?/6m)e—* de eee 
Ae a ee 
de v—i(w+ 7p) 
ne” 


III. DIFFUSION 
(a) General 

The agitational motion of the electrons or ions operates to diminish 
inequalities in their concentration n and to disperse a group of electrons through- 
out the gas. Across an elementary geometrical surface dS at a position where 
grad ” is not zero there is, due to diffusion, a net flux of electrons which is a 
function of the components of gradn. In practice | grad n|/n is small and 
this flux is accurately proportional to —grad n-dS when no magnetic field is 
present. The coefficient of proportionality D is called the coefficient of diffusion. 
The flux is therefore —D grad n-dS=nw-dS, where w is an equivalent con- 
vective velocity that would give the same flux across dS were grad n equal to 
zero. The net transport of electrons across dS in time dt, being the difference 
between those that cross in opposite senses, is therefore 


nw: dSdt=—D grad n-dSdt. 


(b) Formula for the Coefficient of Diffusion D 
Let the direction of grad » be that of the coordinate axis +O so that 
Nw=nw,=—Don/dx, w= —(D/n)en/ox. Let w(e) and n, refer to the electrons 
with speeds between ¢ and e+de. Then w(c)=—(D(e c)/n,)On,| Ox. 
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Consider a volume bounded by surfaces of unit area normal to Ox and a 
distance dv apart. The mean momentum of the electrons (or ions) with speeds Cc 
within this volume is n,mw(c)da and, according to equation (3), the rate at which 
momentum is destroyed within the volume is n,muw(c)g/lda, where 


L=l)/[M/(M +m) —a«]=(4+S8), 


where 8 is defined in equation (5) and g is the mean velocity of an electron (or 
ion) with speed ¢ relative to a molecule. When w(c) is constant this momentum 
is restored by transport of momentum across the boundaries. Let the unit 
boundary surfaces lie at positions # and #+dz on Or. The mean momentum 
transported in time dt by electrons with speeds ¢ in the direction +Oz into the 
volume, across the unit boundary at a, is 4n,mu*dt=jnme*dt. The same 
quantity of momentum leaves the volume in time dt in the direction —Ow and 
by Newton’s third law there is an equal gain of momentum to the volume in the 
direction +Ox. Thus the total gain of momentum is 4n,me*dt in the direction 
+0Ox. Similarly, the gain of momentum in the sense +Oz across the boundary 
at w+-dxw is —4me7dt{n,+(dn,/dx)dz}. 
The total gain of momentum in the sense -++Oz is therefore 


—ime?(dn,/dx)dadt = —(dp,/dx)dadt, 


where p, is the partial pressure of the electrons (or ions) with speeds c. Thus 
the condition that w(c) should not change with time is 
—n,mw(e)g/l—4me7dn,/dx=0, 
or 
—D(c) Lee __ = Wall dn, 
nm, de => Socn, dx? 
from which it follows that 


D(c)=4(le?/g) and D=D(c)=4(le*/g). 
With electrons, g=c and D=Hile). 


When a magnetic field B is present, the free paths of the ions or electrons are 
changed from straight lines between encounters to helices whose axes are parallel 
to B and about which they move with angular velocity w——Be/m. The effect 
of the field is to reduce the coefficient of diffusion in directions normal to B to 
some value D,;<D, whereas that parallel to B retains the value D. It is 
necessary to obtain a formula for D,. 

Let | grad n|=dn/dw, and consider the situation where B is parallel to 
+Oy. The Lorentz force on an electron or ion moving with velocity w across 
a magnetic field B is F=ew x B, consequently the direction of w in the presence 
of a magnetic field B does not in general coincide with that of —grad n. In the 
present instance w has components w, and w, whereas grad n has the single 
component dn/dz. 


The equations of dynamic equilibrium in this case become 


nt —mw (c)g/l+-w,(c)eB} me 


—n,mw,(c)g/l—n w(e)eB —kme?dn,/da=0. 
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Write, as before, o=—eB/m, then 


le? 1 dn 
eres i Ct we = aa Sd 
wc) w,(¢) lg and n,W,(C) 3g [1 + w?l?/g?] dx’ 
whence 

_|e? 1 Mel er 

~ 3g [1+a/g2] 3 (lta) 


D,(¢) 


where T=1l/g=1/v. The coefficient of diffusion is therefore 
D6 tA EN A) Ria ee ee Oe ee (24) 


With electrons g=c and T=l/c=1/v. 
Also, wees tar 
AD em OVE IU tT SILI OL ow ici eae ns a's Sia 5 eee | ategs (25) 


The more general case of diffusion in the presence of a magnetic field with B 
directed along +Oy and grad» arbitrarily directed is expressible in matrix 
notation as follows (#=—Be/m) : 


Ww, D,; 0 wTD,;) (anda 
Mh Wy p=, 0 D 0 On| Oy +. aoeey 2.26) 
w, =—Ee Ds 0 D, ) \@n| ez 


IV. THE DISTRIBUTION FUNCTION f(c) 
Consider first the interchange of energy in a collision between an electron 
(or ion) and a molecule of the gas. 


Fig. 4 


(a) Losses of Energy in Collisions 
Let c and C (Fig. 4) be the velocities respectively of an electron (or ion) 
and a molecule before the encounter and c’ and C’ the velocities following 
the encounter. The velocity G=mc+MC/(m+M)=me'+MC'/(m+M) of the 
centroid X of m and ¥ is unchanged by the encounter but the velocities rand R 


relative to X become r’ and R’. 
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The following relations hold : 


c=G+-4rYy, c’=G+r’, 
whence 
c= G2+r219G-r, c%=—G2+r242G-r’, fp +++ ee (27) 
ce? —¢c’"2= 2G: (r—r’) +7? —1r”?, | 


Further progress requires special assumptions. 

(i) Perfectly Elastic Collisions between Smooth Rigid Spheres.—r?=r? ; all 
directions of r’ about X are equally probable and the mean value of G-r’ is zero. 
Then, 


CC A= 2G-1=2(" or ) a yle—O 
which, after reduction, gives 
go gs A Gira eS aes Ce 
(m+M)? 


Since c:C is zero in the mean, all directions of C being equally probable, it follows 
that the mean loss of energy in an encounter is, when ¢ is given, 


2Mm 
(M +m) 


4me? —4ime?= 


[yme?—4MOF}. oo... (28) 


(ii) Scattering Elastic but not Isotropic.—In this case r’=r, but with c and C 
given the mean value of the projection of r’ on r is not zero but equal to Br 
where § is a function of 7. 


Thus, c?—c’?=2G-:(r—r’)=2(1—8)G-r, and it follows as above, that 


Ge ~<a —B)|me?—M0?—(M —m)c:C] 


so that when C is given all directions and magnitudes 


p —~ 2Mm Ts 
4me*—ime Bares Gramal se gerry Vee al ie Res (29) 


in which 8 is the mean value of the projection of the velocities r’ upon c, which, 
according to Section II (e), is the same as «=cos 0. 


(b) Derivation of the Distribution Function f£(c) 


According to equations (2) and (3) the general form of the distribution 
function with a uniform and constant electric field E is, 


Eel u af 


F (oj) = f(0) eS Ne. ts sa (30) 


In a steady state of motion the number of electrons with speeds between ¢ and 
c+de is 4xF(c,u)e’de and the mean population of this group is constant. The 


MOTIONS OF IONS AND ELECTRONS IN GASES (ail 


mean rate at which energy is supplied to the group is nHewF(c,u) -4nc2de averaged 
over all values of w with ¢ constant and is equal to 


dete Git. 
-- oan act de. 


The group transfers energy to the molecules at the same rate, otherwise its 
population would change. 


Consider first the case in which the molecules are at rest (0?=0) and their 


encounters with electrons (or ions) resemble those between smooth rigid spheres 
(x=0). The rate at which the group loses energy in encounters is 


2Mm g 


AN olla 
(47e7fde)n (+m) 4me i, 


with [,=M1/(M-+m) (Section II (e)). 
It follows that 


Seek das Ve he me ieee 
*\ mg} de 3 de (M+my’ 
whence 

f(e)=Const. exp— reel 7 os (a2) 
where V =Eel/mg. 
In order to infer the form of f(e) when 4MC? is not zero use is made of the fact 
that f(c) reverts to Maxwell’s distribution f(c)=exp (—3mc?/2M0?) when V->0. 


In this condition ~ 
Cid C=— (SMC CAF sw nirsny oie 2: (33) 


and the appropriate form of equation (31), when 4MO? is not zero, may be 
inferred to be 


M aa\f__ _me (34 
(ena itm)" SCS OS (34) 


Thus 


3m e cde 
f(¢)=const. “= arFm| (ree) 


emeriee |. eas 
4(MC?+(M +m)V?) 

¢ mede 
=const. exp = | xf +4(M +m)V™ 


c 
=const. exp — | 


For electrons M+m~M, g->c, V=Eel/me, and the expression for f(c) is 
equivalent to that given by Chapman and Cowling (1952, p. 350). 

The problem of the distribution function f(c) has been considered by many 
investigators and references to their work will be found in the treatises of Chapman 
and Cowling (1952, p. 346) and of Loeb (1955, Ch. IV). 
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Equation (31) may also be derived as follows. The acceleration Ee/m in 
ordinary space corresponds to a constant velocity Ee/m in velocity space with 
a radial component (He/m)(u/e). This radial component produces in the distribu- 
tion given by equation (2) outward flux of representative points over the spherical 
surface with radius c, of amount 


‘He? 1\df 
2 = 5 
ue ee 3) ae ie 


and in a permanent distribution this outward flux is balanced by an inward 
flux brought about by losses of energy in encounters. Let the encounters be 
similar to those between rigid smooth spheres, then the mean loss of speed 
Ac=c—c’' in an encounter is to be obtained from c?—c’?=2Mme?/(M +m)?, 
it being supposed that the molecules are at rest. When Ac/c<1 it follows that 
Ac~Mme|(M-+m)?. The inward flux over the surface with radius ¢ is equal to 
the number of collisions in unit time that are made by all electrons in the velocity 
range ¢ to c+Ac, namely, 


Mme 
(M +m)? 


"fy -47c?- =nft “4707 


me 
(M-+my)’ 
from which it follows that 


sd yen me f 


m g)de M+m’ 


which is equation (31). 

Equation (34) may also be derived with greater rigour as follows. As 
discussed above, the acceleration Ee/m in ordinary space becomes a velocity 
Ee/m of a representative point in velocity space with a radial component Heu/me 
at the surface of a sphere with radius c. This radial component is associated 
with an outward flux of representative points over the sphere equal to 
4nc’nF(c,u)Heu/me where the average is taken with respect to uw with ¢ constant. 
Since 


Pew) =f (ora) 6 = 


it follows that the outward flux of points (in unit time) is 


—A4ne? (He\? 1 df 
—{ n|—] - =. 
3 &) g de 

In a steady distribution this outward flux is cancelled by an equal and 
opposite inward flux that arises from the losses of speed in encounters. First 
suppose that the molecules are at rest. In Figure 5, OA=c is the velocity of 
an electron (mass m) that collides with a molecule at rest (mass M). OX is 
the velocity G=mc/(M-+m) of the centroid of m and M , and G<c. Let the 
velocity of the electron after the collision relative to the centroid be r’=XQ, 
then the total velocity is c°-=OQ. Since r'=XA, it follows that c’<e and all 
representative points c’ lie within the sphere with radius ¢. If the initial velocity 
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is OB with a speed c+, where x—=AB, then the representative points after 
collision lie on a sphere whose centre is at Y ’, where 
m me 
Ox 9p) ~ 
xX Mame +) M-+m’ 
and whose radius r=X’B. This sphere intersects the Sphere with centre O 
and radius ¢ in a circle on which lies the point P. Thus the representative 
points that lie on the spherical cap PBP’ lie outside the sphere with centre O 
whereas those on the complementary arc lie within. Let c2dw be an element 
of surface on which A lies. The number of velocity points within an element of 
volume ¢’dwdx at B is nfc?dwdx and the corresponding number of encounters 


(a) (b) (c) 


Fig. 5 


in unit time is p=c?dwnf(g/l,)dx. Let the angle D‘OP=6, then if the scattering 
is isotropic the number of encounters that give velocity points within the sphere 
of radius ¢ is 4p(1—cos 9) in unit time. When the scattering is not isotropic 
this number becomes 
pe) with a (8) sin 401. 
From the figure, c2=r?+0X’2—20X".r cos 0 with OX’ ={me/(M +m)}(c+2) 
and r={M/(M-+m)}(e+«a), and it can be deduced that 
Mme Mme .. 
Se eceatedss da= — ——_, sin 0d0. 
UH Mm oe 0), da (M +m) 
The number of encounters in unit time for which the initial velocity points lie 
within the element of velocity space ¢?dw.4D and whose final points fall within 


the sphere ¢ is therefore 
2 7 A 
nferdeose | ¢(0)(1 —cos 6)de—" edogau ne | (1 —cos 0)¢@(0) sin 0d0 
OJ A 0 


2,  (M-+m)? 
(1—cos 6), Mme 
i ° (f+m 

= I an te2deo, 


~T(M +m) 


Bue oe 
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where, in accordance with equation (3), l=1)(M-+m) /M(1—cos 6). The 
inward flux of points across unit area in unit time is therefore {gme/l(M +m)}\nf. 
To estimate the influence of the molecular motions, consider Figure 5 (b). When 
the molecules are at rest the velocity of the centroid is OX=mc/(M-+m)=G ; 
G<c. An encounter transforms c into c’=G"r’. When the molecule has a 
velocity C, G=OX becomes G=OX’, c’>c’=G’+r’, with 7r’=r’. As CG 
assumes all possible directions, the end points of ¢” lie on the surface of a sphere 
about P and with radius {M/(M/+m)}C. Thus, in velocity space (Fig. 5 (c)) 
the points which originally lay within an element of volume now become dis- 
persed over a sphere centred on the volume and with radius MC/(M-+m). The 
spherical dispersal of velocity points about such elements of volume lying between 
two spherical surfaces with radii c+ _MC/(M-+m) give a flux of points inwards 
across the spherical surface with radius ¢. If | y|<MOC/(M-+¢m) is the distance 
of an elementary volume from the surface with radius c, then it is readily shown 
that the flux of points inwards across unit surface of the sphere ¢ is 


. MO|(M +m) 

1—cos 6)n- iz df 

: 215 di eee ly p(r-+9gf)au 
(1—cos0)M gMC? df gMC? af 


= 1 ee 


~ 3(M +m), (M-+m)"de 3M —+-m)"ae" 


The mean flux for all values of C is 


ng MC? af 
31 (M +m) de’ 


In a steady distribution of speeds ¢ the total flux across any sphere in velocity 
space is zero, thus, 


= (Be gare 
mn) 3g de 1 (M-+m) 3l(M +m) doa? 
or 
M —)\df me 
ail pe ri he oe 
A(V +n Vee taien? ee (34) 


which is equation (34) whose solution is equation (35). 


The general case in which electrons make both elastic encounters and 
inelastic encounters in which a large proportion of their energy is lost, does not 
appear to lead to a simple formula for the distribution. Neverthalen: when 
electrons move in diatomic gases and their mean energies are not proatce than 
about five times that of the gas molecules, the inelastic losses of energy in 
encounters are associated with changes in the rotational states of the molecules 
Moreover, these transitions are produced by those electrons with energies appreals 
ably greater than the mean energy. Let AQ be the energy absorbed by a molecule 
in changing from one rotational state to a neighbouring state. Suppose that 7 
such collisions AQ/imc?<1. Since c?2 —c?=2AQ/m, AC~AQ/me. Let a pro- 
portion of encounters in which the speed is ¢ and losses AQ occur be ie It 
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follows that the inward flux of points over the sphere of radius ¢ in velocity 
space now becomes, with ,=I1M/(M-+m), 


ress nant (1—a# g/l Ae | Arcetdo—nj| (1 o+(= ae a me Aneede, 


m Mc? | (M +m) 
from which it follows that 
Hel\? df _ eae M-+m)\? xAQ] me 
—3 mg} de m MC?| M ee 
il me 


2 (M aay? 


where a? is the reciprocal of the quantity in the square bracket. Which, when 
modified to allow for the influence of the agitational motions of the molecules 
to give Maxwell’s formula when H=0, becomes 


MC? \d mc 
bic 
whence 
f=const. exp— iat Sat a acsate © (36) 
(M +m)(aV)?+ MC? 
In practice ~<1; m/M<1; 1/a?=1+MxAQ/(mc), 
3mede 
=const. exp — oes eee 37) 
f Po) Mave +e 
Suppose that «(c)=0 when c<c,, then 
c 3mede 
= t. —| — OL amen (a) 
f=const. exp : M(V24.03) 1 ; 8, 
f=const. ex i eziene c>e (d) ii 
eee sieht 0G y) 6. + 


when c,><¢, the distribution, except for the groups with large energies, is the 
same as if all collisions were elastic (a=1, 7=0). 


When ¢ exceeds ¢,, (M/m)xAQ soon exceeds 1 and a—me/(4MAQ)?, 


3mAQ [(¢ «ede 

f-const. exp — (Ee)? aie 
(c) Limiting Form of Equations (35) and (38) 

When V2>C? and c,><¢ the speeds of most electrons are distributed according 


to the law 


3n ff 
f(c)=const. exp — lier bie oh 4 Se it eke eee (40) 


where V =Lel/me. 
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When ([(c)cce™ equation (40) reduces to the form 


4c? big 
ALLO ray regis exp (—c/a")de, ........-- (41) 


where a is a speed and n=4—2r. Also 
ioe) 
cal f(ejede=1. 
0 


It follows from the standard integral 


reo 1 1 
| exp (—yyynay = 0") 


5 n 
that the mean value of the sth power of the speeds ¢ is 


aoe (>) Pin. 1 eee (42) 


n 
The form of the distribution function in a high frequency field is obtained if V? 
in equations (37) and (40) is given the value 


EHe\? alt il 
Visa) Foreeraretcr of 


as follows from equation (14). 


V. Ratio W/D 
From the expressions for the velocity of drift W and the coefficient of 
diffusion D (equations (7) and (22)) it follows that 


W_Be, € (1) | (i ain 
v= ‘aralleal Ppt Beige Sea (43) 


In the special cases in which the speeds e¢ are distributed according to 
Maxwell’s formula (equation (1)) the factor 


d (tc 4 o0 d /Ic? 
= = See od peal 
rae oll es alg) 


Consequently, in this case, equation (43) becomes 


Wes Ke __ He _ EN e ca : ay 
D2 = sip oe ee eee fel (44), 


dme? 
where N, is Avogadro’s number (per mol), x Boltzmann’s constant, and R the 
gas constant. at 


ony 
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The same formula results in the special case in which loc g (inverse fifth power 
law of interaction). 


When ions move in gases their agitational speeds are distributed according 
to Maxwell’s jaw (unless H/p is very large) and equation (44) is satisfied, but the 
motions of electrons in an electric field in a gas do not conform to equation (43) 
in general. The complete expression in this case becomes (with g=c) 


W Ee 
iD (Ga 
with Oe Coes are Hee os (45) 


F=c?.c7 © (le2)/24U6) 


The value of F is determined both by the law of distribution of the speeds ¢ 
and the dependence /=l(c). When loce equation (44) is valid, but when the 
law of distribution assumes the more general form of equation (43) the values*of 
F that correspond to values of n=2, 4, and 6 are F=3/2, 1-312, and 1. 
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GRANULATION NEAR THE EXTREME SOLAR LIMB* 
By R. E. LoucuHeapy and R. J. BRAYt{ 


Rosch (1957) found that on good photographs the photospheric granulation 
remains visible to within less than 10 sec of are from the limb, and sometimes 
to less than 5 sec. Recently, however, these results have been contradicted by 
Edmonds (1960) who, on the basis of an examination of photographs obtained 
in “Project Stratoscope”’ (Schwarzschild 1959), finds that the granulation 
disappears almost completely at 0=75° (33 sec from the limb). Poor focus 
obliterates the last vestiges of the pattern closer to the limb; nevertheless, 
Edmonds concludes that the granulation disappears completely at 6=78° (21 sec 
from the limb). Since a determination of the distance from the limb at which 
the granulation disappears provides an estimate of the height of the top of the 
convection zone (Plaskett 1955 ; de Jager 1959), it is important to resolve this 
contradiction. 

An examination of photographs taken with the 5-in. photoheliograph 
described by Loughhead and Burgess (1958) shows that the granulation remains 
visible very close to the limb, thus confirming Résch’s results. Plate 1 shows 
overlapping regions of the Sun in the neighbourhood of the west limb: both 
enlargements were made from the same original negative, using intermediate 
negatives of slightly different densities in order partially to compensate for limb 
darkening. The white line on Plate 1 (b) indicates the position of the actual 
limb, derived from the original negative. Apart from the much brighter facular 
granules, a number of granules can be seen less than 10 see from the limb ; 
in fact, one rather bright granule can be seen only 4 sec from the limb. Even 
in regions where individual granules are hard to distinguish, the photograph 
gives the impression of a low contrast, foreshortened picture of the ordinary 
1-2 sec granulation. No granules are visible on the original negative in the last 
4sec to the limb; this region does not appear in Plate 1 (b).t However, 
observations of even higher resolution would be required to detect any possible 
continuation of the granulation so close to the limb. § 


These results have been confirmed by other good limb photographs taken 
during the past three years. 


* Manuscript received May 27, 1960. 

} Division of Physics, C.S.I.R.O., University Grounds, Chippendale, N.S.W. 

{ Plate 1 was obtained with the aid of a “ seeing monitor ” (Bray, Loughhead, and Norton 
1959), which is now used to trigger automatically the photoheliograph shutters when the seeing 
signal falls to some predetermined, acceptable value. 


§ On the other hand, we have often obtained photographs showing individual facular granules 
only 1-2 sec from the limb. 


LOUGHHEAD AND BRAY 
A Ae 


GRANULATION NEAR THE EXTREME SOLAR LIMB 


75 
BO 85 80 85 902 


West limb photographed on February 27, 1959, at 11% 44™ K.A.S.T. Both prints were made 
from the same original negative. The white line indicates the position of the actual limb, derived 


from the negative. The upper scale gives the value of 0, the angle between the line-of-sight 
and the normal to the photosphere at the point of observation ; the lower scale gives the distance 
from the limb in seconds of arc. 


Aust. J. Phys., Vol. 13, No. 4 
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Edmonds states not only that the granulation vanishes some considerable 
distance from the limb but that it is replaced by a larger-scale pattern of brightness 
fluctuations with sizes in the range 2-5 sec. Its most prominent features are 
dark areas, which begin to appear at §6=60° and are still easily seen out to 0=86° 
(2 sec from the limb). We have occasionally observed a similar pattern on our 
own limb photographs and had previously made an independent assessment 
of its reality based on an examination of photographs selected from 23 limb films 
obtained during the last three years. However, we concluded that the effect 
is spurious for the following reasons : 


(1) The effect is best seen on mediocre photographs, taken at times of only 
fair seeing. For example, it is particularly prominent on photographs 
affected by a type of seeing sometimes encountered which tends to shear 
the detail by 1 or 2 sec of are—thus destroying the fine detail but leaving 
coarser structures more or less unaltered. On all good photographs 
individual granules can still be distinguished close to the limb ; the dark 
areas no longer seem to dominate the pattern and merely form part of 
the network of intergranular material. 


(2) The value of 0 at which the effect first becomes prominent varies from 
photograph to photograph. 


The effect may be partly due to the fact that on photographs affected by 
poor definition due to seeing or other causes the ordinary granulation is smeared 
out, whereas the larger of the dark areas between the granules (which are 
characteristic features of the ordinary granulation—cf. Bray and Loughhead 
1958, Plate 1) remain visible. Near the limb the crowding together of these 
areas due to foreshortening accentuates the impression of a coarse structure of 
dark areas on a grey background. 


The authors wish to thank Mr. D. G. Norton for help in securing the 
observations and Mr. H. Gillett for processing the films and making the 
enlargements. 
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AN ATTEMPT TO DETECT LINEAR POLARIZATION IN THE 
GALACTIC BACKGROUND RADIATION AT 215 Mc/s* 


By J. L. PawsEey} and EH. HARTINGT 


Current theories of the mode of origin of galactic radio-frequency radiation 
assume the main component to be due to ‘“‘ synchrotron ” emission by relativistic 
electrons in interstellar magnetic fields. Such emission is almost completely 
linearly polarized at the point of origin. The received radiation could, however, be 
substantially depolarized owing to its origin in extended regions of inhomogeneous 
magnetic field, or to effects associated with the rotation of the plane of polarization 
in ionized regions with magnetic fields along the line-of-sight. The detection of 
polarization is a most important observation which could substantiate the 
‘“‘ synchrotron”? emission hypothesis and provide direct evidence on magnetic 
fields in interstellar space. 

Observations show that the received radiation, both from the discrete sources 
and the background, is substantially randomly polarized. Linear polarization 
has been detected at a level of about 5 per cent. at centimetre wavelengths, for 
one only of the discrete sources (the Crab nebula). Careful observations of the 
polarization of the background have been reported by two observers: Razin 
(1958), who used a method based on the rotation of a plane polarized aerial, and 
Thomson (1957), whose method was based on the detection of correlation 
between the signals received on aerials polarized in mutually perpendicular 
planes. Razin worked at wavelengths of 1-45 m (203 Me/s) and 3-3 m (91 Me/s), 
and, in an attempt to eliminate certain sources of error, recorded the difference 
in apparent polarization observed with bandwidths of 200ke/s and 5 Mc/s 
respectively. He reported a difference in apparent polarization on the two 
bandwidths of 2-4 °K at 1-45 m over all parts of the sky observed, except in the 
vicinity of the Milky Way where the effect was zero. He interpreted this result 
as a measure of the polarization of the radiation received on the narrower band, 
that on the wider being assumed substantially zero owing to depolarizing effects 
mentioned above. The 3-3m results were substantially less on a percentage 
basis and were uncertain. Thomson, observing a strip of the sky at declination 
22° N. on 159-5 Me/s with a bandwidth of 4 Mc/s, set an upper limit of 1 per cent. 
to the possible polarization, except in two particular regions where he states, 
‘linear polarization of the order of 1 per cent. may have been detected”. The 
brightness temperature was of the order of 500 °K in this case. 


Razin’s results, if correct, are of major importance since they imply the 
existence of magnetic fields and high energy electrons widely distributed through- 
out the galactic corona. Independent confirmation would be most valuable. 


* Manuscript received September 7, 1960. 
} Division of Radiophysics, C.S.1.R.0., University Grounds, Chippendale, N.S.W. 


SHORT COMMUNICATIONS TAL 


It is unfortunate that Thomson’s negative results for most of the area examined. 
cannot be directly compared with Razin’s since Thomson used a wide bandwidth 
for which Razin claimed that the observed polarization should be inappreciable. 
The present paper describes an attempt to confirm Razin’s observations. 


Observations were made on a frequency of 215 Mc/s with bandwidths of 
1 Me/s, and 700 and 300 ke/s, using a Dicke-type receiver connected to a fixed 
60 ft paraboloid directed towards the zenith. A linearly polarized feed aerial 
at the focus was rotated back and forth through two revolutions and the resulting 
changes in output recorded. 

Spurious results can arise in several ways. If the directional diagram is not 
circularly symmetric about the axis of the aerial, output fluctuations may arise 
through directional effects if there are intense sources off the axis of the aerial. 
We attempted to minimize this effect (1) by using for the feed aerial a 2-dipole 
array having nearly identical directional patterns in the E and H planes and 
(2) by restricting observations to times when the principal intense sources, the 
Sun and the centre of the Galaxy, were below the horizon. 


Output changes may also arise through systematic aerial impedance changes 
associated with the rotation of the system. Such impedance changes were 
observed and corrective measures employed, but apparently a residual effect 
remained. A periodic fluctuation of output with rotation, simulating a polar- 
ization effect of magnitude about 1-5 °K, was observed, which remained constant 
in phase and amplitude over the whole of the survey. A genuine polarization 
effect would be expected to vary with position in the sky, as well as with time 
owing to changing Faraday rotation in the ionosphere, and we feel confident 
that this effect was spurious. We have therefore subtracted the mean fluctuation 
from the individual observations. 

A series of observations was made in the period December 1959 to February 
1960. The beamwidth of the aerial was 7° and the region observed was a strip 
centred on declination 34° S. and ranging in right ascension from 1" 30™ to 10°. 
The noise figure of the receiver was 3-6 (5-5 dB). We consider we could detect 
reliably a periodic output change ranging from 1-5 °K at a bandwidth of 300 ke/s 
to 1°K at 1 Mc/s. The background temperature over this region in the sky is 
about 150 °K. No output changes which could be attributed to polarization 
were observed. 

These observations set an upper limit of 1 per cent. to the degree of linear 
polarization over the strip of sky concerned at a frequency of 215 Me/s and for 
bandwidths of 1Me/s, and 700 and 300ke/s. This result conforms with 
-Thomson’s negative result for most of the region he examined. It appears to 
contradict Razin’s claim that most of the sky remote from the Milky Way shows 
polarization, but it must be pointed out that Razin’s and our actual observations 
refer to different regions in the sky, Razin’s to a strip about declination 56° N., 
ours to one about 34°S. 

A detailed account of the observations is available in a laboratory report 


(Harting 1960). 
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A SYSTEM FOR RECORDING AND INTEGRATING PHYSICAL 
; MEASUREMENTS* 


By N. E. RIDERT 


Introduction 

Many occasions occur when it is required to make recordings of one or more 
fluctuating physical quantities. When the measurements to be made, either of 
the same quantity at a number of positions or of various quantities at one position, 
exceed quite a small number, continuous recording is not a practical proposition. 
In any event it may not be necessary. The frequency with which it would be 
desirable to record a particular measurement will depend on: (a) the expected 
rate of fluctuation of the quantity concerned having regard to smoothing that 
may be introduced by the primary sensing device and the recorder in use ; (b) the 
magnitude of the expected fluctuations, again almost certainly modified by the 
apparatus; (c) the resolution available in a single measurement; (d) the 
accuracy to which a time mean value is required together with the time over 
which the mean is to be taken. 

Very often the frequency with which individual readings can be made is 
dictated by the number of different observations required together with the 
characteristics of the available recording system. Recordings are normally 
made with appropriate sensing heads working in conjunction with multichannel 
recorders. The record finally obtained is a representation of the spot values of 
the quantities at time intervals depending on the switching cycle of the particular 
recorder. Effective multiplication of the number of recording channels with a 
decrease in the frequency of recording of a single variable can sometimes be 
achieved by external switching which is synchronized to any internal switching 
in the recorder. Various devices have been reported (e.g. McHugo 1959) which 
will provide an integrated value of at least one variable which is working into a 
single-channel recorder. A limit to the number of different channels that may be 
used on one recorder is often set by the chart record itself as this becomes too 
confused for convenient analysis. In all cases the chart has to be subsequently 
measured before the desired numerical values are available. This is often tedious 
and time consuming and provides opportunity for subjective error. The system 


* Manuscript received May 9, 1960. 
{ Division of Plant Industry, C.S.I.R.O., Canberra. 
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to be described is fully automatic and no error of this type is possible, the 
individual measurements and their totals over any desired period being presented 
on counters. Moreover, the initial signals themselves need only be of the order 
of a fraction of a microvolt and, since reflecting galvanometers provide the initial 
signal sensing elements, the system is very flexible. It may be made to cater 
for a wide range of signal magnitudes. To make what is to follow specific the 
system as constructed to record and integrate the dry- and wet-bulb temperatures 
from 20 positions will be described. Modifications for other applications will 
readily suggest themselves. 


The System in Outline 

The basic idea was first developed by House, Rider, and Tugwell (1960) 
for use in their surface energy balance computer in which it was necessary to 
digitize the deflections of a reflecting galvanometer. It has not been described 
previously and the present system is a further development to more general use. 
The usual semitransparent scale of a reflecting galvanometer is replaced by a 
long photovoltaic selenium cell which is covered with a grid consisting of alternate 
opaque and transparent bands of equal width. The width of the image of the 
projection lamp slit at the scale distance is made the same as that of these bands. 
As the image of the slit moves across the photocell scale a series of pulses is pro- 
duced by the cell. The number of pulses is a measure of the distance swept out 
by the slit image. If these pulses are applied to an electromagnetic counter via — 
the necessary pulse amplifier a record of the change in galvanometer deflection is 
obtained. When a current passes through a galvanometer it takes up a position 
defined by the magnitude of the current and its closed circuit zero position. 
If the galvanometer is disconnected from the circuit in which this current 
originates and connected to another circuit, the bias circuit, from which a small 
current flows which is sufficient to move the image of the projection lamp slit off 
one end of the scale, the number of pulses produced will be a measure of this 
position with respect to the end of the scale concerned. Repetition of the 
operation with the galvanometer connected to its critical damping resistance 
in place of the signal circuit yields a measure of the zero position of the galvano- 
meter with respect to the same end of the scale. This position is recorded 
on a second counter. The difference in the counter readings is then a measure 
of the original current which it is desired to measure and record. Suitable 
switching enables a sequence of signals to be applied to the one galvanometer 
and the resulting pulses to be passed to individual counters. The accuracy of 
recording will be determined by the galvanometer v. signal characteristics and 
the size of the grid spacing used. If the galvanometer and its associated circuits 
are allowed to pass through a number of cycles of operation the integrated 
deflections are obtained on the counters, which can then be photographed, or 
the indications can be printed out, at desired intervals. 


Application to Temperature and Humidity Recording 

In our present research programme we wished to record a two-dimensional 
grid of temperature and humidity. Absolute accuracies were relatively 
unimportant, but temperature and humidity differences in the vertical and 
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horizontal were required to be as accurate as possible consistent with the use of 
a reasonable amount of apparatus. Fine wire thermocouples of copper/copper- 
nickel were used in psychrometer units mounted on masts, there being five such 
units on each of four masts. Since the dry- and wet-bulb circuits are identical 
it is only necessary to consider one set, the addition of the other set merely 
requiring a duplication of the facilities to be described. The arrangement of 
the masts and wiring (dry bulbs only) for masts 1 and 2 is shown in Figure 1. 


MAST NO. 


HEIGHT 1 2 


COPPER Ls 


__ COPPER 
~ NICKEL 


MULTICORE CABLE 


Fig. 1—Schematic diagram of dry-bulb thermocouple circuits 
for two masts. 


On mast No. 1 it will be seen that, as well as the five thermojunctions at heights 
from A to # inclusive, a standard junction (S,) and link junction (L,) are provided 
and that all the copper-nickel leads are connected together in one point. Masts 
Nos. 2-4 are identical and differ from the first in that they have no standard 
junctions. All the link junctions and the one standard junction are maintained 
in melting ice. Selection of any pair of copper leads from one mast provides 
an e.m.f. which is a measure of the small temperature difference existing between 
the corresponding junctions. These dry-bulb circuits use five galvanometers 
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and switching arrangements, to be discussed later, enabling connections to be 
made according to the programme given in Table 1. The time interval of 7 sec 
is allowed in each circuit since the galvanometers used (Tinsley type 4500, 2 sec 
periodic time, resistance and Sensitivity approximately 50Q and 10 cm/pA 
at 33 cm scale distance respectively) required this time to take up full deflection 
with the fairly low resistance circuits employed, themselves necessary to obtain 
the required deflection per degree of temperature difference between junctions. 
The times not accounted for in Table 1, i.e. 8 to 11, 20 to 23 sec, etc., were occupied 
in sweeping the slits of the galvanometers off the ends of the scales to obtain 
records on the deflections at 7, 19, 31, 43, and 55 sec. With this schedule a 
Spot value of the temperatures at each position of exposure is obtained once per 
minute. A faster schedule would be possible if (a) shorter period galvanometers, 
or (6) more galvanometers, or (c) higher external circuit resistances were used. 


TABLE 1 
SCHEDULE OF GALVANOMETER OPERATIONS FOR DRY BULB CIRCUITS 


: Galvanometer No. 
Time 

(sec) 

1 Z 3 4 5 

0-7 L', @ATs fas), fie SEL Ds fle T1,¥. Ty, ZeETO 
12-19 1 ve a be Gl Dy Abe Phe De Ohh 4B Be Oki T,v. Sy 
24-31 AS ie dk heh Bf Dh On Ahoy, Io Oe hoe US OS A 
36-43 LIER IEE TT Veeeao DRE Li he DBs Os Whe Ts OE ss 
48-55 Zero Zero Zero Zero | ; SMA Dy, Uhre 


However, the present arrangement compares favourably with what could be 
achieved with a reasonable number of multichannel recorders, both from the 
point of view of the resolution available in a single record and in the frequency 
of measurement at one point. It will be noted that one operation of each 
galvanometer is occupied in recording its zero position and that the total of 25 
operations requires the use of the same number of electromagnetic counters. 
The link junctions are required in order to avoid the use of long lengths of copper- 
nickel wire between masts which would otherwise be necessary.so that galvano- 
meter No. 5 could execute its last three operations in the schedule. Since the 
accuracy required in the 7, v. 7, T, v. T,, and T, v. T,, temperature difference 
is the same as in the 7, v. T,, 7, v. 7, etc. determinations the additional resistance 
that would be included in the former circuits by long lengths of copper-nickel 
wire would not be admissible. The link junctions would not be necessary if all 
- measuring points were in close proximity. The resistances which have been 
shown in Figure 1 are provided to allow adjustment of the temperature difference 
». galvanometer deflection characteristics of the various galvanometers and 
circuits. Actually they are not mounted at the masts but ab the recording 
apparatus. These resistances are adjusted so that the sensitivity of all circuits 
is the same, with the exception of the 7, v. 8, circuit which has a sensitivity of 


about one-fifth of the other circuits. 
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Details of Circuit Arrangements 

The basic timing device which controlled the whole cycle of operations 
was a geared down synchronous motor which operated a 1 sec cam, this cam in 
turn operating a mercury-in-glass switch which energized the coil of a miniature 
uniselector. The uniselector had 30 output contacts per level, there being seven 
levels with provision for inter-level switching. This device provided the means 
whereby relays could be energized in time intervals down to 1 sec in the cycle of 
operations lasting 1min. By grouping the uniselector output contacts any 
relay could be made to remain energized for any desired period during the cycle. 


AMPLIFIER 
NO.1 b 


ee 


COUNTERS H> ve K 
A 


INPUTS T3 
FROM 
MAST + 


NO 1 3 
rs 


BIAS 
CIRCUIT 


Fig. 2.—Switching circuits for galvanometer and amplifier No. 1. The contacts are 
shown in their normal positions. The times when they are not in the positions shown 
are given in Table 2. 


We require that the switching arrangements should connect the correct 
circuits to the galvanometers at predetermined times, apply the currents required 
to generate the galvanometer sweeps, and select the correct electromagnetic 
counters to receive the outputs of the pulse amplifiers. Several other functions 
are also required and were arranged quite easily. The shutter of the automatic 
35 mm camera used to record the counter indications must be operated at the 
correct time in the cycle and after the desired number of cycles have been 
completed, and the field illumination must be similarly controlled. If, as in 
the present case, the counters are fitted with reset coils these must be energized 
immediately following the film exposure. To cater for these operations a second 
uniselector of the same type as the first is driven from the latter in such a way 
that it moves on one step per half minute. Many arrangements are possible 
with these uniselectors according to the requirements of the particular application. 
In our case an integrated total was required after 20 cycles of operation and 
arrangements were made for this. We also wished to obtain on occasions an 
individual record of each cycle of operation and change over to this type of 
operation was provided for. Further, times occurred when, owing to some 
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emergency, it was necessary to stop a run and restart. This facility was incor- 
porated by the use of homing arcs on the uniselectors. So many arrangements 
are possible that it is thought not to be worth while to detail the uniselector 
and relay energizing wiring and we shall confine ourselves to giving some account 
of the type of switching that is used in the circuit of one galvanometer (G,) which 
is illustrated in Figure 2. Referring to this figure it will be seen that the copper 
lead from the T, junction is connected to one side of G, and that the other side of 


TABLE 2 
TIMES OF OPERATION OF RELAY CONTACTS SHOWN IN FIGURE 2 


Contact Time of Operation (sec) 
Mis ies, Se 7-11, 19-23, 31-35, 43-47, 55-59 
Dae 0-11 
oo 12-23 
S Pre 24-35 
i, _K, 36-47 
M,, M, 48-59 
U 8-11, 20-23, 32-35, 44-47, 56-59 


the galvanometer is connected via the change-over contact S, to one side of a 
series of relay contacts marked F,, H,, J,, K,, and M,, the other sides being 
connected directly to T,, T,, T,, T;, and R,,, the galvanometer critical damping 
resistance, respectively. The times during which these and all other contacts 
are not in the positions shown in Figure 2 are given in Table 2. We see, for 
example that contact F, is closed in the period 0-11 sec so that G, is in the T, v. T, 
circuit until S, changes over at 7 sec, thus connecting the galvanometer into its 


SK 10KQ 


Fig. 3.—The galvanometer bias circuit. 


bias circuit which produces the sweep of the slit across the cell. During the 
period up to 7 sec contact S, is open so that although F, is closed and the 7, 0. ar. 
counter is in circuit no pulses which may be produced by the cell and amplifier 
due to fluctuations in the deflection of the galvanometer will be passed to the 
counter, only those arising in an actual sweep being recorded. Similar switching 
circuits are provided for each galvanometer and its associated amplifier and 
counters and each has its own bias circuit, illustrated in Figure 3, which should 
be considered in conjunction with Figure 2. The bias circuit consists of a 1-5 V 


748 SHORT COMMUNICATIONS 


dry cell and a delay line which is designed to limit the rate at which the bias 
current rises through the galvanometer, thus ensuring that the sweep commences 
slowly and later speeds up. On first making contact S, all the current flows 
into the condenser and the current rises through the galvanometer and its parallel 
resistances as the condenser charges. This arrangement by itself was found 
to be insufficient to exercise enough control on the sweep speed and so two 
parallel resistances were used as shown, having values of around 3 and 8 Q. 
At first these are both in parallel across the galvanometer but 1 sec after the start 
of a sweep contact U opens to increase the parallel resistance from about 2 to 
8 Q. The exact values of these two small resistances are adjusted to suit the 
particular galvanometer. All these bias circuit arrangements were necessary 
since the photovoltaic cells appear to have an initial lag to changes in light level 
but once they commence to operate they may be pulsed satisfactorily at speeds 
up to afew ke/s. The highest speed here is governed by the maximum counting 
rate of the electromagnetic counters used. 


SIEMENS 


COUNTER 
1542 


@ © | 
x,0Cl6, 


@ ©.— + 24v 


— 24V 


7 
47K2 S47K2Q 


Fig. 4.—The pulse amplifier circuit. 


One pulse amplifier was used with each galvanometer. Although photo- 
voltaic cells do not lend themselves very readily to incorporation in electronic 
circuits (owing to their low internal impedance and small output at the low light 
levels used here) a satisfactory amplifier design was finally produced and is 
shown in Figure 4. Very briefly it consists of two EF86’s connected as a long 
tailed pair, followed by a double triode (12AX7) to give further amplification, 
and the last stage is a Schmitt trigger circuit (Schmitt 1938) which operates a 
high speed relay. A voltage gain approaching 104 is obtained. RV, is the gain 
control and the relay firing point is set with RV,. The frequency range over 
which the amplifier will operate is from less than 1 c/s up to a maximum deter- 
mined by the relay used, in our case about 80 c/s. However, the Sweep speed of 
the galvanometer slit must, as has already been indicated, be kept below the 
maximum working speed of the counters which in our case was around 
35 counts/sec. Precautions have to be taken in the amplifier against oscillation 
at mains frequency. The valve heaters are connected in series and supplied 
with 24 V d.c., this supply being also required for the uniselectors, relay coils, 
and projection jae The 300 V h.t. supply must be taken from a good quality 
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power pack. The contacts of the high speed relay are conveniently protected 
by the power transistor which is shown in Figure 4. Here, with 24 V counters, 
the current passed by the relay contacts is completely non-inductive and is 
less than 10 mA, while virtually the full 24 V is dropped across the counter coil 
when the relay makes. 

The photocells used had the dimensions 15 by 2-5 cm and the metal grids 
passed five pulses per 1 cm of light movement. Since the cells tend to have 
small weak areas it is convenient to replace the single band of light from the 
projector lamp by a double slit or small grid of light the image of which has the 
Same spacing as the grid covering the cell at the scale distance. 

It will be clear that the records obtained are independent of long-term 
galvanometer zero drift. It is only necessary for the zeros to remain fixed for 
the period occupied by one cycle of operation to eradicate any error which might 
arise on this account. 


Conclusion 

The particular form of the system described has been in use in a mobile 
laboratory for some time where it has been subjected to unusually rough treatment 
in travelling. All faults that have occurred can be traced to damage in transit 
and as a static installation the system should prove as reliable as one based on 
conventional recorders. If the latter had been used it is estimated that several 
days’ work would be involved in chart analysis before a grid of temperature and 
humidity could be constructed for a single 30min period. Here the basic 
information is immediately available in numerical form and such work as is 
necessary is limited to the use of the psychrometric equation to arrive at values 
of vapour pressure. It is not necessary to list other possible uses for the system— 
any application to which photographic, multichannel potentiometric or thread 
recorders, or pen recorders have hitherto been appropriate could be suitably 
served. The cost of the system is less than that of potentiometric recorders 
which would provide the same facilities. 
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NOISE SUPPRESSION IN PULSE RECEIVERS* 
By E. ©. McLAUCHLANT 


In the course of measuring Southern Hemisphere meteor rates by radar 
methods, difficulty has been experienced in combating the effect on the meteor 
rate of a variable background, due to noise of both solar and man-made origin. 

The equipment, operating at 69-5 Mc/s, radiates pulses of 23 usec duration 
at a pulse recurrence frequency of 150 per sec, and a peak power of 90-100 kW. 
The echo is detected by a very low noise receiver and applied to the grid of a 
cathode-ray tube, the trace of which is completely blacked out in the absence of 
signal. Echoes are recorded on 35 mm film which moves continuously at the 
rate of 1 ft/hr past the face of the tube. Alternate sweeps of the time-base are 
shifted slightly in order to present echoes with more certainty as double dots. 

A serious shortcoming of this fairly standard system is that any appreciable 
rise in background noise results in excess darkening of the film, and hence lowered 
recognition of echoes. This effect is somewhat aggravated by the inherent 
short grid base of the average cathode-ray tube. The incidence of total black-out 
due to man-made interference has been greatly reduced by the inclusion of a 
compression amplifier in the video section of the receiver. 


Receiver Modifications 

The compression amplifier extends the principle of a constant volume 
audio amplifier described by G. J. Pope (1952). 

The essential control element is a cathode follower used as an anode load of a 
pentode voltage amplifier. Referring to Figure 1, V1 is arranged to give sufficient 
amplitude of video signal from the receiver, and to be of such a phase that the 
final video amplifiers V8 and V9 can be operated with positive-going signals 
from a near cut-off point. The gain control R1 is adjusted so that V2 will not 
be overloaded by the strongest video input signals. 

The resistive anode load of V2 is shunted by the cathode-anode impedance 
of the cathode follower V5 biased towards cut-off. As the signal applied to V2 
is essentially positive-going, V2 is also biased towards cut-off. The anode of 
V2 feeds the voltage amplifier V3, whose output is rectified by the diode V4 to 
provide the positive voltage necessary to offset the standing negative bias on V5. 

The restricted bandwidth of the V3 circuit, together with the very long 
time-constant of the rectifier (V4) load circuit, serves to make the bias on V5— 
and hence the gain of V2—a function of the steady background noise, while it is 
unaffected by the signal pulses. The tendency of the cathode choke (L1) to 
oscillate on the receipt of very strong pulses was adequately suppressed by the 


* Manuscript received May 30, 1960. 
} University of Canterbury, Christchurch, New Zealand. 
{ Porr, G. J. (1952).—Electron. Engng. 24: 464. 
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diode V6. The stage V7 acts as a conventional screen-controlled limiter, the 
limit level being adjusted to avoid defocusing of the cathode-ray tube under 
strong signal conditions. 


System Efficiency 

Normal conditions were synthesized for the purpose of assessing the 
advantage of this additional amplifier over a conventional video system. The 
_ output of a diode noise generator was injected into the receiver simultaneously 
with artificial echoes, of 0-2 sec duration, having the same pulse recurrence 
frequency and pulse width as the operational transmitter. Both the echo 
simulator and the noise generator were of equal output impedance. 

Five echoes decreasing in amplitude from 4 to 1 uV and spaced 30 sec 
apart, at constant range, were recorded during each setting of the noise generator 
output. As the test was intended to be comparative rather than absolute, no 
attempt was made to match the similar impedance sources to the receiver input. 
(It is probable, therefore, that the minimum detectable signal power for actual 
echoes will be below the figures obtained in the present test.) The range of 
available noise power from the noise generator was less than that expected from 
variations of aerial noise temperature. Consequently, it was necessary to increase 
the receiver gain part way through the test. The results of the tests are shown 
in Plate 1, first for the normal video amplifier, and then with the compression 
amplifier added. The conditions existing at each of 11 consecutive 3-min intervals 
were as follows : 


Step 1— The noise generator was at zero output, and the receiver gain 
was adjusted to give the same detector current as when it was 
connected to the aerial under quiet conditions. Under these 
conditions the total noise power is essentially that generated by 
the receiver itself, which is 3:94 x10-16 W. 

Steps 2- 6.—The noise generator output was increased in five steps of 
326 101° We 

Step 7.— The noise generator output was returned to zero and the receiver 
gain increased to give the same receiver detector current as had been 
obtained in step 6. The noise generator output was then increased 
Oo 6 LOS eave 

Steps 8-11.—The noise generator was increased in further steps as for the stages 
3-6. 

Thus by adding the noise power derived from the noise generator to that generated 

by the receiver the total noise increased from 3-94 10-16 W for step 1 to a 

final value of 39-94 x10-1° W for step 11. The plate shows clearly the increased 

readability of echo rates during periods of high noise, while those obtained during 
periods of low noise are unaffected. 


This work was carried out as part of Contract AF64(500)-6 from the U.S. Air 
Force Air Research and Development Command. The author is indebted to 
Dr. ©. D. Ellyett for his assistance in preparing the paper. 
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